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PRECEDING PAGE BLANK NOT FEMES 

FOREWORD 

STEAP II is a series of three computer programs developed by the 
Martin-Marietta Corporation for the mathematical analysis of the navigation 
and guidance of lunar and interplanetary trajectories. STEAP is an acronym 
for Space Trajectory Error Analysis Programs. The first series of programs 
under this name was developed under contract NAS 1-8745 for Langley Research 
Center and was documented in two volumes (STEAP User's Manual, STEAP Analytical 
Manual) as NASA Contract Report 66818. Under contract NAS 5-11795 the STEAP 
series was extensively modified and expanded for Goddard Space Flight Center. 
This second generation series of programs is referred to as STEAP II. 

STEAP II is composed of three independent yet related programs: NOMNAL, 
ERRAN, and SIMUL. All three programs require the integration of n-body 
trajectories for both interplanetary and lunar missions. The virtual mass 
technique is the scheme used for this purpose in all three programs. 

The first program named NOMNAL is responsible for the generation of 
n-body nominal trajectories (either lunar or interplanetary) performing a 
number of deterministic guidance events. These events include initial or 
injection targeting, midcourse retargeting, and orbit insertion. A variety 
of target parameters are available for the targeting events. The actual 
targeting is done iteratively either by a modified Newton- Raphs on algorithm 
or by a steepest descent-conjugate gradient scheme. Planar and nonplanar 
strategies are available for the orbit insertion computation. All maneuvers 
may be executed either by a simple impulsive model or by a pulsing sequence 
model . 

ERRAN, the second program of STEAP II, is used to conduct linear error 
analysis studies along specific targeted trajectories. The targeted trajectory 
may however be altered during flight by retargeting events (computed either 
by linear or nonlinear guidance) and by an orbit insertion event. Knowledge 
and control covariances are propagated along the trajectory through a series 
of measurements and guidance events in a totally integrated fashion. The 
knowledge covariance is processed through measurements using an optimal 
Kalman- Schmidt filter with arbitrary solve- for/consider augmentation. 

Execution errors at guidance events may be modeled either by an impulsive 
approximation or by a pulsing sequence model. The resulting knowledge and 
control covariances may be analysed by the program at various events to 
determine statistical data including probabilistic midcourse correction sizing 
and effectiveness, probability of impact, and biased aimpoint requirements. 

The third and final program in the STEAP II series is the simulation 
program SIMUL. SIMUL is responsible for the testing of the mathematical 
models used in the navigation and guidance process. An "actual" dynamic 
model is used to propagate an "actual" trejectory. Noisy measurements from 
this "actual" trajectory are then sent to the estimation algorithm. Here 
the actual measurement, the statistics associated with that measurement, 
and an "assumed" dynamical model are blended together to generate the filter 
estimate of the trajectory state. This process is repeated continually 
through the measurement schedule. At guidance events corrections are computed 
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based on the estimate of the current state. These corrections are then 
corrupted by execution errors and added to the "actual" trajectory. The 
statistics and augmentation of the filter, the mismatches in the "actual" 
and "assumed" dynamics, and the execution errors and measurement biases 
may then be varied to determine the effects of these parameters on the 
navigation and guidance process. 

The documentation for STEAP II consists of three volumes: the Analytic, 
Programmer's and User's Manuals. Each of these documents is self-contained. 

The Analytic Manual consists of two major divisions. The first section 
provides a unified treatment of the mathematical analysis of the STEAP II 
programs. The general problem description, formulation, and solution ate 
given in a tutorial manner. The second section of this report supplies the 
detailed analysis of those subroutines of STEAP II dealing with technical 
tasks . 

The Programmer's Manual provides the reader with the information he 
needs to effectively modify the programs. Both the overall structure of the 
programs as well as the computational flow and analysis of the individual 
subroutines is described in this manual. 

The User's Manual contains the information necessary to operate the 
programs. The input and output quantities of the programs are described 
in detail. Example cases are also given and discussed. 
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1 • INTRODUCTION 


The User's Manual provides the user of the STEAP II programs with all 
the information necessary to input these programs and interpret the output. 

Chapter 2 presents a summary of the four programs constituting STEAP II 
These are the trajectory propagation package, the targeting program NOMNAL, 
the error analysis program ERRAN, and the simulation program SIMUL. 

Chapters 3 and 4 describe the detailed input and output, respectively, 
for the three programs NOMNAL, ERRAN, and SIMUL. Restrictions on the input 
procedure for these programs are also presented. 

In Chapter 5 are discussed actual sample cases which were run using 
the STEAP II programs. These sample cases are presented primarily to demon- 
strate the operation and versatility of NOMNAL, ERRAN, and SIMUL and to 
assist the user in the input/output procedure for these programs. Selected 
pages from these sample cases are found in the Appendix of the User's Manual 
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SUMMARY OF MODES 


The Space Trajectory Error Analysis Programs (STEAP) consist of 
four subprograms or operational modes. The first mode, used as a 
subroutine by each of the other three programs, is the trajectory mode 
VMP by which an n-body trajectory (lunar or interplanetary) is propa- 
gated by the virtual mass technique. The second mode is the nominal 
trajectory generator or targeter (NOMNAL) by which a lunar or inter- 
planetary trajectory meeting specified conditions is determined. The 
third mode is the error analysis program ERRAN in which the navigation 
and guidance characteristics of a nominal trajectory are analyzed by 
linearly propagating knowledge and control covariances along the 
trajectory. Finally the simulation mode SIMDL tests the mathematical 
models used in the navigation and guidance processes by modeling the 
tracking and correction of an "actual" trajectory. In this chapter 
a general description of each of these modes will be provided. 

2. 1 The Virtual Mass Propagator VMP 

The dynamic model used by STEAP is supplied by the trajectory 
propagation package. The only external forces acting upon the space- 
craft are assumed to be the gravitational forces of the celestial 
bodies considered in the integration. Both the spacecraft and the 
gravitational bodies are assumed to be point masses so neither space- 
craft attitude nor planet asphericities are considered. 

The celestial bodies to be in the integration are specified by 
the user and may include the sun, any of the nine planets, and the 
earth* s moon. The motion of the planets about the sun and the moon 
about the earth are modeled by using mean ecliptic elements of date. 

If the user desires, each of the planets can be set in a fixed ellipse 
referenced to some epoch for speedier computation. 

The coordinate system used in the integration is also specified 
by the user. The options available are either heliocentric ecliptic 
or barycentric ecliptic (nominally for lunar trajectories). 

The actual scheme used in the propagation of the trajectory in 
the virtual mass or varicentric technique (see reference 15) . No 
actual integration is performed by the trajectory mode; the key idea 
of the virtual mass technique is to build up an n-body trajectory by 
using a sequence of conic sections around a moving effective force 
center called the virtual mass. At each instantaneous moment along 
the trajectory, the combined effects of all the gravitational bodies 
can be viewed as resulting from a fictitious body of unique magnitude 
and position which is called the virtual mass. The computational pro- 



cedure then assumes that over a small time interval, the motion 
of the spacecraft can be represented by a two-body conic section 
arc relative to this virtual mass. The complete trajectory is thus 
generated by a series of small arcs pieced together in steps while 
updating the position and magnitude of the effective force center. 
The main advantage of the virtual-mass technique is that the tedi- 
ous numerical integration of the differential equations is avoided. 

Another significant feature of the virtual-mass technique is 
its flexibility. By varying a simple parameter called the "accur- 
acy level" related to the true anomaly increment of each step, tra- 
jectories ranging from a sequence of relatively few conic section 
arcs corresponding to a very approximate solution to those requir- 
ing a large number of arcs corresponding to highly accurate solu- 
tions may be generated. 

2.2 The Nominal Trajectory Targe ter NOMNAL 

NOMNAL is responsible for the generation of a nominal trajec- 
tory for either lunar or interplanetary missions. The method or 
propagation in either case is the virtual-mass n-body integrator. 
The trajectory may be processed through a series of deterministic 
maneuvers including initial or injection targeting, subsequent 
retargeting, miniprobe targeting, and orbit insertion. A variety 
of target parameters are available for the targeting events. Both 
coplanar and nonplanar strategies are permitted in the orbit in- 
sertion maneuver. 

If an initial state for the problem is known, this may be read 
in to start the trajectory. Otherwise NOMNAL generates its own 
zero iterate. In interplanetary missions this involves solving 
the Lambert time-of-f light equation for the massless planet tra- 
jectory that connects the desired initial and final positions in 
the specified time interval. Four options are available in de- 
scribing these reference points. 


Initial Point 

Final Point 

Launch Planet 
Launch Planet 
Specified Point 
Specified Point 

Target planet 
Specified Point 
Target Planet 
Specified Point 


If the initial point is referenced to the launch planet, a launch 
profile is consulted to generate a realistic set of injection con- 
ditions consistent with the heliocentric trajectory. 




For lunar trajectories a slightly different procedure is used. 
The required data for the lunar zero iterate includes specification 
of the desired semimajor axis with respect to the moon, radius and 
time of closest approach to the moon, and inclination to the lunar 
equator. Then the zero iterate is generated by first targeting a 
patched conic trajectory and then a multiconic trajectory to the 
desired conditions. 

A targeting event may be processed immediately after obtaining 
a zero iterate state or at any point along the nominal trajectory. 
At a targeting event the current velocity is refined to yield a 
trajectory satisfying target parameter constraints. The possible 
target parameter are: 


1) 

TPS 

5) 

B • T 

9) 

SMA (Lunar) 

13) 

DCP 

2) 

TSI 

6) 

B • R 

10) 

XF 

14) 

RAP 

3) 

TCS 

7) 

RCA 

ID 

YF 

15) 

TPR 

4) 

TCA 

8) 

INC 

12) 

ZF 




The targeting method to be used is specified by the user. Either 
a modified Newton-Raphson algorithm or a steepest descent/conju- 
gate gradient technique may be used. 

Orbit insertion events are also available in NOMNAL. At a 
specified time the spacecraft state relative to the target body 
is computed. The resulting conic trajectory relative to the tar- 
get body is then compared with the desired orbit to determine the 
optimal time to make the insertion and the required correction. 

At the proper time the velocity correction is then implemented. 

Two strategies are permitted in the orbit insertion computation: 

1) Coplanar - The desired semimajor axis, eccentricity, and 

periapsis shift of a coplanar orbit are speci- 
fied; 

2) Nonplanar - The desired plane of the postinsertion state 

is specified along with nominal values of the 
orbit elements. 

The targeted correction, orbit insertion correction, or an ex- 
ternally supplied correction may be executed if desired. Two mod- 
els are available for this implementation — a simple impulsive ad- 
dition or a more complex multiple pulse model. 
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NOMNAL is also capable of targeting a set of three miniprobes 
to three specified target sites. Since achieving impacts at three 
specified points on the planet surface constitutes a six-degree- 
of-freedom constraint while only four miniprobe release controls 
are available, any targeting process can, at most, achieve a mini- 
mum-miss solution. NOMNAL uses as its miss-index a weighted sum 
of the squares of the distances between the respective actual and 
desired B-plane asymptote pierce points. The weighting factors, 
which are supplied by the user, indicate the relative importance 
of securing nearby impacts at the respective target sites. NOMNAL 
computes its weighted least-squares solution by a hybrid pseudo- 
inverse and steepest-descent algorithm. The initial control it- 
erate is constructed by approximately targeting the first mini- 
probe to one of the target sites using a single Nev/ton-Raphson 
step. 

Finally the program integrates and records all segments of 
the nominal trajectory between guidance events from injection at 
the launch planet until the appropriate termination condition 
input by the user. For a conglomerate vehicle NOMNAL records 
the separate branches of the trajectory belonging to the main 
probe and miniprobes as well as to the bus. 



2.3 The Error Analysis /Generali zed Covariance Analysis Program 
ERRAN 

The error analysis/generalized covariance analysis program 
ERRAN is a preflight mission analysis tool that is used to deter- 
mine how selected error sources influence the orbit determination 
process for interplanetary or lunar missions. 

In the error analysis mode, ERRAN provides three primary quan- 
titative results: (1) knowledge covariance matrices, which pro- 

vide a measure of how well the actual trajectory is known, (2) 
control covariance matrices, which when propagated forward to the 
target provide a measure of how well the nominal target conditions 
will be satisfied by the actual trajectory, and (3) statistical 
midcourse AVs, which provide a measure of the amount of fuel re- 
quired for a successful mission. 

In the generalized covariance analysis mode, ERRAN provides 
all of the above information plus corresponding "actual" statis- 
tical information. The three results discussed in the previous 
paragraph are all computed on the basis of statistical distribu- 
tions assumed by the navigation filter to describe the significant 
error sources. In the generalized covariance analysis mode, "actual" 
knowledge covariances, control con variances , and statistical mid- 
course AVs are computed on the basis of statistical distributions 
that actually describe both error sources acknowledged by the 
navigation filter and the error sources ignored. The primary use 
of the generalized covariance analysis program is to study the 
sensitivity of filter performance to off-design conditions. 

ERRAN allows for employing gain generators for user-specified 
linear recursive navigation filters. Two gain generators are cur- 
rently available in ERRAN: (1) Kalman-Schmidt filter, and (2) 

equivalent recursive consider mode weigh ted- least-squares filter. 

State transition matrices are required to propagate covariance 
matrices over an arbitrary interval of time. Three methods are 
available for computing the 6x6 position/velocity state transition 
matrix. The first two methods, which are analytical methods, are 
analytical patched conic and analytical virtual mass. The third 
method uses numerical differencing to compute the state transition 
matrix. To increase the accuracy of the analytical techniques over 
long intervals, a state transition matrix cascading option is also 
available. Augmented parameter state transition matrices are al- 
ways computed using numerical differencing. 
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Up to 23 dynamic and measurement parameters may be solved-for 
or considered by the navigation filter. Parameters not acknowledged 
in design of the filter may be treated as ignore parameters when 
ERRAN is run in the generalized covariance analysis mode. The dy- 
namic parameters include biases in the gravitational constants of 
the sun and the target planet and biases in the six orbital ele- 
ments of the target planet. Measurement biases include biases in 
the locations of the three earth-based tracking stations, and biases 
in all measurements. Available measurement types are range, range- 
rate, star-planet angles, and apparent planet diameter measurements. 
Measurement noise for each measurement type is assumed to be con- 
stant. 

The computational procedure in ERRAN is divided into basic 
cycle computations and event computations. Basic cycle computa- 
tions are concerned with the propagation of covariances forward 
to a measurement time and processing the measurement. Events refer 
to a set of specialized computation, not directly concerned with 
measurement processing, that can be scheduled to occur at arbitrary 
times along the trajectory. 

The four events available in ERRAN are eigenvector, prediction, 
guidance, and probe release. At an eigenvector event the position 
and velocity partitions of the knowledge covariance matrix are di- 
agonalized to reveal geometric information about the size and ori- 
entation of the position and velocity navigation uncertainties. 
Associated hyperellipsoids are also computed. At a prediction 
event the most recent covariance matrix is propagated forward to 
some critical trajectory time to determine predicted navigation 
uncertainties in the absence of further measurements. 

The guidance event is the most complex event and yields much 
useful information for preflight mission analysis. Several types 
of guidance events are available in ERRAN. At a midcourse guidance 
event the user can choose from three midcourse guidance policies. 

The midcourse guidance event can also be constrainted to satisfy 
planetary quarantine requirements. At an orbital insertion guidance 
event the user can choose from two insertion policies. Options 
are also available for changing target conditions in midflight 
and retargeting the trajectory using nonlinear techniques, or for 
simply applying an externally supplied or precomputed AV at some 
arbitrary trajectory time. Two thrust models are available — impulse 
and impulse series. Execution error statistics are generated using 
an error model defined by a proportionality error, a resolution 
error, and two pointing angle errors. At a midcourse guidance 
event in ERRAN we also compute a statistical AV and the target 
condition covariance matrix both before and after the midcourse 
correction. 



Probe release events provide the capability to study missions 
employing multiprobe spacecraft. The multiprobe spacecraft is 
modeled as (1) a primary vehicle, or bus, with thrusting capabil- 
ity, (2) a main probe, with no thrusting capability, and (3) three 
miniprobes located symmetrically on booms attached to the bus, with 
no thrusting capability, and released simultaneously with AVs pro- 
vided by spinning the bus. Probe release events currently operate 
only in the error analysis mode of ERRAN. All measurement types 
and solve-for or consider parameters described previously are de- 
fined for all probes. Separate measurement schedules can be defined 
for the bus and the main probe. An additional measurement schedule 
can also be defined for all three miniprobes. Knowledge and con- 
trol covariances are propagated for each probe in sequential fash- 
ion. 


2.4 The Simulation Program SIMUL 

The simulation program SIMUL is the most complex program in 
the STEAP set of programs. In SIMUL the validity of the naviga- 
tion and guidance process is examined by simulating an actual mis- 
sion. Spacecraft state estimates are generated in SIMUL, as well 
as knowledge covariance matrices . The results given by the error 
analysis program ERRAN become meaningful only when SIMUL shows that 
the estimated spacecraft trajectory converges, within reasonable 
bounds specified by the covariance matrix, to the simulated actual 
trajectory. 

All state transition matrix, parameter augmentation, and meas- 
urement options described in section 3.3 are also available in 
SIMUL. As in ERRAN, the computational procedure in SIMUL is divided 
into basic cycle computations. The SIMUL basic cycle is concerned 
with the generation of state estimates and an actual trajectory, 
together with all quantities generated in the ERRAN error analysis 
basic cycle. Eigenvector and prediction events in SIMUL involve 
all computations performed in the corresponding ERRAN events. In 
addition, the SIMUL prediction event propagates state estimates 
forward to the time to which we are predicting. 


All options available in the ERRAN guidance event (see section 
3.3) are also available in the SIMUL guidance event. The treatment 
of the midcourse guidance event, however, is different in several 
respects. First, since an estimated spacecraft state is generated 
in SIMUL, an actual midcourse AV can be computed rather than a 
statistical AV as in ERRAN. Also, all linear midcourse AVs com- 
puted in SIMUL can be recomputed using nonlinear techniques. 
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Finally, since an actual trajectory is generated in SIMUL, actual 
target errors after the midcourse correction are also computed. 

Probe release events are also available in SIMUL. In addition 
to propagating knowledge and control covariance matrices for each 
probe, SIMUL also generates state estimates for each probe. 



3. 


INPUT DESCRIPTION 


3.1 NOMNAL Input Description 

The input for NOMNAL is transmitted via the namelist TAR1N 
and read in subroutine PRELIM. Each of the variables of TARIN will 
be described in full in this section. Many namelist variables will 
be specified by the program if they are not set by the user in the 
namelist input; these assumed values are the quantities enclosed 
in parentheses following the variable definition. 

Namelist TARIN: 


a) Nominal trajectory parameters 


ALNGTH 

TM 

ACKT 

NBOD 

NB(ll) 


Body 

Sun 

Mercury 

Venus 

Earth 

Mars 

Jupiter 

Saturn 

U ranus 

Neptune 

Pluto 

Moon 


- Length units per AU for output (ALNGTH 
= 149598500 : kilometers) 

- Time units per day for output (TM = 86400 
seconds) 

- Integration accuracy level used in nomi- 
nal trajectory propagation 

- Number of gravitation bodies included in 
integration 

- Array of codes (defined below) of gravi- 
tating bodies to be included in integra- 
tion 

Code 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 


NLP 


- Code of launch planet 


JO 



NTP - Code of target body 

IBARY - Flag determining inertial coordinate 

system (IBARY = 0) 

= 0 - Inertial system is heliocentric 
ecliptic 

= 1 - Inertial system is barycentric 
ecliptic 

NCPR - Number of integration increments between 

printouts (NCPR = 100) 

TMPR - Number of days between printouts 

(TMPR = 500.) 

SSS(3) - Array of direction cosines of spacecraft 

spin axis (SSS = (0,0,1)). 

b) Zero iterate parameters 

IZERO - Flag determining type of zero iterate 

generation (IZERO = 0) 

= 0 - Initial state read in through ZDAT 
( 1 - 6 ) 

= 1 - Launch planet at initial time to 
target body at final time 
= 2 - Launch planet at initial time to 

arbitrary point (specified by ZDAT 
(4-6)) at final time 

= 3 - Arbitrary initial point (specified 
by ZDAT (1-3) at initial time to 
target body at final time 
= 4 - Arbitrary initial point (specified 
by ZDAT (1-3) at initial time to 
arbitrary final point (specified 
by ZDAT (4-6)) at final time 
= 10 - Lunar targeting 

Calendar date of initial time (year, month, 
day, hour, minute, second) 

Specified to fractional second for 
IZERO = 0, 3, 4 

Specified only to day for IZERO = 1, 2 
Not required for IZERO = 10 (lunar tar- 
geting) 

Example: KALI = 1974, 5, 16, 0, 0, Si - 0.) 




KALT(5 ,10) , 
TS (10) 


ZDAT(6) 


RP 


- Calendar date of final time (year, month, 
day, hour, minute, second) in zero iterate 
computation (t^ for lunar zero iterate 
generation) 

- Vector defining zero iterate computation 
for interplanetary targeting 

ZDAT(l-3) = arbitrary initial position 
(inertial ecliptic) 

ZDAT(4-6) = arbitrary final position 
(inertial ecliptic) 

for lunar targeting 

ZDAT(l-3) = desired semimajor axis, 

radius of closest approach, 
and equatorial inclination 
of lunar conic 

- Parking orbit radius in launch model 
(RP = 6560 km) 


SIGMAL - Nominal value of launch aximuth which 

may be varied if necessary (SIGMAL = 90°) 

KOAST - Flag determining launch plane (KOAST = 1) 

=1 - long coast time orbit 

= -1 - short coast time orbit 


The following parameters are used to define the launch profile for 
interplanetary trajectories: 

FI - True anomaly at injection (FI = 3.7°) 

PSI1 - Angle of first burn (PSI1 = 17.°) 

PS I 2 - Angle of second burn (PSI2 = 8.°) 

TIM1 - Duration of first bum (TIM1 = 500.) 

TIM2 - Duration of second burn (TIM2 = 100.) 

THELS - Longitude of launch site (THELS = 279.457°) 

PHILS - Latitude of launch site (PHILS = 28.317°) 

THEDOT - Rotation rate of launch planet 

(THEDOT *= 15.041) 

RPRAT - Inverse parking orbit rate (RPRAT = 14.689) 
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Up to 10 guidance events are permitted in any run of NOMNAL. The 
guidance event parameters are defined as arrays; parameters cor- 
responding to the same event occur in the same component of all 
arrays. The indices of the events are not required to be consecu- 
tive nor must they be in chronological order. 


c) General guidance event parameters 

KTYP(10) - The types of guidance events 


= -1 


= 1 


= 2 


= 3 


= 4 


= 5 


- Termination event (each run must 
have such an event to terminate 
the integration of the nominal 
trajectory) 

- Targeting event (the current ve- 
locity will be directly refined 
to yield a trajectory satisfying 
input target conditions) 

- Retargeting event [the current 
position and target state will be 
used to compute a zero iterate 
velocity (see IZER below) that 
will then be refined to meet pre- 
scribed target conditions] 

- Orbit insertion event (using the 
current state, the velocity cor- 
rection needed at a later time to 
insert into a desired orbit is 
computed and executed at that la- 
ter time) 

- Main probe propagation event (the 
current state is stored so it can 
be returned to after the event. 

Then it is propagated to a stop- 
ping condition generating a printed 
time history. Finally the original 
current state is restored in prep- 
aration for the next event. This 
event type is used in treating 
branched trajectories) 

- Miniprobe release event (the cur- 
rent state is stored. Then the 
release controls are calculated to 
apply at the current time to tar- 
get three miniprobes to three re- 
spective target sites characterized 
by input values of declination and 



KTIM(IO) 


TIMG(IO) 


KALG(5 ,10) 
GS(10) 


IZER(IO) 


right ascension. Next each of the 
miniprobes is propagated from re- 
lease to impact using the minimum- 
miss controls. Time histories of 
these trajectories are printed dur- 
ing their generation. Finally the 
current state is restored for the 
next event) 

(KTYP = 1,1,1,...) 

- Codes defining epochs to which times of 
guidance events are referenced 

= 0 - Event is not to be processed 
= 1 - Epoch is initial time 
= 2 - Epoch is time of intersection of 

sphere of influence of target body 
= 3 - Epoch is time of closest approach 
to target body 

= 4 - Time is read in as calendar date 
in KALG, SG 
(KTIM = 0,0,0, . . .) 

- Time intervals (days) of guidance events 
past reference epochs (KTIM) 

(TIMG = 0. ,0. ,0. , . . .) 

- Dates of guidance events (year, month, 
day, hour, minute, seconds) used only 
for KTIM(i) = 4. 

Example: KALG(l,i) = 1975,10,12,2,14, 

GS(i) = 10.234 

- Flag determining method of computing 
zero iterate velocity for retargeting 
event 

= 3 - Current position to target body 
at target time 

= 4 - Current position to target posi- 
tion (KTAR = 10,11,12) at target 
time 

(IZER = 0,0,0, . . .) 



KMXQ(IO) - Compute/execute modes of events 
(KMXQ = 3,3,3,...) 

= 1 - Compute velocity correction only 
- 2 - Execute velocity correction only 
= 3 - Compute and execute velocity cor- 
rection at same time 
= 4 - Compute velocity correction but 
execute at a later time (orbit 
insertion) 

DELV(3,10) - The impulsive velocity increment in 

ecliptic coordinates for each event i 
such that KMXQ(i) = 2 (km/s) . 

Example: DELV(l,i) = .050, + .115, - .007 

MDL(10) - Execution models for velocity cor- 

rections (MDL = 1,1,1,...) 

= 1 - Impulsive model 
= 2 - Multiple pulsing arc 

The following parameters need be set only if the multiple pulsing 
arc execution model is used in the run. 


PULMAG - Thrust magnitude (T) of pulsing engine 

PULMAS - Nominal mass in of spacecraft during 

pulsing arc 

DUR - Time duration (At) of single pulse 

DTI - Time interval (At^) between pulses (days) 

The units of the first three variables must be such that the ve- 
locity increment imparted on a single pulse Av = T/m At is in km/s. 


KLP(10) 

KTP(10) 


- Flags used to change the launch and tar- 
get bodies at a guidance event 

= 0 - Use previous codes NLP and NTP 
= K, - Replace previous code by K. 


Example: KLP = 0,0,5, KTP = 0,0,6, deter- 

mines that at the third event Mars be- 
comes the launch planet and Jupiter the 
target planet. 

(KLP = 0,0,..., KTP = 0,0,...) 


15 



KTAR(6 , 10) 


- Targeting event: Target parameter codes 

for each targeting event using the defi- 
nitions listed below. Example: KTAR(l,i) 

= 7,8,3, triggers target parameters r , 

it LA 

CA* CS* 


Code 

Symbol 

Def ini tion 

1 

TPS 

Time at impact in Julian days 
epoch 1900 extrapolated conically 
from SOI conditions.* 

2 

TSI 

Time at SOI 

3 

TCS 

Time at CA extrapolated from SOI 
conditions 

4 

TCA 

Time at CA 

5 

BDT 

B*T at target time 

6 

BDR 

B*R at target time 

7 

RCA 

Radius of closest approach com- 
puted at target time 

8 

INC ' 

Inclination (planet equatorial) 
at target time 

9 

ASI 

Semimajor axis (computed from SOI 
conditions) 

10 

XRF 

X-component of final position 
(inertial ecliptic) 


If 1 is among the target codes, 13 and 14 must be the remaining 
two. Then all the target values, like the 1 value, are extrapolated 
from SOI conditions. 

•j* 

The inclination must be specified as follows. Let a be the 
magnitude of the inclination with 0° < a i 90°. Then (see IMPACT 
for further details) 

6=0 posigrade motion 6 = +1 northern approach 

= 1 retrograde motion = -1 southern approach 

i “ 6 180° 4- 6 a 

M A 
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Code 

Symbol 

Definition 

11 

YRF 

Y-component of final position 
(inertial ecliptic) 

12 

ZRF 

Z-component of final position 
(inertial ecliptic) 

13 

DCP 

Declination in degree of probe at 
impact relative to the probe co- 
ordinate system specified by 1PCS 

14 

RAP 

Right ascension in degree of probe 
at impact relative to the probe 
coordinate system specified by 
IPCS 

15 

TPR 

Time at impact in Julian days 


epoch 1900 obtained completely by 
integration 


Orbit insertion event : The flag deter- 

mining whether the coplanar (=* 1) or 
nonplanar (= 2) option is to be used 
on an orbit insertion event. Example: 
KTAR(l,i) = 1 triggers the coplanar 
option 

TAR(6,10) - Targeting event : The desired target 

parameter values for each targeting 
event in the same order as the target 
parameter codes. Example: TAR(l,i) 

= 7000. , 50. , sets the desired values 
r = 7000 km, i = 50° if KTAR(l,i) = 7,8 

vA 

Coplanar orbit insertion : (KTYP(i) = 3, 

KTAR(l,i) = 1). The desired values of 
semimajor axis (km) , eccentricity, and 
periapsis shift for each coplanar in- 
sertion. Example: TAR(l,i) = 20000., 

.75, 5 

Nonplanar orbit insertion : (KTYP(i) = 3, 

KTAR(l,i) = 2). The desired values of 
semimajor axis (km) , eccentricity, ar- 
gument of periapsis (equatorial) , in- 
clination (equatorial) , and longitude 
of ascending node (equatorial) for each 
nonplanar insertion. Example: 

TAR(l,i) = 20000., .75, 30., 15., 210 
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KTAR(l.lO) 


KTAR(2 , 10) 


TAR(j.lO) 


TAR( j+3» 10) 


TOL( j , 10) 


Miniprobe targeting ; The contents of 
the arrays KTAR, TAR, TOL, and WGHTJi 
are redefined as follows for this event 
type. 

- Flag specifying the spin axis orientation 
mode 

= 11 - Ecliptic declination and right 
axcension of spin axis are free 
controls 

* 12 - Spin axis is aligned with space- 
craft velocity vector at release 
= 13 - Spin axis is perpendicular to 

spacecraft sun line, parallel to 
ecliptic plane, and within 90° of 
spacecraft velocity vector at re- 
lease 

= 14 - Ecliptic declination and right 
ascension of spin axis are each 
fixed at values input by the user 

- Flag specifying the mode of miniprobe 
propagation 

= 11 - Conic model alone 
= 12 - Conic initial iteration followed 
by virtual-mass refinement 

- Desired declination in degree of impact 
point of jth miniprobe in coordinate 
system opted by IPCS(i) and for planet 
radius specified by RPS(i) 

- Desired right ascension in degree of 
impact point of jth miniprobe in the co- 
ordinate system and at the planet radius 
described above for TAR(j,10) 

- Weighting factor applied to the B*T and 
B*R errors of the jth target site in the 
miss-minimizing algorithm. The least 
important site should be assigned a unity 
weight, while the more important ones 
should be given progressively larger 
values 



T0L(4, 10) 


TOL(5 , 10) 


TOL(6 , 10) 


WGHTM( 10) 


RPS(10) 


IPCS(IO) 


KALT(5 , 10) , 
TS(10) 


- Tolerance, e, on release control con- 
vergence. For the actual convergence 
criterion, consult the STEAP Analytical 
Manual under miniprobe targeting. A 
suggested value for e is 1 

- Declination in degrees of the spin axis 
for the fixed inertial-orientation spin 
axis mode 

- Right ascension in degrees of the spin axis 
axis for the fixed inertial-orientation 
spin axis mode 

- Length of maximum pseudoinverse step in 
the control space. Longer steps are de- 
ferred in favor of a best-step steepest- 
descent correction. A suggested value 
is 0.5 radians or decameters 

- The radius of the probe impact sphere in 
km. It must be loaded for each targeting 
event to planet impact parameters, each 
main probe propagating event, and each 
miniprobe targeting event 

- The flag specifying the planetocentric 
coordinate system for the various types 
or probe events. It must be loaded for 
the same cases as is RPS ( XP£> — O, 0j Oj ' ’ “ 

0 = equatorial 

1 = subsolar orbit plane 

- Desired value of target time in calendar 
dates for targeting event regardless of 
whether the target time is t t , or 

O J- vb 

t . Thus, if KALT(l,i) = 1975,9,18,0,0, 

TS(i) = 0., and KTAR(3,i) =2, the de- 
sired time at the target planet SOI is 
September' 18, 1975. If KTAR(3,i) «= 4, 
that date is to be the desired time of 
closest approach. If IZERO =1, 2, 
KALT(5,0) and TS(1) must correspond to 
the final time of the zero iterate 
computation 



T0L(6 , 10) 


NPAR(10) 


Targeting Scheme 
METH(IO) 


WGHTM(IO) 

PERV(IO) 


■ Tolerances on desired values for tar- 
geting events in same order as the tar- 
get parameter codes. Thus TOL(l,i) = 100., 
1. , .005 sets the tolerances of 100 km 

in r , 1° in i, and .005 days in t 

V-»A L*o 

if KTAR(l,i) = 7,8,3 

■ The number of targeting parameters to 
be targeted for each targeting event. 

If NPAR(i) = 3, all three velocity com- 
ponents will be refined to meet the 
three velocity components will be re- 
fined to meet the three target parameters. 
If NPAR(i) = 2, the X- and Y-components 
will be refined to meet the first two 
target parameters (NPAR = 3,3,3,...) 

Parameters 

- The method to be used in the targeting 

= 0 - Use Newton-Raphson targeting matrix 
method 

= n - Use cycle of n-1 conjugate gradient 
steps and one steepest descent step 
repeatedly . 

Example: METH = 0,1,10 specifies Newton- 

Raphson technique for first targeting 
event, steepest descent for second event, 
and a conjugate gradient technique recti- 
fied by steepest descent every 10 steps 
for third event (METH = 0,0,0,...) 

- The weighting factor used for the time 
variable in assigning a scalar loss func- 
tion for the auxiliary parameters for both 
the bad-step logic and gradient compu- 
tations (WGHTM = 10 5 , 10 5 , 10 5 ,...) 

- The velocity perturbations (km/s) used 
for each targeting event to compute 
either the targeting matrix or the gra- 
dient (PERV = .00005, .00005,...) 



DVMAX(IO) 

SPHFAC(IO) 

NOIT(IO) 

MAXB(IO) 

IB ADS (10) 

1ST ART 


PHI(3,3) 


- The maximum step allowed (km/s) on any 
iterate (DVMAX = 0. 1,0. 1,0. 1, . . .) 

- The factors by which the target planet 
sphere of influence is to be reduced at 
each targeting event. Thus if SPHFAC(i) 
= .5, the SOI will be reduced by half 

to do the targeting at the ith event 
(SPHFAC = l.,l.,l ) 

- The number of total iterations allowed 
at the first and last level of the tar- 
geting events (NOIT = 8,8,8,...) 

- The number of bad steps allowed during 
any targeting event (MAXB = 4,4,4,...) 

- The bad-step flags for each targeting 
event 

= 1 - Never use bad-step check 
= 2 - Use bad-step check at final level 
only 

= 3 - Use bad-step checks throughout 
(IBADS =3,3,3,...) 

- Stage of first targeting event 

= 0 - Compute targeting matrix on first 
iteration 

= 1 - The first phase of the targeting 
has been started and a valid tar- 
geting matrix for the first phase 
will be read in as PHI 
= 2 - The second phase of the targeting 
has been started and a valid tar- 
geting matrix for the second phase 
will be read in as PHI 
(ISTART = 0) 

- The targeting matrix to be used re- 
peatedly, defined by the value of ISTART 



MAT (10) 


- Targeting matrix computation code for 
each targeting event 

= 1 - Compute targeting matrix only at 
first level 

= 2 - Compute targeting matrix at each 
step 

(MAT = 1,1,1,...) 

AC(5,10) - The accuracy levels used for each event. 

The final accuracy level at each guidance 
event should be identical to the tra- 
jectory accuracy level ACKT. Thus if 
AC ( 1 , i ) - l.E-4, 2 .5E-5 , 5.E-6, the ith 
guidance event will be targeted at those 
progressive levels 

- The number of accuracy levels used for 
each targeting event (LVLS = 3,3,3,...) 

- The flag used in designating the controls 
to be used in calculating the sensitivity 
matrix in a targeting event 

= 1 0 - Intertial x, y, and z spacecraft 
velocity components 

= 2 ®- Magnitude of the spacecraft veloc- 
ity relative to the launch planet, 
and in-plane and out-of-plane ro- 
tation angles from the current 
relative velocity* 

(CON = 2,^2,,...) 


See the analysis section of the subroutine KTROL for a de- 
tailed description of these controls. 


LVLS (10) 

-X 

-€0N(1O) 
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3.2 ERRAN Input Description 


The input of the error analys is /generalized covariance analysis 
program consists of: 

a) A card containing the variable IRUNX (110 field) 
that indicates how many different runs are to be 
made and is read only once; 

b) A card containing the problem identification variable 
(110 field) that precedes each set of input data; 

c) An error analysis namelist section entitled ERRAN; 

d) Three successive measurement schedules for the pri- 
mary vehicle, main probe, and miniprobes in that order 
(see namelist variables NENT, NENT1, and NENT2) ; 

e) A generalized covariance analysis namelist section 
entitled GENRAL , that must appear only if a generalized 
covariance analysis is to be performed. 

Most namelist variables are preset by the program; these preset 
values are the quantities enclosed in parentheses in the namelist 
definitions. Unless otherwise indicated, input units correspond 
to the internal units defined by the variables ALNGTH and TM. 
Unspecified angular units are assumed to be radians. 

3.2.1 Namelist ERRAN 

1. Nominal trajectory variables 

XI (6) - Initial position/velocity state of 

spacecraft; not specified if 
IC00R =3 (XI = 6 * 0.) 

IC00R - Code that specifies coordinate system 

of initial spacecraft state. 

(IC00R = 2) 

= 0, heliocentric ecliptic 

= 1, geocentric equatorial 

= 2, geocentric ecliptic 

= 3, JPL conditions: RDS , PHIT, 

THETA, VEL, CAMMA, SIGMA 
=" 4, planetocent ric ecJIptic (target 

planet) 



= 5, planetocentric equatorial 

= 6, planetocentric orbital elements: 

semimajor axis, eccentricity, 
inclination, longitude of the 
ascending node, argument of 
periapsis, and true anomaly 


The following six variables define the JPL conditions: 


RDS 

PHIT 

THETA 

VEL 

GAMMA 

SIGMA 

LM0 

LDAY 

LHR 

SECL 

LYR 

IM0 

I DAY 

IHR 

IMIN 

SECI 

IYR 

ALNGTH 


- Earth - centered injection radius 

- Injection declination (degrees) 

- Injection right ascension (degrees) 

- Injection velocity relative to the 
Earth 

- Injection flightpath angle (degrees) 

- Injection azimuth (degrees) 

- Launch month (integer) 

- Launch day (integer) 

- Launch minute (integer) 

- Launch second (floating) 

- Launch year (integer) 

- Month of final computation (integer) 

- Day of final computation (integer) 

- Hour of final computation (integer) 

- Minute of final computation (integer) 

- Second of final computation (floating) 

- Year of final computation (integer) 

- Length units per AU (ALNGTH = 149598500. 
Kilometers) 
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- Time units per day. (TM = 86400. 
seconds) 

- Initial trajectory time. (TRTMI =0.) 

- Injection trajectory time. (TINJ ■ 0.) 

- Nominal trajectory module code. 

(NTMC = 2) 

= 1 - patched conic (not supplied 

with this program) 

= 2 - virtual mass 

- Number of celestial bodies considered 
in the generation of the nominal tra- 
jectory (NB0D = 3) 

- Array of codes of celestial bodies 
considered in the generation of the 
nominal trajectory 

= 1 - Sun 

= 2 - Mercury 

= 3 - Venus 

= 4 - Earth 

= 5 - Mars 

= 6 Jupiter 

= 7 - Saturn 

= 8 - Uranus 

= 9 - Neptune 

- 10 - Pluto 
= 11 - Earth's Moon 


NLP 

- Launch planet code 

NTP 

- Target planet code 

IEPHEM 

- Ephemeris code (IEPHEM = 1) 


= 0 - Place each planet in an ellipse. 

The date at which this ellipse 
is calculated is determined by 
reading in a variable entitled 
as the first six letters of 
the name of the planet consid- 
ered. This variable should 
contain six integers specifying 
the month, day, hour, minute, 
second, and year. Example: 

EARTH - 7, 24, 6 15, 38, 1973. 

= 1 - Calculate orbital elements for 

each planet at each time interval 


TM 

TRTMI 

TINJ 

NTMC 


NB0D 

NB(ll) 



SSS(3) 


ACC 

ISP2 

IB ARY 

State transition 
ISTMC 

DTMAX 

NDACC 


- Array of direction cosines of space- 
craft spin axis (SSS = 0.0.1) 

- Nominal trajectory accuracy figure 
(ACC = 1.0 x 1(T 6 ) 

- Code of virtual mass trajectory 
(ISP2 =0) 

= 0 - Continue integrating to final 
time 

- 1 - Stop integrating when target 

planet sphere of influence is 
encountered 

- Reference coordinate system option 
(IBARY = 0) 

= 0 - Reference coordinate system is 
heliocentric ecliptic 
= 1 - Reference coordinate system is 
barycentric ecliptic 

matrix variables 

- State transition matrix code 
(ISTMC = 1) 

= 1 - Analytical patched conic 

- 2 - Analytical virtual mass 

= 3 - Numerical differencing using 
virtual mass 

- Maximum time interval for which analy 
tical computation of the state transi 
tion matrix is considered valid 
(DTMAX = 8. days) 

- Accuracy code for numerical differ- 
encing (NDACC =0) 

= 0 - Same accuracy as is employed 

in the computation of the nomi- 
nal trajectory 

= 1 - Accuracy = ACCND S described 
next 



ACCND 


- Accuracy to be used in the calculation 
of the state transition matrix by nu- 
merical differencing (ACCND = 2.5 x 10 -5 ) 

ISTM1 - Cascaded state transition matrix code 

(ISTM1 = 1) 

(ISTM1 should be set to 3 for lunar 
missions since cascading option has 
not been defined for such missions) 

= 1 - Patched conic Danby method 
= 2 - Virtual-mass Danby method 
= 3 - Numerical differencing (CASCAD 
not called) 

DTSUN - Integration interval when sun is cen- 

tral body and ISTM1 = 1. 

(DTSUN =2.0 days) 

DTPLAN - Integration interval when target planet 

is central body and ISTM1 = 1. 

(DTPLAN =0.25 days) 

FACP - Position factor for numerical differ- 

encing (FACP = 1.) 

FACV - Velocity factor for numerical differ- 

encing (FACV = 1. x 10 -1+ ) 

The following eight variables are used to compute the augmented 
state transition matrix by numerical differencing: 

DELMUS - Sun gravitational constant factor; 

need be specified only if IAUGIN(IO) 
is nonzero (DELMUS = 1. x 10 7 ) 

DELMUP - Target planet gravitational constant 

factor; need be specified only if 
I AUGIN ( 11) is nonzero (DELMUP =0.1) 

DELAXS - Target planet semimajor axis factor; 

need be specified only if IAUGIN(12) 
is nonzero (DELAXS = 100.) 
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DELECC - Target planet eccentricity factor; 

need be specified only if IAUGIN(13) 
is nonzero (DELECC = 1. x 10 -5 ) 

DELICL - Target planet inclination factor; 

need be specified only if IAUGIN(14) 
is nonzero (DELICL = .0000484814 
radians) 

DELN0D - Target planet longitude of the as- 

cending node factor; need be specified 
only if IAUGIN(15) is nonzero 
(DELN0D » .0000484814 radians) 

DELW - Target planet argument of periapsis 

factor; need be specified only if 
IAUGIN(16) is nonzero (DELW = .0000484814 
radians) 

DELMA - Target planet mean anomaly factor; 

need be specified only if IAUGIN(17) 
is nonzero (DELMA = .0000484814 radians) 


Parameter augmentation variables 


IAUGIN(24) - Array of augmented parameter codes; 

unspecified elements are assumed to 
be zeros. Up to 12 solve-for parameters 
may be augmented; up to 8 dynamic- 
consider parameters; up to 15 measurement- 
consider parameters; and up to 12 
ignore parameters 


IAUGIN(I) = 0 - neglected parameter 

= 1 - consider parameter 
= 2 - solve-for parameter 
= 3 - ignore parameter (generalized 
covariance only) 


I = 


1 

2 

3 

4 

5 

6 

7 

8 
9 


Radius error of station 1 
Latitude error of station 1 
Longitude error of station 1 
Radius error of station 2 
Latitude error of station 2 
Longitude error of station 2 
Radius error of station 3 
Latitude error of station 3 
Longitude error of station 3 


measurement 

parameters 



I = 10 
11 


12 

13 

14 

15 

16 
17 


Sun gravitational constant bias 
Target planet gravitational 
constant bias 

Target planet semimajor axis 
bias 

Target planet eccentricity bias 
Target planet inclination bias 
Target planet longitude of 
ascending node bias 
Target planet argument of 
periapsis bias 

Target planet mean anomaly bias 


dynamic 

parameters 


18 

19 

20 
21 
22 

23 

24 


Range bias of station 1 
Range-rate bias of station 1 
Star-planet angle 1 bias 
Star-planet angle 2 bias 
Star-planet angle 3 bias 
Apparent planet diameter bias 
Undefined 


I 


measurement 

parameters 


4. Measurement variables 

NENT - Number of entries (cards in primary 

vehicle measurement schedule (NENT = 0) 


NENT1 - Number of entries (cards) in the main 

probe measurement schedule (NENT1 = 0) 

NENT2 - Number of entries (cards) in the sin- 

gle measurement schedule for all three 
miniprobes (NENT2 = 0) 


NST - Number of tracking stations (at most 3) 

on the rotating earth (NST = 3) . If » 
no tracking station information is read 
in, the following three stations will 
be assumed: 




Altitude 

Latitude 

Longitude 

1. 

Golds tone 

1.031 km 

35.384 N 

116.833 W 

2. 

Madrid 

.050 km 

40.417 N 

3.667 W 

3. 

Canberra 

.050 km 

35.311 S 

149.136 E 


2 7-2 



If different tracking stations are desired, their locations must 
be specified by the following three arrays. 

SAL(3) - Array of altitudes of each tracking 

station 

SLAT(3) - Array of latitudes of each tracking 

station in degrees north 

SL0N(3) - Array of longitudes of each tracking 

station in degrees east 

Direction cosine arrays of three 
reference stars. If not specified, 
the three stars and their direction 
cosines are as follows: 

Canopus Betelgeuse Rigel 

UST -.061351 .028986 .201963 

VST .237886 .960388 .831343 

WST -.969355 -.277141 -.517784 

5. Eigenvector and prediction event variables 

NEV1 . - Number of eigenvector events (NEV1 = 0) 

NEV2 - Number of prediction events (NEV2 = 0) 

Tl(20) - Array of times at which eigenvector 

events occur; specified only if NEV1 
is nonzero. Chronological order 
required 

T2(20) - Array of times at which prediction 

events occur; specified only if NEV2 
is nonzero. Chronological order 
required 

TPT2(20) - Array of times to which one wishes 

to predict. The elements of the TPT2 
array must correspond to the elements 
of the T2 array and must be specified 
only if the T2 array has been specified 


UST(3) 
VST( 3) 
WST( 3) 
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IEIG 


- Eigenvector code (IEIG = 1) 


= 0 - Only position eigenvectors will 
be calculated 

= 1 - Both position and velocity eigen- 
vectors will be calculated 

IHYP1 - Hyperellipsoid sigma level code 

(IHYP1 = 2) 

= 1 - Sigma level equals one 
= 2 - Sigma level equals three 
= 3 - Sigma levels of both one and 
three 

F0P - A value to be used as an off-diagonal 

annihilation element in subroutine 
JAC0BI for position eigenvalues and 
eigenvectors (F0P = 1. x 10~ 15 ) 

F0V - A value to be used as an off-diagonal 

annihilation element in subroutine 
JAC0BI for velocity eigenvalues and 
eigenvectors (F0V = 1. x 10 -25 ) 

6. Covariance variables (filter, or assumed, covariances) 

P(6,6) - Initial P (position and velocity) co- 

variance matrix. Referenced to in- 
ertial frame (diag P = l.,l.,l., 

1. x 10 -4 , 1. x 1CT 4 , 1. x 10~ 4 ) 

The structure of the following eight parameter covariance matrix 
partitions must correspond to the structure of the solve-for, 
dynamic-consider, and measureraent-consider parameter vectors. 

PS(12,12) - Initial P (solve-for parameter) co- 

b 

variance matrix (PS = identity matrix) 

U0(8,8) - Initial U q (dynamic-consider param- 

eter) covariance matrix (UO = identity 
matrix) 

VO(15,15) - Initial V q (measurement-consider 

parameter) covariance matrix 
(VO = identity matrix) 
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CXXS (6 , 12) 


- Initial C 


covariance matrix 


(CXXS = 0) 

CXU(6,8) - Initial C covariance matrix 

xu 

(CXU = 0) 

CXV(6,15) - Initial C covariance matrix 

xv 

(CXV = 0) 

CXSU(12,8) - Initial C covariance matrix 

x u 
s 

(CXSU = 0) 

CXSV(12,15) - Initial C covariance matrix 

x v 
s 

(CXSV = 0) 

IGAIN - Filter gain generator code (IGAIN = 1) 

= 1 - Kalman-Schmidt filter 
= 2 - Equivalent recursive weighted 
least-squares consider filter 
(available only if subroutine 
GAIN2 has been loaded) 

IDNF - Dynamic noise flag (IDNF = 0) 

= 0 - Dynamic noise is zero 
= 1 - Dynamic noise is not zero 

DNCN(3) - Array of constants used to calculate 

dynamic noise covariance matrix; must 
be specified if IDNF equals 1 

IMNF - Measurement noise flag (IMNF = 0) 

= 0 - Measurement noise is constant 
= 1 - Measurement noise is not con- 
stant (option is not available 
with this program) 
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MNCN( 12) 


SIGRES 

SIGPR0 

SIGALP 

SIGBET 

PSIGS 

PSIGK 

PS IGA 


- Array of variances for each type of 
measurement. If not specified, the 
following values are assumed: 


MNCN(l) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


1. x 10 -8 Range (idealized 
station) 

1. x 10~ 12 Range rate 
(idealized station) 

1. x 10“ 8 Range (station 1) 

1. x 10~ 12 Range rate 
(station 1) 

1. x 10 -6 Range (station 2) 

1. x 10” 12 Range rate 
(station 2) 

1. x 10“ 6 Range (station 3) 

1 / x 10" 12 Range rate 
(station 3) 

2.5 x 10 -9 Star-planet angle 1 
2.5 x 10“ 9 Star-planet angle 2 
2.5 x 10 -9 Star-planet angle 3 
2.5 x 10 -9 Apparent planet 
diameter 


- Variance of resolution execution error 
(SIGRES = 4. x 10- 8 ) 

- Variance of proportionality execution 
error (SIGPR0 = .0001) 

- Variance of pointing angle alpha exe- 
cution error (SIGALP = .0043625 radians 2 ) 

- Variance of pointing angle beta execu- 
tion error (SIGBET = .0043625 radians 2 ) 

- Variance of resolution execution error 
for pulsing engine (PSIGS = 4. x 10~ 8 ) 

- Variance of proportionality execution 
error for pulsing engine (PSIGK = .0001) 

- Variance of pointing angle alpha exe- 
cution error for pulsing engine 
(PSIGA = .0043625 radians 2 ) 
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PSIGB 


I GEN 


Print codes 
DELTP 

INPR 

IPRINT 


KPRINT 


IPRT 


- Variance of pointing angle beta exe- 
cution error for pulsing engine 
(PSIGB = .0043625 radians 2 ) 

- Code that indicates if a generalized 
covariance analysis is to be performed 
(IGEN =0) 

= 0 - No generalized covariance analysis 
= 1 - Generalized covariance analysis 


- Trajectory print interval in days 
(DELTP = 1. x 10 50 ) 

- Trajectory print interval in incre- 
ments (INPR = 100000) 

- Measurement print interval; measure- 
ment information printed every IPRINT 
measurements (IPRINT = 1) 

- Correlation matrix print code 
(KPRINT = 0) 

= 0 - Print out P and P correlation 

u 

matrices and standard deviations 
at a measurement 

= 1 - Print out all correlation ma- 
trices and standard deviations 
at a measurement 

- Array of print codes (IPRT = 1, 1,1,1) 

IPRT(I) = 1 - Print out information 

= 0 - Do not print out infor- 
mation 

1=1 Ephemeris data 

2 Spacecraft trajectory relative 
to planets 

3 Virtual-mass data 

4 Navigation parameters (TRAPAR 
called) 



8. Guidance event variables 


NEV3 
T3( 10) 


ICDQ3(10) 


- Number of guidance events (NEV3 = 0) 

- Array of times at which guidance 
events occur; specified only if NEV3 
is nonzero. Chronological order 
required 

- Array of codes for guidance events 
to determine how the execution error 
covariance matrix is to be calculated. 
These codes must correspond to the 
elements of the T3 array and need be 
specified only if the T3 array has 
been specified (ICDQ3 = array of l's) 

= 0 - Calculated directly from ve- 
locity correction covariance 
matrix 

= 1 - Calculated from the eigenvector 
corresponding to the maximum 
eigenvalue of the velocity cor- 
rection covariance matrix 


IGUID(5, 10) 


IGUID(ljI) 


IGUID(2,I) 


- Array of guidance event codes. This 
array is defined more explicitly be- 
low for the I-th guidance event oc- 
curring at time T3(I) 


= 1 
2 

3 

4 

5 

6 
7 


Fixed-time-of-arrival midcourse 
guidance event 

Two-variable B-plane midcourse 
guidance event 

Three-variable B-plane midcourse 
guidance event 

Planar orbital insertion deci- 
sion event 

Nonplanar orbital insertion 
decision event 

Externally supplied velocity 
change event 
Retargeting event 


= 0 - Linear guidance 
1 - Nonlinear guidance 
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IGUID( 3, I) 


IGUID(4,I) 


IGUID(5,I) 


= 0 - Planetary quarantine constraints 
not in effect 

1 - Planetary quarantine constraints 

in effect; use linear guidance 
to achieve biased aimpoint 

2 - Planetary quarantine constraints 

in effect; use nonlinear guidance 
to achieve biased aimpoint 

= 1 - Single impulse thrust model 

2 - Pulsing thrust model 

3 - Finite burn (not available) 

= 1 - Execute event only 

2 - Compute event only 

3 - Compute and execute event 

4 - Compute, but execute event later 


The IGUID(J,I) array cannot be chosen arbitrarily, but must con- 
form to the following table of permissible combinations. 


IGUID(1,I) 

IGUID(2,I) 

IGUID( 3,1) 

IGUID(4,I) 

IGUID(5 , I) 

1, 2, 3 

0 

0, 1, 2 

1, 2 

2, 3 

4, 5 

— 

0 

1, 2 

2, 4 

6 

— 

0 

1, 2 

1 

7 

— 

o 

1, 2 

2, 3 


Additional restrictions on the IGUID(J,I) array are listed below: 

a. ICDQ3(I) must be set to 1 if IGUID(3,I) + 0. 

b. Only one orbital insertion event may occur, and it 
must be the final guidance event. 

c. An externally supplied velocity change event may not 
be preceded by a midcourse guidance event if the ori- 
ginal nominal does not pierce the target planet sphere 
of influence. Note also that nominal target conditions 
TN0MB and TN0MC are not altered when an externally 
supplied velocity change event occurs. 
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No additional guidance variables are required if the guidance 
event is a linear impulsive midcourse guidance event subject to 
no planetary quarantine constraints. Other types of guidance 
events require that some of the following guidance variables be 
specified: 


LKTAR(6 ,10) 
XTAR(6 , 10) 
XPERV(IO) 

XDVMAX(IO) 
LNPAR(IO) 
LLVLS(IO) 
TGT3(10) 
DELV(3, 10) 
TN0MB ( 3) 
TN0MC( 7) 

PR0BI 

IDENS 

PULMAG 


- Array defining target parameters, 
identical to KTAR in NOMNAL 

Desired target value; identical to 
TAR in N0MNAL 

- Velocity perturbation used to compute 
targeting matrix; identical to PERV 
in N0MNAL 

- Maximum allowable velocity correction; 
identical to DVMAX in N0MNAL 

- Number of target parameters desired; 
identical to NPAR in N0MNAL 

- Number of integration accuracy levels 
used; identical to LVLS in N0MNAL 

- Desired target times; referenced to 
initial trajectory time 

- Array of externally supplied velocity 
changes (DELV = array of zeros) 

- Nominal B-plane target conditions: 

B*T, B*R, and tg^. 

- Nominal closest approach target con- 

ditions: R^, V^, and t^ 

- Allowable probability of impact 
(PR0BI = 1.) 

- Code defining method of treating 
probability density function in sub- 
routine BIAIM. Option not available . 
at present (IDENS = 1) 

- Magnitude of pulsing engine thrust 
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PULMAS 


- Nominal mass of spacecraft 


DUR - Duration of single pulse 

DTI - Time interval between pulses 

List of required variables for guidance events other than linear 
impulsive midcourse guidance events subject to no planetary quar- 
antine constraints: 

a. Retargeting: LKTAR, XTAR, XT0L, XAC, LNPAR, XPERV, 

XDVMAX, LLVLS , TGT3. 

b. Nonlinear guidance: XT0L, XAC,, XPERV, XDVMAX, LLVLS, 

TN0MB, TN0MC. 

c. Orbital insertion: XTAR. 

d. Biased aimpoint guidance (planetary quarantine) : 

PR0BI , TN0MB, TN0MC, TINJ , I DENS. 

. e. Pulsing thrust model: PULMAG, PULMAS, DUR, DTI, PSIGS, 

PSIGK, PSIGA, PSIGB. 

f. Externally supplied velocity change: DELV. 

9. Probe release events 


T6 

T7 

PMN(12) 


IPMN 


- Time of main probe release event 

- Time of miniprobe release event 

- Array of main probe measurement noise 
variances for each type of measurement. 
PMN(I) refers to same measurement type 
as MNCN(I) and is preset to the same 
values 

- Main probe measurement noise code 
(IPMN = 0) 

= 0 - PMN array will be set equal to 
MNCN array 

= 1 - PMN array will be specified in 
the namelist 
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SMN( 12) 


- Array of miniprobe measurement noise 
variances for each type of measurement. 
SMN(I) refers to same measurement 
type as MNCN(I) and is preset to the 
same values 

ISMN - Miniprobe measurement noise code 

(ISMN = 0) 

= 0 - SMN array will be set equal to 
MNCN array 

= 1 - SMN array will be specified in 
the namelist 

See section 4 for definitions of NENT1 and NENT2. 

XEE(5) - Miniprobe release execution error 

variances 

XEE(l) : spin rate variance 
XEE(2) : boom length variance 
XEE(3) : spin axis right ascension 
variance 

XEE(4) : spin axis declination 
variance 

XEE(5) : release angle variance 

YYL - Miniprobe boom length 

TIMPCT - Approximate probe impact (trajectory) 

time. Can be obtained from a N0MNAL 
run 

RPS - Probe sphere radius; equals planet 

radius plus altitude of probe above 
planet surface at entry 

IUTC - Miniprobe targeting code (IUTC = 0) 

= 0 - Target controls computed inter- 
nally 

= 1 - Target controls supplied by user 

The next four variables are user-supplied miniprobe target controls 
and must be specified only if IUTC = 1. 

XPHI - Miniprobe release angle 
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ABW - Miniprobe spin rate magnitude 

ALFA - Right ascension of spin axis relative 

to the ecliptic coordinate system 

DELT - Declination of spin axis relative 

to the ecliptic coordinate system 

The remaining variables are required for internal miniprobe tar- 
geting and must be specified only if IUTC = 0. 

ACTPP - Accuracy level for virtual-mass propa- 

gation for miniprobe targeting 
(ACTPP = 2.5 x 10- 5 ) 

IPR0PI - Trajectory propagation code (1PR0PI = 1) 

= 1 - Conic propagation 
= 2 - Initial conic propagation with 

virtual-mass propagation refine- 
ment 

IPCSK - Code that specifies planetocentric 

coordinate system relative to which 
miniprobe targets are defined 
(IPCSK = 1) 

= 1 - Subsolar orbital-plane coordinate 
system 

= 2 - Equatorial coordinate system 

RATP(3) - Array of right ascensions of the 

three miniprobe targets relative to 
the coordinate system specified by 
IPCSK 

DCTP(3) - Array of declinations of the three 

miniprobe targets relative to the co- 
ordinate system specified by IPCSK 

ISA0 - Spin axis orientation flag (ISA0 = 1) 

= 1 - Spin axis declination and right 
ascension are free controls 
~ 2 - Spin axis is aligned with space- 
craft velocity vector at release 
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= 3 - Spin axis is perpendicular to 

the spacecraft-sun line, parallel 
to the ecliptic plane, and with- 
in 90° of the spacecraft velocity 
vector at release 

= 4 - Spin axis declination and right 
ascension are fixed by user- 
specified values of DCSAF and 
RASAF 

DCSAF - Fixed spin axis declination at release; 

specified only if ISA0 = 4 

RASAF - Fixed spin axis right ascension at re- 

lease; specified only if ISA0 = 4 

SO - Step size upper bound in the control 

space (SO =0.1) 

WFLS(3) - Miniprobe target site weighting fac- 

tors that indicate the relative impor- 
. tance of impacting each of the three 
target sites 


3.2.2 Measurement Schedules 

Three successive measurement schedules must appear immediately 
after the namelist ERRAN section. The first schedule appears on 
NENT cards and defines the measurement schedule for the primary 
vehicle. Immediately following these NENT cards are NENT1 cards 
that define the main probe measurement schedule. Following the 
NENT1 cards are NENT2 cards that define the single measurement 
schedule for all three miniprobes. 

Each card defines an entry in the measurement schedule ac- 
cording to the following format : 

From DAY1 (F10.0) to DAY 2 (F10.0), every X (F10.0) days, 
measurement code ITRK (110) . 

The measurement codes are defined as follows: 

ITRK = 1 - Range rate (idealized station) 

2 - Range and range rate (idealized 

station) 

3 - Range rate (station I) 
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= 4 - Range and range rate (station 1) 

5 - Range rate (station 2) 

6 - Range and range rate (station 2) 

7 - Range rate (station 3) 

8 - Range and range rate (station 3) 

9 - Three star-planet angles 

10 - Apparent planet diameter 

11 - Star-planet angle 1 

12 - Star-planet angle 2 

13 - Star-planet angle 3 

The total number of primary vehicle measurements must not exceed 
1000, and measurement times must not coincide. The total number 
of main probe measurements must not exceed 100. The total number 
of miniprobe measurements must not exceed 100. 


3.2.3 Namelist GENRAL 
GP(6,6) 

GPS (12, 12) 

GU ( 8 , 8) 

GV(15 , 15) 

GW(12,12) 

GCXXS(6 ,12) 

GCXU(6,8) 
GCXV(6 , 15) 

GCXW(6 , 12) 


- Actual spacecraft position/velocity 

covariance matrix P' (GP = P) 

o 

- Actual solve-for parameter covariance 

matrix P* (GPS = PS) 
s 

- Actual dynamic-consider parameter co- 

variance matrix V' (GU = VO) 
o 

- Actual measurement-cons ider parameter 

covariance matrix V' (GV = VO) 

o 

- Actual ignore parameter covariance 
matrix (GW = identity matrix) 

- Actual s tate/solve-f or parameter co- 

variance matrix C' (GCXXS = CXXS) 
xx 

s 

- Actual state/dynamic-consider param- 
eter covariance matrix C' (GCXU = CXU) 

xu 

- Actual state/measurement-consider 
parameter covariance matrix C’ 

vy 

(GCXV = CXV) 

- Actual state/ignore parameter covari- 
ance matrix C' (GCXW = 0) 

xw 
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GCXSU( 12 , 8) 

GCXSV( 12 , 15) 

GCXSW(12, 12) 
GCUV(8, 15) 
GCUW(8, 12) 
GCVW( 15 , 12) 

EXI(6) 

EXSI(12) 

EU(8) 

EV(15) 

EW(12) 


- Actual solve-for parameter/dynamic- 

consider parameter covariance matrix 

C' (GCXSU = CXSU) 
xu 
s 

- Actual solve-for parameter/measurement- 
consider parameter covariance matrix 

C' (GCXSV = CXSV) 
s 

- Actual solve-for parameter/ignore 
parameter covariance matrix C' 

(GCXSW = 0) x s W 

- Actual dynamic-consider parameter/ 

measurement-consider parameter co- 

variance matrix C' (GCUV = 0) 
uv 

- Actual dynamic-consider parameter/ 

ignore parameter covariance matrix 

C' (GCUW = 0) 
uw 

- Actual measurement-consider parameter/ 

ignore parameter covariance matrix C' 
(GCVW =0) W 


- Actual spacecraft position/velocity 

deviation mean x' (EXI = 0) 
o 

- Actual solve-for parameter deviation 
mean x' (EXSI = 0) 

o 

- Actual dynamic-consider parameter 
deviation mean u^ (EU = 0) 

- Actual measurement-consider parameter 

deviation mean v' (EV = 0) 
o 

- Actual ignore parameter deviation 
mean w ' (EW = 0) 
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IGDNF - Actual dynamic noise flag (IGDNF = IDNF) 

= 0 - Actual dynamic noise is zero 
= 1 - Actual dynamic noise is not zero 

GDNCN(3) - Array of constants used to calculate 

actual dynamic noise covariance ma- 
trix; must be specified if IGDNF equals 
1 

IGMNF - Actual measurement noise flag (IGMNF = 0) 

= 0 - Actual measurement noise is 
constant 

= 1 - Actual measurement noise is not 

constant (option is not available 
with this program) 

GMNCN(12) - Actual measurement noise variance for 

each type of measurement. GMNCN(I) 
refers to same measurement type as 
MNCN(I) (GMNCN - MNCN) 


EVK - Actual proportionality execution er- 

ror mean (EVK =0.) 

EVS - Actual resolution execution error 

mean (EVS =0.) 

EVA - Actual pointing angle alpha execution 

error mean (EVA = 0.) 

EVB - Actual pointing angle beta execution 

error mean (EVB =0.) 

VARK - Actual proportionality execution 

error variance (VARK = SIGPR0) 

VARS - Actual resolution execution error 

variance (VARS = SIGRES) 

VARA - Actual pointing angle alpha execution 

error variance (VARA = SIGALP) 

VARB - Actual pointing angle beta execution 

error variance (VARB = SIGBET) 
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3.3 SIMUL Input Description 

The input of the simulation program consists of: 

a. A card containing the variable 1RUNX (110 field) that 
indicates how many different runs are to be made and 
is read only once. 

b. A card containing the problem identification variable 
IPR0B (110 field) that precedes each set of input 
data. 

c. A simulation namelist section entitled SIMUL. 

d. Three successive measurement schedules for the primary 
vehicle, main probe, and miniprobes in that order. 

Namelist SIMUL contains all variables appearing in namelist 
ERRAN except for the variables ICDQ3 and IGEN. No generalized co- 
variance variables may appear in namelist SIMUL and no element of 
the IAUGIN array may be set to 3. The variable IGUID(2,I) may 
also take on the variable 1 for IGUID(1,I) = 1, 2, or 3. Variables 
appearing in namelist SIMUL, but not in namelist ERRAN, are defined 
below. 


1. Actual trajectory variables 

NB0D1 - Number of celestial bodies considered 

in the generation of the actual tra- 
jectory (NB0D1 = 11) 

NBl(ll) - Array of codes of celestial bodies 

considered in the generation of the 
actual trajectory; NB1 codes defined 
identically to NB codes (NB1 =1, 2, 
3, 4, 5, 6, 7, 8, 9, 10, 11) 


ACC1 - Actual trajectory accuracy figure 

(ACC1 - 1.0 x 10“ 6 ) 

ADEVX(6) - Actual initial position and velocity 

deviations from the most recent nominal 
trajectory (ADEVX =0, 0, 0, 0, 0, 0) 
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2. Actual dynamic biases 


DMUSB 

DMUPB 

DAB 

DEB 

DIB 

DN0B 

ro 

DMAB 

TTIM1 

TTIM2 

UNMAC(3 8 3) 


- Actual bias in the gravitational 
constant of the sun (DMUSB =0.) 

- Actual bias in the gravitational 
constant of the target planet 
(DMUPB =0.) 

- Actual bias in the semimajor axis of 
the target planet (DAB = 0.) 

- Actual bias in the eccentricity of 
the target planet (DEB =0.) 

- Actual bias in the inclination of the 
target planet (DIB =0.) 

- Actual bias in the longitude of the 
ascending node of the target planet 
(DN0B =0.) 

- Actual bias in the argument of periapsis 
of the target planet (DWB = 0.) 

- Actual bias in the mean anomaly of 
the target planet (DMAB = 0.) 

- The first time at which the values 
used for the actual unmodeled accel- 
eration will be altered 

(TTIM1 - 1. x 10 50 ) 

- The second time at which the values 
used for the actual unmodeled accel- 
eration will be altered 

(TTIM2 « 1. x 10 50 ) 

- Array of actual unmodeled accelerations 
to be used over a given time interval. 
Each row defines the values used over 

a given time interval; columns define 
the X, Y, and Z components of the un- 
modeled acceleration (UNMAC = array 
of zeros) 
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3. Actual measurement biases and variances 


BIA(12) 


SLB(9) 


- Array of actual biases for each type 
of measurement (BIA = array of zeros) . 
The following association exists: 

BIA(l) = Range bias (idealized station) 

(2) = Range-rate bias (idealized 

station) 

(3) = Range bias (station 1) 

(4) = Range-rate gias (station 1) 

(5) = Range bias (station 2) 

(6) = Range-rate bias (station 2) 

(7) = Range bias (station 3) 

(8) = Range-rate bias (station 3) 

(9) = Star-planet angle 1 bias 

(10) = Star-planet angle 2 bias 

(11) = Star-planet angle 3 bias 

(12) = Apparent planet diameter 

bias 

- Array of actual biases in the loca- 
tions of the three tracking stations 
on the rotating earth (SLB ■ array 

of zeros) . The following association 
exists : 

SLB(l) = Station 1 altitude bias 

(2) = Station 1 latitude bias 

(degrees north) 

(3) = Station 1 longitude bias 

(degrees east) 

(4) = Station 2 altitude bias 

(5) = Station 2 latitude bias 

(degrees north) 

(6) = Station 2 longitude bias 

(degrees east) 

(7) = Station 3 altitude bias 

(8) = Station 3 latitude bias 

(degrees north) 

(9) = Station 3 longitude bias 

(degrees east) 


2 7-20 



IAMNF 


AVABM(12) 


- Actual measurement noise code 
(IAMNF = 0) 

= 0 - Actual measurement noise has 

same statistics assumed by the 
navigation process (as represented 
by the MNCN array and the co- 
variance matrix) 

= 1 - Actual measurement noise has dif- 
ferent statistics from those as- 
sumed by the navigation process. 
These statistics are defined by 
the AVARM array 

- Array of actual variances for each 
type of measurement; specified only 
if IAMNF = 1. The following asso- 
ciation exists: 


AVARM ( 1) 
( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


= Range variance (idealized 
station) 

= Range-rate variance (idea- 
lized station) 

= Range variance (station 1) 

= Range-rate variance 
(station 1) 

= Range variance (station 2) 

= Range-rate variance 
(station 2) 

= Range variance (station 3) 

= Range-rate variance 
(station 3) 

= Star-planet angle 1 vari- 
ance 

= Star-planet angle 2 vari- 
ance 

= Star-planet angle 3 vari- 
ance 

= Apparent planet diameter 
variance 
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4. Actual midcourse velocity correction errors 


ARES (10) 


APR0(1O) 

AALP(10) 


ABET (10) 


- Array of actual resolution execution 
errors. The elements of this array 
and the following three arrays must 
correspond to the elements of the T3 
array and need be specified only if 
the T3 array has been specified 
(ARES = array of zeros) 

- Array of actual proportionality exe- 
cution errors (APR0 = array of zeros) 

- Array of actual pointing angle alpha 
execution errors (AALP = array of 
zeros) 

- Array of actual pointing angle beta 
execution errors (ABET = array of 
zeros) 


5. Quasi-linear filtering event variables 

NEV5 - Number of quasi-linear filtering 

events (NEV5 =0) 

T5(20) - Array of times at which quasi-linear 

filtering events occur; specified 
only if NEV5 is nonzero. Chronologi- 
cal order required 

6. Probe release event variables 


PAVARM(12) 


SAVARM(12) 


- Array of actual main probe measure- 
ment noise variances for each measure- 
ment type. PAVARM(I) refers to same 
measurement type as AVARM(I) and is 
preset to the same values 

- Array of actual miniprobe measurement 
noise variances for each measurement 
type. SAVARM(I) refers to same meas- 
urement type as AVARM(I) and is pre- 
set to the same values 
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PBIA( 12) 


- Array of actual main probe biases 
for each measurement type. PBIA(I) 
refers to same measurement type as 
BIA(l) and is preset to the same 
values . 

SBIA(12) - Array of actual miniprobe biases 

for each measurement type. SBIA(I) 
refers to same measurement type as 
BIA(I) and is preset to the same 
values 

The next five variables define the actual miniprobe release execu- 
tion errors: 

DW - Actual spin-rate execution error 

(DW = 0.) 

DA - Actual spin axis right ascension 

execution error (DA = 0.) 

DD - Actual spin axis declination execu- 

tion error (DD = 0.) 

DP - Actual release angle execution error 

(DP = 0.) 

DL - Actual boom length error (DL = 0.) 


NQLE(2) - Number of probe quasi-linear filtering 

events. NQLE(l) = number of main 
probe quasi-linear filtering events. 
NQLE(2) = number of miniprobe quasi- 
linear filtering events 

QLTIM(2,10) - Array of probe quasi-linear filtering 

event times. QLTIM(1,I) defines the 
sequence of event times for the main 
probe. QLTIM(2,I) defines the se- 
quence of event times for the mini- 
probe. Event times must be arranged 
in chronological order 
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3,3 SIMUL Input Description 


The input of the simulation program consists of (a) a card containing 
the variable IRUNX (110 field) which indicates how many different runs are 
to be made and is read only once, (b) a card containing the problem identi- 
fication variable IPROB (110 field) which precedes each set of input data, 

(c) a namelist entitled SIMUL, and (d) a measurement schedule. The measure- 
ment schedule is treated exactly the same way as it is treated in the error 
analysis program. Namelist SIMUL contains all variables appearing in 
name list ERRAN except for the variable ICDQ3 . In addition, IGUID(2,I) 
can alsc take on the value 1 for IGUID(1,I) * 1,2, or 3. 

Variables appearing in namelist SIMUL, but not in namelist ERRAN, 

are defined below. 

a. Actual trajectory variables 

NBODl - Number of celestial bodies considered in the 

generation of the actual trajectory. (NBODl * 11) 

NBl(ll) - Array of codes of celestial bodies considered in 

the generation of the actual trajectory; NBl codes 
defined identically to NB codes. (NBl •» 1, 2, 3, 

4, 5, 6, 7, 8, 9, 10, 11) 

ACC1 - Actual trajectory accuracy figure. (ACCl = 1.0 x 

10 " 6 ) 

ADEVX(6) - Actual initial position and velocity deviations 

from the most recent nominal trajectory. (ADEVX = 

0, 0, 0, 0, 0, 0) 

b. Actual dynamic biases 

DMUSB - Actual bias in the gravitational constant of the 

Sun. (DMUSB =0.) 

DMUPB - Actual bias in the gravitational constant of the 

target planet. (DMUPB =0.) 

DAB - Actual bias in the semi -major axis of the target 

planet. (DAB = 0.) 

DEB - Actual bias in the eccentricity of the target planet. 

(DEB - 0.) 

DIB - Actual bias in the inclination of the target planet. 

(DIB - 0.) 
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DN0B - Actual bias in the longitude of the ascending 

node of the target planet. (DN0B =0.) 

DWB - Actual bias in the argument of periapsis of 

the target planet. (DWB =0.) 

DMAB - Actual bias in the mean anomaly of the target 

planet. (DMAB = 0.) 

TTIMl - The first time at which the values used for the 

actual unmodeled acceleration will be altered. 
(TTIMl = 1. x 10 50 ) 

TTIM2 - The second time at which the values used for the 

actual unmodeled acceleration will be altered. 
(TTIM2 = 1. x 10 50 ) 

UNMAC(3,3) - Array of actual unmodeled accelerations to be used 
over a given time interval. Each row defines the 
values used over a given time interval; columns 
define the X, Y, and Z components of the unmodeled 
acceleration. (UNMAC = array of zeros) 

c. Actual measurement biases and variances 


BIA(12) - Array of actual biases for each type of measurement. 

(BIA = array of zeros) The following association 
exists : 


BIA(l) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 
( 11 ) 
( 12 ) 


Range bias (idealized station) 
Range-rate bias (idealized station) 
Range bias (station 1) 

Range-rate bias (station 1) 

Range bias (station 2) 

Range-rate bias (station 2) 

Range bias (station 3) 

Range-rate bias (station 3) 
Star-planet angle 1 bias 
Star-planet angle 2 bias 
Star-planet angle 3 bias 
Apparent planet diameter bias 


SLB(9) - Array of actual biases in the locations of the 

three tracking stations on the rotating earth. 
(SLB = array of zeros) The following association 
exists : 
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IAMNF 


AVARM (12) 


Actual midcourse 
ARES (10) 


APRO (10) 


SLB(l) = Station 1 altitude bias 

(2) = Station 1 latitude bias (degrees north) 

(3) = Station 1 longitude bias (degrees east) 

(4) = Station 2 altitude bias 

(5) = Station 2 latitude bias (degrees north) 

(6) = Station 2 longitude bias (degrees east) 

(7) = Station 3 altitude bias 

(8) = Station 3 latitude bias (degrees north) 

(9) = Station 3 longitude bias (degrees east) 

- Actual measurement noise code. (IAMNF = 0) 


= 0 - Actual measurement noise has same statistics 
assumed by the navigation process (as repre- 
sented by the MNCN array and the R^ covariance 
matrix) . 

= 1 - Actual measurement noise has different statistics 
from those assumed by the navigation process. 
These statistics are defined by the AVARM 
array. 

Array of actual variances for each type of measure- 
ment; specified only if IAMNF = 1. The following 
association exists: 


AVARM ( 1 ) 
( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 
(ID 
( 12 ) 


Range variance (idealized station) 
Range-rate variance (idealized station) 
Range variance (station 1) 

Range-rate variance (station 1) 

Range variance (station 2) 

Range-rate variance (station 2) 

Range variance (station 3) 

Range-rate variance (station 3) 
Star-planet angle 1 variance 
Star-planet angle 2 variance 
Star-planet angle 3 variance 
Apparent planet diameter variance 


velocity correction errors 

Array of actual resolution errors. The elements 
of this array and the following three arrays 
must correspond to the elements of the T3 array 
and need be specified only if the T3 array has 
been specified. (ARES = array of zeros) 

Array of actual proportionality errors. (APRO = 
array of zeros) 



AALP(IO) 


- Array of actual pointing angle alpha errors. 

(AALP = array of zeros) 

ABET(IO) - Array of actual pointing angle beta errors. 

(ABET = array of zeros) 

e. Quasi-linear filtering event variables 

NEV5 - Number of quasi-linear filtering events. (NEV5 = 0) 

T5(20) - Array of times at which quasi-linear filtering 

events occur; specified only if NEV5 is non-zero. 
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4. OUTPUT DESCRIPTION 


4.1 NOMNAL Output Description 

The output from the nominal trajectory generator NOMNAL is 
conveniently divided into seven sections — initial data, inter- 
planetary zero iterate, lunar zero iterate generation, nominal 
trajectory data, nonlinear guidance, orbit insertion, execution of 
a correction, and miniprobe targeting. Each of these sections will 
be discussed individually. 

a. Initial Data 


The namelist TARIN is first printed out as the data are input. 
Then much of the data with useful extensions are recorded In 
unified sections. More specific definitions of the input data 
are given in the Input Description Section. The trajectory 
parameters are first recorded. 


ALNGTH 


Length units per AU for trajectory output 


TM 

NBOD 

NB 

IB ARY 

ICOORD 

NCPR 

TMPR 


Time units per day for trajectory output 

Number of bodies used in integration 

Array of codes of gravitational bodies 

Flag indicating whether inertial system 
is heliocentric ecliptic (IBARY = 0) or 
barycentric ecliptic (IBARY = 1) 

Flag not currently used 

Integration increments between printouts 
of nominal trajectory 

Days between printouts of nominal tra- 
jectory 


Next the zero iterate data are printed, including the initial 
state itself and the parameters necessary to compute it. 

IZERO Flag designating which option of zero iter- 

ate generation was used 

ZDAT(6) Zero iterate state (inertial ecliptic) 

generated 
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RP 

FI 

A1 

A2 

T1 

T2 

LAT 

LON 

THD 

RAT 

AZI 


Parking orbit radius 
Injection true anomaly 
Angle of first bum 
Angle of second burn 
Time interval of first burn 
Time interval of second burn 
Latitude of launch site 
Longitude of launch site 
Rotation rate of launch planet 
Inverse rate in parking orbit 
Launch azimuth 


Then for each guidance event to be processed (KTIM i- 0 ) , the 
following data are recorded. 


EVENT TYPE 


REF TIME 
REF CODE 


Value of KTYP flag designating whether 
event is targeting, retargeting, orbit 
insertion, main-probe propagation, mini- 
probe targeting, or termination 

Value of TIMG, giving time interval be- 
tween event and reference epoch 

Value of KTIM, designating the epoch to 
which event times are referenced (initial 
time, SOI, CA) 


EVENT (TARGET) The calendar date of the event (target 
CALENDAR DATE time) 


EVENT (TARGET) The Julian date (epoch 1900) of the gui- 
JULIAN DATE dance event (target time) 


EVENT (TARGET) The time of the guidance event (target 
TRAJ DAY time) referenced to the initial time 



CONTROL 


IMP 


MOD 

TAR KEY 

TAR1, TAR2, 
TAR3 


Code indicating which of the two possible 
sets of velocity control variables are to 
be used 

Implementation code giving the compute- 
execute mode KMXQ 

Model to be used in execution (MDL) 

Target parameter keys (KTAR) 

Desired values of target parameters (TAR) 


TOLl, T0L2, 
T0L3 


Allowable tolerances of target parameters 
(TOL) 


DVX, DVY , DVZ Velocity components (inertial ecliptic) 

of correction supplied externally (KMXQ = 2) 

MAT, BAD, Values of targeting matrix option (MAT), 

ITS, BIT bad-step option (IBADS) , iterations (NOIT) 

and bad iterations (MAXB) 


Interplanetary Zero Iterate 

The massless planet (or point-to-point) approximation is used 
as a zero iterate for interplanetary trajectories. Since 
these data are printed in SPARC, it is convenient to record 
the zero iterate information in the same format as SPARC. The 
launch date, arrival date, and flight time are self explanatory. 
The heliocentric parameters are listed. 


RL, RP 

Heliocentric radius of initial and final 
points (10 6 km) 

LAL, LAP 

Heliocentric ecliptic latitude of initial 
and final points 

LOL, LOP 

Heliocentric ecliptic longitude of ini- 
tial and final points 

VL, VP 

Heliocentric speeds at initial and final 
points (km/s) 

GAL, GAP 

Flightpath angles at initial and final 
points 



AZL, 

AZP 

Azimuth at initial and final points 

HCA 


Heliocentric central angle of transfer 

TAL, 

TAP 

True anomaly at initial and final points 

SMA 


Semimajor axis of heliocentric conic 
(10 6 km) 

ECC 


Eccentricity of heliocentric conic 

INC 


Heliocentric ecliptic inclination of conic 

RCA 


Perihelion of conic (10 6 km) 

APO 


Aphelion of conic (10 6 km) 

VI 


Heliocentric speed of launch planet (if 
applicable) 

V2 


Heliocentric speed of target planet (if 
applicable) 


The data defining the launch planet conic are then listed: 


C3 

The launch energy (= VHL 2 ) 

VHL 

The hyperbolic excess velocity (km/s) 

DLA 

The declination of the departure asymptote 

RAL 

The right ascension of the departure 
asymptote 

RAD 

The injection radius with respect to the 
launch planet (km) 

VEL 

The injection velocity with respect to the 
launch planet (km/s) 

PTH 

The injection path angle 

VHP 

The hyperbolic excess velocity at the 
target planet 

DPA 

The declination of the approach asymptote 
(heliocentric ecliptic) 



RAP 


The right ascension of the approach asymptote 
(heliocentric ecliptic) 


ECC The eccentricity of the launch planet 

conic 

LNCH AZMTH The launch azimuth 

LNCH TIME The time of launch on the launch date 

L-I TIME The time between launch and injection 

(seconds) 

INJ LAT The injection latitude 

INJ LONG The injection longitude 

INJ RT ASC The right ascension at injection 

INJ TIME The time at injection on the launch date 

PO CST TIM The parking orbit coast time (sec) 

c. Lunar Zero Iterate Generation 

The zero iterate generation for lunar trajectories proceeds 
in two stages. The first stage determines iteratively a tar- 
geted patched conic; the second stage generates a targeted 
multi-conic trajectory. The following information is recorded 
for each iterate and for each perturbed trajectory in the 
patched conic trajectory: 

ITR Iteration counter 

ALPHA, DELTA, Controls used on current trajectory (see 

THETA LUNCON, LUNTAR) 

SIGMA Launch azimuth 

SMA, RCA, INC Semimajor axis, radius of closest ap- 
proach, and equatorial inclination of 
current iterate 

VSI (1, 2, 3) Hyperbolic approach velocity of current 
iterate 

B*T, B • R Impact plane parameters (earth equatorial) 

of current iterate 
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Following the perturbed trajectories a summary of the targeting 
data for that iterate is listed 


PHI MATRIX 

Sensitivity matrix from numerical differ- 
encing 

PHI INVERSE 

Targeting matrix (inverse of sensitivity 
matrix) 

ERRORS 

The current errors in semimajor axis, B*T 
and B*R from the desired values 

TARGETS 

The desired values of semimajor axis, 
B • T and B* R 

CORRECTION 

The correction to be added to the controls 
for the next iterate 

In the initial 

targeting to semimajor axis the numerical partial 


of semimajor axis to the alpha control is recorded as PARTA. 


During the multi-conic stage of targeting slightly different 
information is recorded. Initially the targeting scheme is 
described: 


TARGETS Desired values of semimajor axis, incli- 

nation, radius of closest approach, and 
time of closest approach 


TOLERANCES 

PERTURBATIONS 


MAX STEPS 

MAX ITERS 

MULTI-CONIC 

STEP 


Acceptable tolerances of semimajor axis, 
B*T, B*R, and t 

v>A 

Perturbations used in constructing the 
first targeting matrix (later iterations 
use a perturbation which for each com- 
ponent would null the time error if the 
previous sensitivity matrix were still 
valid) 

Maximum correction allowed on first iter- 
ate (later iterations use 100 times the 
perturbation size) 

Maximum number of iterations allowed in 
multi-conic targeting 

Step size in hours of multi-conic propa- 
gator 
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For each iterate and for each perturbed trajectory the fol- 
lowing data is recorded: 

ITER Iteration counter 


JULIAN DATE Julian date of injection (referenced to 

1950) 


x, y, z, vx, Injection state in earth ecliptic coor- 
vy, vz dinates 


SMA, B*T, B*R, Target parameter values achieved on cur- 
TCA rent trajectory (impact plane parameters 

in earth ecliptic coordinates) 


Following the perturbed trajectories a summary of the current 
trajectory data is given. 


SENSITIVITY Sensitivity matrix computed from numerical 

MATRIX partials 


TARGETING Targeting matrix (inverse of sensitivity 

MATRIX matrix) 


ERRORS 

TARGETS 

PREDICT 


Errors in target parameters (a, B*T, B*R, 
t CA^ 

Target values of (a, B*T, B* R, t„.) 
Predicted corrections 


CORRECT 


Actual correction added to controls after 
applying constraints 


After obtaining an acceptable trajectory, a summary is given 
listing the injection Julian and calendar date and the injec- 
tion state in both earth-centered and barycenter ecliptic co- 
ordinates . 


Nominal Trajectory Data 

The nominal trajectory can consist of up to five branches if 
the spacecraft is a conglomerate vehicle made up of a bus, main 
probe, and three miniprobes. The first branch corresponds to 
the bus trajectory, the second to the main probe, and the last 
three to the miniprobe's. The propagation and concurrent print- 
out of the bus trajectory between guidance events is directed 



by the subroutine TRJTRY. The virtual-mass propagation (VMP) 
initiation data are provided, together with the trajectory 
status information at intervals of NCPR integration steps and 
TMPR days. 

A time-history of the main-probe trajectory is initiated by 
scheduling a main-probe propagation event. When the trajectory 
time reaches that at which the propagation event is scheduled, 
control is transferred to MPPROP. This subroutine then calls 
VMP to propagate and simultaneously record the trajectory 
of the main probe from its current state, which is identical 
to that of the bus, to its final stopping condition. The re- 
sulting trajectory history is preceded for identification pur- 
poses by the title "Main-Probe Propagation Event" followed 
by the heading "Main-Probe Approach Trajectory." It contains 
the usual VMP propagation-initiation and trajectory status data 
at fixed intervals of 100 integration steps and 5 days. 

Finally histories of the three miniprobe trajectories for the 
minimum-miss release controls are provided without user sched- 
uling. At completion of the miniprobe targeting, the subroutine 
TPRTRG directs VMP via TPPROP to provide the usual VMP time 
history for each of the miniprobe trajectories in succession, 
starting from the release state. The trajectory corresponding 
to the ith miniprobe is identified by the title "Miniprobe I 
Minimum-Miss Approach Trajectory." The intervals between tra- 
jectory state printouts are the same as for the main-probe his- 
tory. 

The standard trajectory information printed during a given call 
to the VMP, with the print flag set, is as follows: 

1) Initiation data (provided only at the start of the 

propagation) 

a) The output length units per AU and time units per 
day, 

b) The number of gravitational bodies, 

c) The initial date of integration, 

d) The final date of integration, 

e) The initial trajectory time, 

f) The integration accuracy level. 
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g) The true anomaly increment in radians correspond- 
ing to the accuracy level; 

2) Trajectory status data (provided at intervals of NCPR 


integration steps and TMPR days) , 


a) Block 1 - Spacecraft inertial state. 


(1) 

The current trajectory time, 


(2) 

The cumulative number of integration incre- 
ments used to this point. 

(3) 

The current spacecraft state (position 
radius, velocity vector, and speed) in 
(heliocentric or barycentric) ecliptic 
dinates. 

vector , 

inertial 

coor- 


b) Block 2 - Gravitational bodies states, 

(1) The current calendar date, 

(2) The current Julian date (reference year 0) , 

(3) The current state of all gravitating bodies 
in inertial ecliptic coordinates, 

c) Block 3 - Spacecraft relative state, 

(1) The current state of the spacecraft relative 
to all gravitational bodies in inertial 
ecliptic coordinates, 

d) Block 4 - Virtual mass data, 

(1) The current state of the virtual mass in in- 
ertial ecliptic coordinates, 

(2) The current state of the spacecraft relative 
to the virtual mass in inertial ecliptic co- 
ordinates, 

(3) The Kepler (angular momentum) vector in iner- 
tial ecliptic coordinates, 

(4) The eccentricity vector in inertial ecliptic 
coordinates , 



(5) The virtual mass magnitude and magnitude rate, 

e) Block 5 - Virtual mass relative positions, which 
are the current position vector and radius of the 
virtual mass relative to all gravitational bodies, 


f) Block 6 - Navigation parameters. 



(1) 

Flightpath angle, 


(2) 

Angle between relative velocity and plane of 
the sky , 


(3) 

Geocentric declination. 


(A) 

Earth/spacecraft/target planet angle. 


(5) 

Antenna axis/ear th angle, 


(6) 

Antenna axis/limb of sun angle, 


(7) 

Spacecraft occultation ratio for the sun and 
gravitational bodies. 

Nonlinear Guidance 

At a targeting or retargeting guidance event, the following 
general information is first printed: 

1) 

The trajectory time, calendar date, and Julian date 
(reference 1900) of the current event; 


2) The event codes of the current event including its 
index KUR, the type KTYP, the compute/execute flag 
KMXQ, and the model MDL; 

3) The current spacecraft state in inertial ecliptic 
and all-gravitational body-centered ecliptic coordin- 
ates ; 

A) A listing of the targeting parameter key definitions; 

5) The targeting specifications, including the target 
parameter keys, desired values, and tolerances; 

6) The targeting scheme parameters, including the accuracy 
levels to be used, the maximum velocity movement al- 
lowed, and the bad-step flag. 



For each trajectory generated during the course of a nonlinear 
guidance event the following data are recorded. 


ACCURACY 

The integration accuracy level used on 
the integration 

VX, VY, 

VZ 

The velocity components (inertial eclip- 
tic) used on the integration 

TAR(i) , 
TAR(k) 

TAR(j) , 

The target parameter codes (i,j,k) and 
the corresponding values achieved on the 
integration 

AUX(Jt) , 
AUX(n) 

AUX(m) , 

The auxiliary parameter codes (£,m,n) and 
the corresponding values achieved on the 
integration 

INCR 


The number of integration increments in 


the integration 


If the method used to target the trajectory is the Newton- 
Raphson technique, the nominal trajectory and the three per- 
turbed trajectories used to generate the targeting matrix are 
recorded in the above format. If the velocity controls used 
are those relative to the current spacecraft launch-plane to- 
centric velocity rather than to the heliocentric velocity, 
additional printout is supplied from the subroutine KTROL as 
follows . 

STATE Launch-planetocentric state of spacecraft 

in km and km/s 

IOPT Code signaling to KTROL which of the ve- 

locity controls C^, C2* and are to be 

perturbed 

CON Velocity control components. First com- 

ponent has units of km/s; the last two 
components have units of radians 

DV Control correction to the spacecraft launch 

planetocentric velocity vector in km/s 



The following targeting information is printed after the last 
perturbed trajectory data of the current iteration and before 
the nominal data of the next. 


SENSITIVITY 

MATRIX 


TARGETING 

MATRIX 


AUX ERROR 


VEL COR 


DES AUX VAL 
DES TAR VAL 
TAR TOL 


The Jacobian matrix giving the sensitivi- 
ties of the auxiliary targets to the ve- 
locity controls 

The inverse of the Jacobian matrix giving 
the sensitivities of the velocity controls 
to the auxiliary targets 

The negatives of the errors of the current 
auxiliary target values from the desired 
values 

The predicted ecliptic cartesian velocity 
correction (after the DVMAX constraint) 
in km/s 

The desired auxiliary target values 
The desired actual target values 
The actual target tolerances 


Then the four trajectories of the next iterate are given and 
so on through the targeting. 

If the targeting method is the steepest descent technique, 
five trajectories are printed out for each iterate. The first 
trajectory is the nominal or iterate trajectory. The next 
three trajectories are the perturbed trajectories used to 
compute the gradient vector. If the velocity controls are 
those relative to the launch-planetocentric state, the same 
additional information is printed by the subroutine KTROL with 
each trajectory as was described for the Newton-Raphson tar- 
geting scheme. The following data unique to the descent tar- 
geting technique are printed out for each iterate. 

WEIGHTS The weighting factors for the auxiliary 

parameters for the scalar error 

PER-ERRS The errors corresponding to each of the 

perturbed trajectories 

The gradient based on those pertu-rbed 
trajectories 


GRADIENT 



CON -GRAD 

The conjugate gradient 

DIRECTION 

The unit vector in the direction of 
the correction 

EN 

The error of the nominal trajectory 

DD 

The directional derivative of the error 
in the correction direction 

HB 

The linearly predicted step size 

EB 

The error corresponding to the linearly 
predicted step size 

HS 

The predicted optimal step size 

HH 

The step size to be actually used 

DELTAV 

The actual correction to be made on the 
next iterate 

AUX ERR 

The auxiliary errors of the current 
iterate 

D-AUX 

The desired auxiliary parameter values 

D-TAR 

The desired target parameter values 

following diagnostic messages may be printed out during 
targeting. 

ENTER OUTER 
TARGETING 

(The current iterate missed the target 
planet SOI) 

RCA 

The radius of closest approach to the 
target planet 

ARTIFICIAL SOI 

The SOI set up for outer targeting 
(= 1.2 RCA) 

TCA 

The time of closest approach to the tar 
get planet 

ORIG TSI 

The original target time at the SOI 
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MODIFIED TSI 


The adjusted target time to hit the 
artificial SOI 



EXIT OUTER (The outer targeting has been success- 

TARGETING fully accomplished; return and target to 

the original values) 

BAD STEP - (Since the current iterate missed the SOI 

ITERATE when previous iterates intersected it, 

MISSED SOI reduce the correction) 

BAD STEP (The current iterate has a larger error 

than the previous) 

REDUCTION • The fractional reduction made in the ve- 

locity correction for a bad step 

PREVIOUS ERROR The previous error 

CURRENT ERROR The current error 

f. Orbit Insertion 

At the start of an orbit insertion event, the following general 
information is given: 

1) The trajectory time, calendar date, and Julian date 
(referenced 1900) of the current event; 

2) The event codes of the current event including its 
index KUR, the type KTYP, the compute/execute flag 
KMXQ, and the model MDL: 

3) The current spacecraft state in inertial ecliptic and 
all-gravitational body-centered ecliptic coordinates; 

4) The orbit insertion option selected (coplanar or non- 
planar) and the desired values of the target parameters 
(semimajor axis, eccentricity, and periapsis shift 

for coplanar, equatorial conic elements for the non- 
planar option. 

The program then records a detailed description of all impor- 
tant states generated during the analysis of the orbit in- 
sertion. In this section all states refer to target planet- 
centered equatorial coordinates. The Cartesian state (X, Y, 

Z, R, VX, VY, VZ, V) provides the position vector, radius, 
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velocity vector, and speed of the spacecraft with respect to 
the target body at the given epoch. The conic parameters (A, 

E, W, TA, I, N, RP, RA, TIME) supply the semimajor axis, ec- 
centricity, argument of periapsis, true anomaly, inclination, 
longitude of the ascending node, periapsis radius, apoapsls 
radius, and time from periapsis of the target planet-centered 
conic at the given epoch. The terms used in the summary of 
both coplanar and nonplanar insertions are defined as follows: 

1) The "decision state" is the Cartesian and conic state 
(on the approach hyperbola) of the spacecraft at the 
time at which the insertion event is computed; 

2) The "target orbit" is the unmodified desired orbit. 
Obviously only the conic parameters may be given; 

3) The "preinsertion state" is the predicted Cartesian 
and conic state on the hyperbola at the instant be- 
fore a candidate impulsive correction; 

4) The "postinsertion state" is the predicted Cartesian 
and conic state on the ellipse immediately following 
the candidate impulsive correction; 

5) The "insertion velocity" is the impulsive velocity 
(equatorial) of a candidate solution; 

6) The "errors" are the weighted scalar loss functions 
associated with each candidate solution; 

7) The "selected correction" is the impulsive velocity 
(equatorial) having the minimum loss function and 
therefore chosen for execution. 

In the coplanar insertion, the target orbit automatically lies 
in the plane of the approach asymptote. The shape of the tar- 
get ellipse is determined by the desired values of semimajor 
axis and eccentricity. Its orientation is fixed by the de- 
sired periapsis shift from the approach hyperbola. There are 
two possibilities for solutions: either the hyperbola and el- 

lipse intersect or they do not. In the former case there are 
two candidate solutions and the one with minimum delta velocity 
is chosen for execution. In the latter case three modifications 
of the target orbit leading to tangential solutions are analy- 
zed: (1) vary periapsis while holding apoapsls constant, (2) vary 

apoapsis holding periapsis constant, and (3) vary semimajor 
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axis while holding eccentricity constant. The weighting fac- 
tors used are identical to the indices. When a candidate modi- 
fication cannot lead to a tangential solution, a message to 
that effect is printed out. 

In the nonplanar case the two points of intersection of the 
approach hyperbola and the desired plane are computed and 
recorded. If one of the points lies in the impossible region, 
that information is recorded. The candidate modifications of 
the target orbit discussed in the coplanar case are made to 
determine the optimal impulsive correction. The weights are 
the same as above with the extension that all corrections on 
the departure ray are doubled. Invalid candidate modifica- 
tions are recorded. 

g. Execution of Correction 

At an execution event, the following information is recorded 
regardless of the model used. 

1) The title "Execution Event" is written; 

2) The impulsive velocity correction and its magnitude 
are recorded in inertial ecliptic coordinates; 

3) The conic elements (semimajor axis, eccentricity, 
argument of periapsis, inclination, longitude of 
the ascending node, and true anomaly) of the space- 
craft with respect to the dominant body in ecliptic 
coordinates before and after an impulsive addition 
of the correction is recorded; 

4) If the dominant body is not the sun, the same infor- 
mation is listed in dominant body equatorial coordi- 
nates . 

If the model used is the multiple thrusting arc, the following 
data are also written: 

1) The thrust magnitude, the nominal mass of the space- 
craft during the pulsing arc, the duration of a single 
pulse, and the resulting magnitude of the velocity in- 
crement are recorded; 

2) The ecliptic components of the total velocity incre- 
ment to be imparted, the nominal impulse of the series, 
and the final (partial) pulse of the series are listed; 
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3) The time information of the pulsing arc, including the 
calendar and Julian dates of the initiation, midpoint, 
and termination of the arc, are provided; 

4) The terms of the f and g series used to propagate the 
launch and target bodies through the pulsing arc are 
given; 

5) The pulse-by-pulse listing of the inertial ecliptic 
state following each of the pulses is given (propa- 
gation between pulses is perturbed conic); 

6) The state computed by adding the correction impul- 
sively at the nominal time of the correction and then 
propagating by the perturbed conic to the time at the 
end of the arc is provided for comparison; 

7) The two final states of (5) and (6) are converted to 
conic elements for a final comparison of the two tech- 
niques . 

h. Miniprobe Targeting 

When a miniprobe targeting event is begun, the following 
general information is given: 

1) The trajectory time, calendar date, and Julian date 
(epoch 1900) of the current event; 

2) The event codes for the current event, including its 
index KUR, type KTYP, compute/execute flag KMXQ, and 
model MDL; 

3) The current spacecraft (bus) state in heliocentric- 
ecliptic and gravitational-body-centered coordinates 
for all gravitating bodies (units are km and km/s) ; 

4) Planetocentric ecliptic state of spacecraft (bus) at 
impact based on an n-body propagation from the release 
state (units are km and km/s) ; 

5) Equivalent conic planetocentric ecliptic state of 
spacecraft (bus) at release in km and km/s. 
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Next, information from the Gauss least-squares routine is printed. 
This begins with a summary of procedure parameters. 

Gauss Least-Squares Parameters 


N 

M 

DELTA 


Cl 


C2 


EPS 


SO 


ITLIM 


Number of control variables (preset at 4 
for miniprobe targeting) 

Number of constraint variables (preset at 
6 for miniprobe targeting) 

Perturbation size used for all control 
variables in approximating the Jacobian 
matrix by divided differences (preset at 
10 -J for miniprobe targeting) 

Weighting factor applied to length change 
of control vector in convergence criterion 
(preset at 10 1 * for miniprobe targeting) 

Weighting factor applied to change in 
magnitude of miss index in convergence 
criterion (preset at 1 for miniprobe tar- 
geting) 

Upper limit on weighted sum of change in 
length of control vector and change in 
magnitude of miss index in the convergence 
criterion (preset at 1 for miniprobe tar- 
geting) 

Upper bound on length of control correction 
above when a steepest-descent rather than 
a pseudo inverse step is taken (input 
by user for miniprobe targeting) 

Maximum number of iterations permitted 
before the least-squares procedure is ter- 
minated (preset at 20 for miniprobe tar- 
geting) 


Then a series of descriptions of the individual iterations 
is printed. 
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Individual Iteration Description 


XN1 Nth iterate value of first control vari- 

able — release roll angle of the first 
miniprobe in radians 

XN2 Nth iterate value of the second control 

variable — tangential release velocity 
in decameters/s 

XN3 Nth iterate value of the third control 

variable — inertial ecliptic declination 
of spacecraft spin axis at release in 
radians 

XN4 Nth iterate value of the fourth control 

variable — inertial ecliptic right as- 
cension of spacecraft spin at release in 
radians 

GRADN Magnitude of the gradient of the miss 

index at the nth iterate 

YN Value of miss index at the nth iterate 

in km 2 

PHI(j) Jth component of constraint (miss) vector 

printed at the iteration point itself as 
well as at those points obtained from it 
by successively perturbing each compo- 
nent by DELTA (units are km) 

JACOBIAN Jacobian matrix of constraint vector with 

MATRIX respect to control vector as obtained by 

divided differencing 

PROJECTION Pseudoinverse of Jacobian matrix, i.e., 

MATRIX ( jT J} -i / 

LAMBDA Initial estimate of the step size in the 

search direction necessary to bracket the 
minimum of the miss index as required in 
the descent procedure 

ALPHA Fraction of the bracketing step at which 

miss index is evaluated for cubic inter- 
polation 
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YP(O) Slope of miss index in the search direction 

at the current iterate 

XK Length of a trial step in the search di- 

rection as used in bracketing the minimum 
of the miss index and fitting it with a 
. cubic polynomial 

Y(XK) Value of the miss index at the trial step 

XK 

The following messages may be printed out when the iteration 
process is terminated. A brief explanation of each is pro- 
vided. 

Iteration Termination Messages 

ADEQUATE CONVERGENCE OCCURRED ON PREVIOUS STEP - This 
is the normal termination after a successful miniprobe 
targeting 

CONVERGENCE DID NOT OCCUR - After ITLIM iterations the 
current and previous iterations are still too far apart 
to satisfy the convergence criterion 

PERFORMANCE INDEX DECREASES MONOTONICALLY IN SEARCH 
DIRECTION. SEARCH HAS BEEN TERMINATED - A minimum 
of the miss index in the direction of search could 
not be bracketed in 10 trial steps 

To conclude the least-squares printout, an "Iteration History" 
is supplied. It contains the control vector, the miss index 
and the gradient magnitude of the miss index for each iterate 
in sequential order. 

The "Minimum-Miss Release Controls for Conic Propagation" are 
next printed after being converted from the unusual dimensions 
required by the iteration process to the standard ones indi- 
cated in the output itself. 

Then the "Conic-Model Probe Impact Data" are printed. For 
each miniprobe and the bus (numbered miniprobe 0) the impact 
declination and right ascension in the probe-sphere system, 
the impact speed, and the impact flightpath angle are recorded. 
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The information for the miniprobes is based on the conic model 
while that for the bus is derived from the initial n-body pro- 
pagation of the bus from the release state to impact. In ad- 
dition to the above data, an angle of attack is printed for 
each miniprobe assuming its longitudinal body axis at impact 
remains parallel to the spacecraft spin axis at release. 

If n-body release controls are not requested by the user, the 
n-body time histories of the respective miniprobe trajectories 
for minimum-miss conic release controls are printed next as 
described in the nominal trajectory section. 

If n-body controls are requested, the output from the applica- 
tion of the least-squares routine to the n-body propagation 
model is printed next. Its format is identical to that for 
the conic model. After the miss-minimizing procedure is com- 
pleted, the "Minimum-Miss Release Controls for the N-Body 
Propagation" are printed. They too are identical in format 
to the conic case. Then the respective n-body miniprobe tra- 
jectory histories are printed for the minimum-miss n-body re- 
lease controls as described in the nominal trajectory section. 

Finally the "N-Body Model Probe Impact Data" are printed. 

The same quantities are presented as in the corresponding 
section for the conic model described above, but they are now 
calculated from the recently computed n-body miniprobe trajec- 
tories . 
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4.2 ERRAN Output Description 


The printed output of the error analysis mode is described in this 
section according to the following groups: input data, measurement output, 

additional trajectory output, eigenvector event output, prediction event 
output, guidance event output, and summary output. 

a. Input data 

The initial output consists of the following input data: 

(1) Namelist ERRAN. 

( 2 ) Calendar date and Julian date at launch. 

( 3 ) 1'inal calendar date and Julian date. 

(4) Initial trajectory time in days (TRTM1). 

(5) Lists of solve-for, dynamic consider, and measurement consider 
parameters augmented to the position/velocity state vector. 
Definitions of names appearing in this list are given below: 

RADIUS 1 Radius error of station 1 

LAT 1 Latitude error of station 1 

LONG 1 Longitude error of station 1 

RADIUS 2 Radius error of station 2 

LAT 2 Latitude error of station 2 

LONG 2 Longitude error of station 2 

RADIUS 3 Radius error of station 3 

LAT 3 Latitude error of station 3 

LONG 3 Longitude error of station 3 


MU-SUN 

MU-PLN 

A 

E 

I 

NODE 

OMEGA 

M 

RANGE 
R-RATE 
ST ANG 1 
ST ANG 2 
ST ANG 3 
APP DIAM 


Sun gravitational constant bias 

Target planet gravitational constant bias 

Target planet semi-major axis bias 

Target planet eccentricity bias 

Target planet inclination bias 

Target planet longitude of ascending node bias 

Target planet argument of periapsis bias 

Target planet mean anomaly bias 

Range bias of station 1 

Range-rate bias of station 1 

Star-planet angle 1 bias 

Star-planet angle 2 bias 

Star-planet angle 3 bias 

Apparent planet diameter bias 




(6) Definition of inertial frame. 

(7) Initial spacecraft position/velocity state vector in both 
input coordinates (ICOOR) and inertial coordinates (IBARY). 

(8) Nominal trajectory code (NTMC) . 

(9) List of celestial bodies assumed in the computation of the 
nominal trajectory (NB array). 

(10) Target planet (NTP). 

(11) Length units per A.U. (ALNGTH) and time units per day (TM) . 

(12) If the orbital elements of the assumed celestial bodies are to 
be computed at every time interval, a message to this effect 
will be printed. Otherwise, the orbital elements that will be 
used throughout the trajectory will be printed. 

(13) If output at the initial and final steps of the virtual mass 
trajectory is to be suppressed, a message to this effect will 
be printed. 

(14) If the virtual mass program will integrate only until the 
sphere of influence of the target planet is reached, a message 
to that effect will be printed. However, if the trajectory 
will continue until reaching a normal stopping condition, the 
appropriate message will be printed. 

(15) Trajectory accuracy figure (ACC). 

(16) Trajectory print intervals in days (DELTP) and increments (INPR). 

(17) Measurement schedule; measurement codes defined in section 
dealing with input description. 

(18) Schedule of eigenvector, prediction, and guidance events. 

(19) Sigma level of hyperellipsoid computed at an eigenvector event 
(IHYPl) . 

(20) Initial P, C ,C ,C ,P,C , U , and V covariance matrix 

xx xu xv s xu o o 

s s 

partitions; defined in section dealing with input description. 

(21) Definition of structure of augmented state transition, observation, 
and covariance matrices and their dimensions. 
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(22) State transition matrix code (ISTMC). If the state transition 
matrix is to be computed using numerical differencing, the 
position and velocity factors are also printed. 

(23) Dynamic noise constants used to compute the dynamic noise co- 
variance matrix if dynamic noise is non-zero. 

(24) Measurement noise for range, range-rate, star-planet angle, and 
apparent planet diameter measurements. 

(25) Tracking station locations. 

Measurement output 

Measurement information is printed every IPRINT measurements. At 
such a time the following information is printed: 

(1) Measurement number and corresponding trajectory time. 

(2) Type of measurement. 

(3) Trajectory time t 1 at most recent measurement or event (initial 
trajectory time). 

(4) Trajectory time t^ at present measurement (final trajectory time). 

(5) Initial and final spacecraft ecliptic position/velocity components 
and magnitudes relative to inertial space, the Earth, and the 
target planet. 

(6) If IPRT(4) = 1, all navigation parameters at the present measurement 
time will be printed. 

(7) State transition matrix partitions J", 9 , and 9 over the time 

XX XU 

interval frk-l* fc k] ’ relating deviations S in spacecraft position 

and velocity, solve-for parameters, and dynamic consider parameters 
at time t^ ^ to spacecraft position and velocity deviations at time 

t^. Note that transposed matrices are printed. 

(8) Diagonal of dynamic noise covariance matrix Qj represents unmodeled 
accelerations over the time interval I^k-l* * 



(9) Observation matrix partitions H, M, G, and L relating deviations 

in spacecraft position and velocity, solve-for parameters, dynamic 
consider parameters, and measurement consider parameters at time 
t^ to deviations in the observables at time t^. Note that trans- 
posed matrices are printed. 

(10) Measurement noise correlation matrix and standard deviations 
(covariance matrix R) . 

(11) Measurement residual correlation matrix and standard deviations 
(covariance matrix J). 

(12) Kalman gain matrix partitions. The K matrix is used in the fil- 
tering equations to compute the P, C , C , and C covariance 

XX XU XV 

s 

matrix partitions. The S matrix is used in the filtering equa- 
tions to compute the P , C , and C covariance matrix partitions. 
r s xu xv 

s s 

(13) Correlation matrix partitions and standard deviations at time t^, 
just before the measurement. The first group of correlation 
matrix partitions represents the correlation between spacecraft 
position and velocity and the variables listed in the left hand 
column: they are obtained by converting P, C , C , and C 

s 

into the corresponding correlation matrices and standard deviations. 
The second group represents the correlation between the solve-for 
parameters and the variables listed in the left hand column; they 
are obtained by converting P , C , and C into the correspond- 

3 S X u xv r 

s s 

ing correlation matrices and standard deviations. 

(14) Correlation matrix partitions and standard deviations at time t^, 
just after processing the measurement. See (13) above for 
definitions of the two groups of matrix partitions. 

c. Additional trajectory output 

If the spacecraft encounters the sphere of influence or closest approach 
during the course of the nominal trajectory, the information related to 
the encounter is printed. Also information normally printed during 
trajectory mode operation is printed during error analysis mode operation 
every DELTP days and INPR increments. This information includes space- 
craft inertial state, planet ephemeris data, spacecraft state relative 
to planets, virtual mass data, and navigation parameters, depending on 
the value of the IPRT vector. 
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d 0 Eigenvector event output 

At an eigenvector event the following information is printed: 

(1) Name of event and event time t 

v ' ev 

(2) Spacecraft position/velocity state vector relative to inertial 

space at event time t 
r ev 

(3) If IPRT(4) = 1, all navigation parameters at event time t will 
be printed. 

(4) State transition matrix partitions JT, 9 , and 9 over the 

XX XU 

s 

time interval It. , , t 1 , where t. . is the time of the most 
L k-1 evj k-1 

recent measurement or event. 

(5) Diagonal of dynamic noise covariance matrix Q; represents un- 
modeled accelerations over the time interval |t, , , t 1 . 

L k-1* ev J 

(6) Correlation matrix partitions and standard deviations at event 
time t gv propagated forward from time t^ See article (13) 

under measurement output for definitions of the two groups of 
matrix partitions. 

(7) Spacecraft position and velocity eigenvalues, square roots of 
eigenvalues, and eigenvectors at event time as specified by 
IEIG code. 

(8) Hyperellipsoids for both position and velocity covariance matrix 
partitions at event time for the sigma level specified by IHYPl 
code. 

e. Prediction event output 

At a prediction event the following information is printed: 

(1) Name of event, event time t , and time t to which prediction 

is Demg made. r 

(2) Articles (2) through (6) under eigenvector event output. 

(3) State transition matrix partitions § , 9 , and 9 over the time 

XX XU 

interval |t , t . 

*- ev pJ 
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(A) Diagonal of dynamic noise covariance matrix Q; represents 

unmodeled accelerations over the time interval [t , t ]. 

ev p 

(5) Correlation matrix partitions and standard deviations at 

time t based on prediction from time t . See article 
p ev 

(13) under measurement output for definitions of the two 

groups of matrix partitions. 

(6) Spacecraft position and velocity eigenvalues, square root 

of eigenvalues, and eigenvectors at time t as specified 
by IEIG code. P 


(7) Hyperellipsoids for both position and velocity covariance 
matrix partitions at time t for the sigma level specified 
by IHYP1 code. P 


(8) If time t^ occurs within the target planet sphere of in- 
fluence, the Cartesian position/velocity correlation 
matrix and standard deviations are transformed to B-plane 
coordinates B*T, B*R, time-of-f light, S*R, S*T, and C^- 

The transformation matrix relating these coordinates to 
Cartesian position/velocity coordinates is printed, fol- 
lowed by the B-plane correlation matrix and standard de- 
viations. The semimajor axis, semiminor axis, and orien- 
tation of the B-plane 1-a uncertainty ellipse are also 
printed. 


Output preceding all types of guidance events 

At a guidance event the following information is printed: 

(1) Articles (1) through (8) under eigenvector event output. 

(2) State transition matrix partitions over the time inter- 
val [t , t ], where t is the time of the previous gui- 

g ev J g r 6 

dance event (t = t if no guidance event has occurred 
previously) . ® 3 

(3) Diagonal of dynamic noise covariance matrix Q; represents 

unmodeled accelerations over the time interval [ t , t ] . 

g ev 

(A) Control correlation matrix partitions and standard devia- 
tions at time t , just before the guidance correction 

is applied. See article (13) under measurement output for 



definitions of the two groups of matrix partitions. 
Eigenvalues, eigenvectors, and hyperellipsoids are also 
printed . 

(5) Description of guidance event: 

(a) Guidance policy 

(b) Linear or nonlinear guidance 

(c) Status of planetary quarantine constraints 

(d) Thrust model 

(e) Guidance event treatment 

g. Linear midcourse guidance event output 

Three midcourse guidance policies are available: fixed-time- 

of-arrival (FTA) , two-variable B-plane (2VBP), and three- 
variable B-plane (3VBP) . 

(1) Time, position, and velocity when spacecraft encounters 
closest approach at target planet if FTA. Time, posi- 
tion, and velocity when spacecraft pierces target planet 
sphere of influence, together with B, B*T, and B*R, if 
2VBP or 3VBP. 

(2) Matrix M relating position/velocity deviations at t to 
deviations in B*T and B*R. 

(3) State transition matrix partitions over [t , t ] if 

FTA; over [t t CT ] if 2VBP . ev CA 

ev SI 

(4) Variation matrix n (or partitions) relating position/ 

velocity deviations at time t to target condition 
, . ev 

deviations . 

(5) Target condition correlation matrix and standard devia- 
tions (covariance matrix W ) immediately prior to gui- 
dance correction, together with eigenvalues, eigenvectors, 
and hyperellipsoid. 

(6) Guidance matrix T used to compute the velocity correction 
required to null out target condition deviations. 
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(7) Velocity correction correlation matrix and standard de- 
viations (covariance matrix S) , together with eigenvalues 
and eigenvectors. The hyperellipsoid is also printed if 
the guidance policy is not 2VBP. 

(8) Expected value of the effective velocity correction. 

(9) Execution error correlation matrix and standard devia- 
tions (covariance matrix Q) . 
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(10) Control (and knowledge) correlation matrix partitions and 
standard deviations just after the guidance correction at 
time t , together with eigenvalues, eigenvectors, and 

hyperellipsoids . 

(11) Target condition correlation matrix and standard deviations 
(covariance matrix W ) just after guidance correction is 

applied, together with eigenvalues, eigenvectors, and hyperellipsoid. 

(12) Targeted nominal trajectory just after guidance correction. 

h. Orbital insertion guidance event output 

(1) Orbital insertion decision output: see NOMNAL output description. 

(2) Orbital insertion execution output: 

(a) Orbital insertion AV in both ecliptic and equatorial 
coordinates . 

(b) Execution error correlation matrix and standard deviations. 

(c) Spacecraft position/velocity state relative to target 
planet immediately after orbital insertion in both ecliptic 
and equatorial coordinates. 

(d) Spacecraft orbital elements relative to the planeto-centric 
equatorial coordinate system immediately after orbital 
insertion. 

i. Nonlinear guidance and retargeting event output is summarized in the 
NOMNAL output description. 

j. Pulsing thrust model output. 

(1) Pulsing arc output: see NOMNAL output description. 

(2) Pulsing arc covariance output: 

(a) Nominal state transition matrix over the time interval separate 
ing pulses. 

(b) Error model variances. 

(c) Nominal pulse components and magnitude. 

(d) Nominal execution error covariance matrix. 
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(e) Final pulse components and magnitude. 

(f) Final execution error covariance matrix. 

(g) Accumulated execution error covariance matrix. 
k„ Biased aimpoint guidance event output. 

(1) Target planet capture radius. 

(2) Matrix ^ used to compute the velocity correction required to 
achieve the biased aimpoint. 

(3) Previously imparted aimpoint bias, 5 m, if any. 

(4) Velocity correction Av^ required to remove the previously 
imparted aimpoint bias. 

(5) Execution error covariance matrix associated with the velocity 
correction required to remove bias and null out nominal target 
errors . 

(6) Covariance matrix the projection of the target condition 
covariance matrix W into the impact plane. 

(7) Probability of impact if the spacecraft is targeted to the 
nominal aimpoint. 

Articles (8) through (13) appear for each iteration which occurs 

within subroutine BIAIM. 

(8) Equation defining the probability constraint ellipse. 

4 . * 

(9) Biased aimpoint m and aimpoint bias &M. 

(10) Velocity correction AV, required to achieve the biased 

» ■ , D 1.3*S 

aimpoint . 

(11) Execution error covariance matrix associated with the velocity 
correction required to achieve the biased aimpoint and null 
out target errors. 

(12) Covariance matrix A . 

(13) Probability of impact if the spacecraft is targeted to the 
biased aimpoint. 
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1. Probe release event output 


The probe release event output for both main probe and mini- 
probe release events is essentially the same. Differences are 
noted in the following probe release event output description: 


( 1 ) 


Type of probe release event and trajectory time t at 
release. J 


( 2 ) 

(3) 


Heliocentric ecliptic probe state at 


If IPRT(4) = 1, all navigation parameters at t will be 
printed. ^ 


(4) State transition matrix partitions over the interval 

[ t, , t.], where t, is the time of the previous measure- 

ment or event. 


(5) Probe knowledge correlation matrix partitions and stand- 
ard deviations at t . . 

3 

(6) Eigenvalues, eigenvectors, and hyperellipsoids of the 
position and velocity partitions of the probe position/ 
velocity knowledge covariance matrix. 

(7) If a miniprobe is being treated and if miniprobe target- 
ing has been performed internally (i.e., user-specified 
target controls are not available) , the following infor- 
mation is printed, 


a) Bus state at probe sphere relative to target planet 
obtained by propagating bus state at t^ forward to 

entry. B-plane coordinates of bus at entry. 


b) Miniprobe target controls, including release roll 
angle of miniprobe No. 1 (radians) , tangential ve- 
locity magnitude of miniprobe at release ( decameters /s) , 
declination of spin axis (radians) , and right ascen- 
sion of spin axis (radians) , 


c) Value of KKWIT. If KKWIT = 0, miniprobe targeting 
was successful. If KKWIT «• 1, targeting failed. 


(8) If a miniprobe is being treated, the execution error 
covariance matrix for the ith miniprobe is printed. 
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( 9 ) 


Probe plane tocentric ecliptic state relative to target 
planet when probe sphere is pierced at entry trajectory 
time t . Probe B-plane coordinates at entry. 

(10) Probe planetocentric ecliptic state (same as item 9)). 

(11) Julian date (epoch 1900) and trajectory time at entry. 
Probe sphere radians (AUs) . 


(12) State transition matrix partitions over [t., t ]. 

J E 

(13) Probe control correlation matrix partitions at t £ . 

(14) Julian date (epoch 1900) and calendar date at t . 

E 

(15) Probe state relative to target planet at t^ in planeto- 
centric ecliptic coordinates. 

(16) Probe state relative to target planet at t^ in subsolar 
orbital-plane coordinates; 


(17) 


Probe communication angle at 


V 


(18) Transformation matrix relating Cartesian position and 
velocity coordinates to (LTR) entry parameters h, v, 
y, and <f> g at entry. 


(19) Entry parameters h, v, y, <j> g , ft g , and i at t £ . 

(20) Entry parameter control covariance matrix at t £ for 
entry parameters h, v, y, and <|> g . 


(21) Measurement output for every IPRINT probe measurement over 

the interval [t. t ]. Format is identical to the standard 
J E 

ERRAN measurement output format. 


(22) Probe planetocentric ecliptic state relative to target 

planet at t„. Probe B-plane coordinates at entry. 

E 

(23) Items 1) through 20) with the words "control correlation" 
replaced with "knowledge coorelation. " 


(24) If probe release event is a miniprobe release event and 

if not all miniprobes have been treated (three miniprobes), 
return to item 8) and continue. 
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m. Summary output 


At a successful conclusion of an error analysis run, the fol- 
lowing error analysis summary is printed: 

(1) Method used to compute nominal trajectory. 

(2) Trajectory accuracy figure (ACC) and true anomaly incre- 
ment. 

(3) Length units per AU; time units per day. 

(4) Method used to compute orbital elements of planets. 

(5) Initial and final trajectory times, calendar dates, and 
Julian dates. 

(6) Inertial coordinates of spacecraft position and velocity 
at initial and final times. 

(7) Spacecraft position and velocity coordinates relative to 
both Earth and target planet at final time. 

(8) Time of closest approach and position and velocity rela- 
tive to target planet at closest approach. 

(9) If the spacecraft did not reach the target planet sphere 
of influence, a message to that effect is printed. Other- 
wise, the time at which the sphere of influence is pierced, 
together with position and velocity relative to the target 
planet at that time B, B*T, and B*R, are printed. 

(10) Method used to compute state transition matrix, together 
with associated information. 

(11) Number of measurements taken. 

(12) Number of events having occurred and the number of each 
type of event. 

(13) Variances used for the resolution, proportionality, and 
pointing angle errors in guidance events. 

(14) Locations of tracking stations. 

(15) Dynamic and measurement noise constants. 


54-1 



(16) Direction cosines of the three reference stars. 

(17) Lists of solve-for, dynamic consider, and measurement 
consider parameters. 

(18) Correlation matrix partitions and standard deviations 
at final time. 
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4.3 SIMUL Output Description 

The printed output of the simulation mode is described in this section 
according to the following groups: input data, measurement output, addi- 

tional trajectory output, eigenvector event output, prediction event output, 
guidance event output, quasi-linear filtering event output, and summary 
output. 

a. Input data 

The initial output consists of the following input data: 

(1) Namelist SIMUL. 

(2) All input data printed out in the error analysis mode. 

(3) Actual trajectory information: 

(a) List of celestial bodies assumed in the computation of 
the actual trajectory (NBl array). 

(b) Trajectory accuracy figure (ACGl). 

(c) Actual measurement biases (BIA array). 

(d) Dynamic biases. 

(e) Description of unmodeled acceleration characteristics. 

(f) Station location biases (SLB array). 

(g) Actual spacecraft position and velocity deviations from the 
most recent nominal trajectory at the initial time (ADEVX 
array) . 

(h) Actual measurement noise variances (AVARM array). 

b. Measurement output 

Measurement output in the simulation mode repeats the measurement output 
in the error analysis mode except for the following differences: 

(1) Spacecraft states are given for three trajectories: targeted 

nominal, most recent nominal, and actual. 

(2) Navigation parameters are based on the actual trajectory. 
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(3) State transition, observation, and covariance matrices are based 
on the most recent nominal trajectory. 

In addition, the following output appears: 

(4) Actual dynamic noise representing effect of the actual un- 
modeled acceleration. 

(5) Actual measurement noise correlation matrix and standard deviations 
(covariance matrix R) . 

(6) Actual measurement noise V . 

(7) Estimated and actual measurements. Measurement residuals c . 

(8) Estimated and actual spacecraft position/velocity deviations 
from the most recent and the targeted nominal trajectories. 

Actual orbit estimation error. 

(9) Estimated and actual solve-for parameter deviations. Actual 
estimation error. 

c. Additional trajectory output appearing in the simulation mode is identical 

to that appearing in the error analysis mode. 

d. Eigenvector event output. 

Eigenvector event output appearing in the simulation mode is identical 

to that appearing in the error analysis mode with the following additions: 

(1) Spacecraft position/velocity states for three trajectories: 
targeted nominal, most recent nominal, and actual. 

(2) Actual dynamic noise representing effect of the actual unmodeled 
acceleration. 

(3) Estimated and actual spacecraft position/ velocity deviations from 
the most recent nominal trajectory. 

(4) Estimated and actual solve-for parameter deviations. 

e. Prediction event output 

Prediction event output appearing in the simulation mode is identical 

to that appearing in the error analysis mode with the following additions: 
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(1) Spacecraft position/velocity states for three trajectories: 
targeted nominal, most recent nominal, and actual. 

(2) Actual dynamic noise representing effect of the actual dynamic 
noise at the time of the event t 

ev 

(3) Estimated and actual spacecraft position/velocity deviations 
from the most recent nominal trajectory at t . 

(4) Estimated and actual solve-for parameter deviations at t . 

(5) Most recent nominal trajectory state and the estimated deviations 
from this state at time t to which the prediction is being made. 

f. Output preceding all types of guidance events. 

(1) All eigenvector event output. 

(2) State transition matrix partitions over the time interval 

where t is the time of the previous guidance event (t = 

§ g 

no guidance event has previously occurred). 

(3) Diagonal of dynamic noise covariance matrix Q. 

(4) Control correlation matrix partitions and standard deviations at 
time t , just before the guidance correction is applied. 

(5) Eigenvalues, eigenvectors, and hyperellipsoids associated with 
the position and velocity partitions of the previous control 
covariance. 

(6) Description of guidance event: 

(a) Guidance policy 

(b) Linear or nonlinear guidance 

(c) Status of planetary quarantine constraints. 

(d) Thrust model 

(e) Guidance event treatment 

g. Linear midcourse guidance event output 

Three midcourse guidance policies are available: fixed-time-of -arrival 

(FTA), two-variable B-plane (2VBP), and three-variable B-plane (3VBP). 


[V C ev] ’ 


g 

t if 
o 
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(1) Closest approach conditions on targeted nominal trajectory if 
FTA. Sphere of influence conditions on targeted nominal if 
2VBP or 3VBP. 

(2) Matrix M relating position/velocity deviations at sphere of 
influence to B»T and B*R deviations. 

(3) Closest approach conditions on most recent nominal trajectory if 
FTA. Sphere of influence conditions on most recent nominal 
trajectory if 2VBP or 3VBP. 

(4) State transition matrix partitions over j^t^, t^J if FTA; over 
£t ev> t gjJ if 2VBP. Based on most recent nominal trajectory. 

(5) Variation matrix i? (or partitions) relating position/velocity 

deviations at time t to target condition deviations. 

ev 

(6) Target condition correlation matrix and standard deviations 
(covariance matrix W ) immediately prior to guidance correction, 
together with eigenvalues, eigenvectors, and hyperellipsoid, 

(7) Guidance matrix F used to compute the velocity correction required 
to null out target condition deviations. 

(8) Velocity correction correlation matrix and standard deviations 
(covariance matrix S), together with eigenvalues and eigenvectors. 
The hyperellipsoid is also printed if the guidance policy is not 
2VBP. 

(9) Estimated and actual spacecraft position/velocity deviations from 
targeted nominal immediately prior to the guidance correction. 

(10) Commanded and perfect velocity corrections Av and AV . 

(11) Magnitude of commanded velocity correction. 

(12) Error in velocity correction, AV g , due to navigation uncertainty. 

(13) Execution error correlation matrix and standard deviations 
(covariance matrix Q) . 

(14) Control (and knowledge) correlation matrix partitions and standard 
deviations just after the guidance correction at time t , together 

with eigenvalues, eigenvectors, and hyperellipsoids. 
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(15) Actual velocity correction execution error 6AV. 

(16) Actual velocity correction AV. 

(17) Target condition correlation matrix and standard devia- 

tions (covariance matrix W ) just after the guidance cor- 
rection is applied, together with eigenvalues, eigenvectors, 
and hyperellipsoid. 

(18) Actual target errors e and e due to navigation un- 

° nav ex 

certainty and execution error. Total target error. 

(19) Most recent and targeted nominal trajectories immediately 
following guidance correction. 

(20) Actual and estimated spacecraft position/velocity devia- 
tions from most recent nominal trajectory immediately 
following guidance correction. 

h. Orbital insertion guidance event output 

(1) Orbital insertion decision output: see NOMNAL output 

description. 

(2) Orbital insertion execution output: 

(a) Actual orbital insertion AV in both ecliptic and 
equatorial coordinates. 

(b) Execution error correlation matrix and standard 
deviations . 

(c) Actual spacecraft position/velocity state relative 
to target planet immediately after orbital insertion 
in both ecliptic and equatorial coordinate. 

(d) Actual spacecraft orbital elements relative to the 
planetocentric equatorial coordinate system immedi- 
ately after orbital insertion. 

i. Nonlinear guidance and retargeting event output is summarized 
in the NOMNAL output description. 

j. Pulsing thrust model output is identical to that appearing in 
the error analysis mode except that pulsing arc information 
is printed for both the estimated and actual trajectories. 
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k. Biased aimpoint guidance event output is identical to that 

appearing in the error analysis mode. 

l . Quasi-linear filtering event output 

(1) Spacecraft position/velocity states for three trajectories 
targeted nominal, most recent nominal, and actual. 

(2) Navigation parameters based on actual trajectory. 

(3) State transition matrix partitions over the interval 

[t, , t 1, where t. is the time of the last event or 
1 k’ ev k 

measurement. 

(4) Correlation matrix partitions and standard deviations 
at the time of the event, together with eigenvalues, 
eigenvectors, and hyperellipsoids. 

(5) Actual dynamic noise representing effect of the actual 
unmodeled acceleration. 

(6) Estimated and actual spacecraft position/velocity devia- 
tions from the most recent nominal trajectory just prior 
to the event. 

(7) Estimated and actual solve-for parameter deivations just 
prior to the event . 

(8) Most recent nominal trajectory just after the event. 

(9) Estimated and actual spacecraft position/velocity devia- 
tions from the most recent nominal trajectory just after 
the event. 

(10) Estimated and actual solve-for parameter deviations just 
after the event. ' 


m. Probe release event output 


Probe release event output in SIMUL is identical to probe re- 
lease event output in ERRAN, except for the following items: 


(1) Targeted nominal, most recent nominal, and actual probe 
states are printed at release time t, instead of just 


the targeted nominal. 


j 
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(2) If a miniprobe release event is being executed, the 
actual miniprobe execution error is printed in addition 
to the execution error covariance matrix. 

(3) Format for probe measurement output is identical to the 
standard SIMUL measurement output format. 

(4) Quasi-linear filtering event output for a probe is iden- 
tical to standard SIMUL quasi-linear filtering event out- 
put. 

Summary output 

(1) Accuracies used in nominal and actual trajectory compu- 
tation. 

(2) Bodies treated in nominal and actual trajectory computa- 
tion. 

(3) Gravitational constant biases used in actual trajectory. 

(4) Ephemeris biases used in actual trajectory. 

(5) Initial trajectory time. 

(6) Final trajectory time. 

(7) Position and velocity of vehicle relative to sun, earth, 
and target planet at initial time. 

(8) Position and velocity of vehicle relative to sun, earth, 
and target planet on targeted nominal, most recent nomi- 
nal, and actual trajectory at final time. 

(9) Time at closest approach plus position and velocity of 
vehicle relative to target planet on all three trajector- 
ies. 

(10) The time at which the vehicle enters the sphere of influ- 
ence of the target planet in addition to the position 
and velocity of the vehicle relative to the target planet 
and B, B*T, and B*Ron all three trajectories. 

(11) Method by which the state transition matrix is computed 
in addition to its limitations. 



(12) Number of measurements taken. 

(13) Number of events plus the number of each type of event. 

(14) Variances of errors used in guidance events. 

(15) Actual errors used in guidance events. 

(16) Station location constants. 

(17) Dynamic noise constants. 

(18) Actual unmodeled acceleration. 

(19) Assumed measurement noise constants. 

(20) Actual measurement noise constants. 

(21) Direction cosines for three star planet angles. 

(22) Initial state vector for both nominal and actual 
trajectories . 

(23) Final state vector for all three trajectories. 

(24) Estimated and actual deviations from most recent nominal 
at final time. 

(25) Estimated and actual deviations from targeted nominal 
at final time. 

(26) Actual orbit determination error at final time. 

(27) Initial correlation matrix partitions and standard 
deviations . 

(28) Final correlation matrix partitions and standard devia- 
tions . 
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4.4 GENC0V Output Description 

Although the generalized covariance program GENC0V is actually 
a part of the error analysis program ERRAN, for the purpose of 
clarity and convenience to the user, the GENC0V output will be de- 
scribed separately in this section according to the following groups 
Input data, measurement output, and guidance event output. Output 
for eigenvector and prediction events will not be described be- 
cause of its similarity to standard ERRAN eigenvector and predic- 
tion event output. The only difference consists in the fact that 
both assumed and actual statistics are printed for GENC0V eigen- 
vector and prediction events. 

a. Input data 

The initial output consists of the following input data: 

(1) Namelist ERRAN. 

(2) Calendar date and Julian date at launch or initial time. 

(3) Final calendar date and Julian date. 

(4) Initial trajectory time in days (TRTM1) . 

(5) Lists of solve-for, dynamic-cons ider, measurement- 
consider, and ignore parameters. Definitions of param- 
eter names can be found in subsection a of the ERRAN 
output description. 

(6) Definition of inertial frame. 

(7) Initial spacecraft position/velocity state vector in both 
input coordinates (ICOOR) and inertial coordinates (IB ARY) . 

(8) Nominal trajectory code (NTMC) . 

(9) List of celestial bodies assumed in the computation of 
the nominal trajectory (NB array) . 

(10) Target planet (NTP) . 

(11) Length units per AU (ALNGTH) and time units per day (TM) . 



(12) If the orbital elements of the assumed celestial bodies 
are to be computed at every time interval, a message to 
this effect will be printed. Otherwise, the orbital 
elements that will be used throughout the trajectory will 
be printed. 

(13) If output at the initial and final steps of the virtual 
mass trajectory is to be suppressed, a message to this 
effect will be printed. 

(14) If the virtual mass program will integrate only until 
the sphere of influence of the target planet is reached, 
a message to that effect will be printed. However, if 
the trajectory will continue until reaching a normal 
stopping condition, the appropriate message will be 
printed. 

(15) Trajectory accuracy figure (ACC). 

(16) Trajectory print intervals in days (DELTP) and incre- 
ments (INPR) . 

(17) Measurement schedule; measurement codes defined in sec- 
tion dealing with input description. 

(18) Schedule of eigenvector, prediction, and guidance events. 

(19) Sigma level of hyperellipsoid computed at an eigenvector 
vent (IHYP1) . 

(20) Initial assumed covariance matrix partitions; defined 
in input description section. 

(21) Definition of structure of augmented state transition, 
observation, and assumed covariance matrices and their 
dimensions . 

(22) State transition matrix code (ISTMC) . If the state 
transition matrix is to be computed using numerical dif- 
ferencing, the position and velocity factors are also 
printed. 

(23) Dynamic noise constants used to compute the assumed dy- 
namic noise covariance matrix if dynamic noise is non- 
zero. 
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(24) Assumed measurement noise variances for range, rante-rate, 
star-planet angle, and apparent planet diameter measure- 
ments . 

(25) Nominal tracking station locations. 

(26) Namelist GENRAL. 

(27) Initial position/velocity, solve-for, dynamic-consider, 
measurement-consider , and ignore parameter deviation 
means . 

(28) Initial actual covariance matrix partitions; defined in 
input description section. 

(29) Definition of structure of augmented actual covariance 
matrix and dimensions of each partition. 

(30) Dynamic noise constants used to compute the actual dy- 
namic noise covariance matrix if dynamic noise is non- 
zero. 

(31) Actual measurement noise variances for range, range-rate, 
star-planet angle, and apparent planet diameter measure- 
ments . 

b. Measurement output 

Measurement information is printed every IPRINT measurements. 

At such a time the following information is printed: 

(1) Measurement number and corresponding trajectory time. 

(2) Measurement type. 

(3) Trajectory time at most recent measurement or event. 

(4) Trajectory time t^ at present measurement. 

(5) Initial and final spacecraft position/velocity components 
and magnitudes relative to inertial space, the earth, 
and the target planet. 

(6) Elevation and azimuth of spacecraft relative to the track- 
ing station if a range or range-rate measurement is being 
processed. 
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(7) If IPRT(4) = 1, all navigation parameters will be printed. 

(8) State transition matrix partitions 8 ,0 , and 0 

V ' XX XU xw 

s 

t^] , relating spacecraft 

position/velocity, solve-for, dynamic-consider, and ignore 
parameters at t^ respectively. Note that transposed 

matrices are printed. 

(9) Diagonal of assumed dynamic noise covariance matrix Q. 

(10) Observation matrix partitions H, M, G, L, and N, relating 
deviations in spacecraft position and velocity and solve- 
for, dynamic-consider paramters, measurement-consider and 
ignore parameters, at time t^ to deviations in the observ- 
ables at time t^. Note that transposed matrices are 
printed. 

(11) Assumed measurement noise correlation matrix and standard 
deviations (covariance matrix R) . 

(12) Assumed measurement residual correlation matrix and stand- 
ard deviations (covariance matrix J) . 

(13) Kalman gain matrix partitions. 

(14) Assumed correlation matrix partitions and standard de- 
viations at time t, just before the measurement. The 

k 

first group of correlation matrix partitions represents 
the correlation between spacecraft position and velocity 
and the variables listed in the left-hand column. The 
second group represents the correlation between the solve- 
for parameters and the variables listed in the left-hand 
columns . 

(15) Assumed correlation matrix partitions and standard devia- 
tions at time t^ just after processing the measurement 

(see 14) for definitions of the two groups of matrix 
partition^) . 

(16) Diagonal of actual dynamic noise covariance matrix Q' . 

(17) Actual measurement noise correlation matrix and standard 
deviations (covariance matrix R’) . 


over the time interval [ t^ 
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(18) Actual measurement residual mean E[e']. 

(19) Actual measurement residual correlation matrix and 
standard deviations (2nd moment matrix J'). 

(20) Actual estimate error means at time t^ just before the 

measurement for both position/velocity and solve-for 
parameter states. 

(21) Actual correlation matrix partitions and standard de- 
viations at time t, just before the measurement. 

k 

The first group of correlation matrix partitions repre- 
sents the correlation between spacecraft position and 
velocity and the variables listed in the left-hand column. 
The second group represents the correlation between the 
solve-for parameters and the variables listed in the 
left-hand column. 

(22) Actual estimation error means at time t^ just after 

processing the measurement, for both position/velocity 
and solve-for parameter states. 

(23) Actual correlation matrix partitions and standard devia- 
tions at time t^ just after processing the measurement. 

See 21) for definitions of the two groups of matrix 
partitions . 

c. Guidance event output 

Generalized covariance analysis information relating to the 
execution of the guidance event is printed immediately after 
the standard ERRAN guidance event information has been printed. 
This standard guidance event output, which is described in 
the ERRAN output description, comprises the assumed guidance 
data in contrast to the actual guidance data generated by the 
generalized covariance analysis. The generalized covariance 
analysis guidance event output for a midcourse guidance policy 
follows. The output for other guidance policies is a subset 
of this output. 

(1) Actual position/velocity and solve-for parameter devia- 
tion means just before the guidance correction. 

(2) Actual control correlation matrix partitions and stand- 
ard deviations just before the guidance correction. 
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( 3 ) 


Eigenvalues, eigenvectors, and hyperellipsoids of the 
position and velocity partitions of the actual position/ 
velocity control covariance matrix. 

(4) Actual target state deviation mean, E[6 t* ], just before 
the guidance correction. 

(5) Actual target condition correlation matrix and standard 
deviations just before the guidance correction (2nd mo- 
ment matrix W* ) . 

(6) Eigenvalues, eigenvectors, and hyperellipsoid of actual 
target condition covariance matrix. 

(7) Actual velocity correlation 2nd moment matrix S', together 
with eigenvalues and eigenvectors. 

(8) Actual velocity correction correlation matrix and stand- 
ard deviations (2nd moment matrix S'). 

(9) Mean of actual commanded velocity correction, E[AV']. 

(10) Mean of magnitude of actual commanded velocity correction, 
E [ | AV T | ] . 

(11) Actual statistical, or effective, velocity correction, 
"E[AV']." 

(12) Actual execution error mean, E[6AV']* 

(13) Actual execution error correlation matrix and standard 
deviation (2nd moment matrix Q' ) . 

(14) Actual position/velocity deviation means just after the 
guidance correction. 

(15) Actual position/velocity estimation error means just af- 
ter the guidance correction. 

(16) Actual control (and knowledge) correlation matrix parti- 
tions and standard deviations just after the guidance 
correction. 

(17) Eigenvalues, eigenvectors, and hyperellipsoids of the 
position and velocity partitions of the actual position/ 
velocity control (and knowledge) covariance matrix. 
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(18) Actual target state deviation mean, E[6 t' + ], just after 
the guidance correction. 

(19) Actual target condition correlation matrix and standard 
deviations just after the guidance correction (2nd 

moment matrix W ,+ ) . 

(20) Eigenvalues, eigenvectors, and hyperellipsoid of actual 
target condition covariance matrix. 



5 . SAMPLE CASES 

5.1 NOMNAL Sample Cases 

Three typical trajectories generated by KOMNAL will be described in this 
section to illustrate the operation and versatility of the nominal trajectory 
generator NGMNAL. The three cases to be discussed are: 

Case N-l. Broken Plane Viking Mars "75 Mission 
Case N-2. Planetary Explorer Venus '78 Mission 
Case N-3. Lunar Viking '76 Mission 

5.1.1 Broken Plane Viking Mars '75 Mission 
a. Sample Data 


NB0D=3,NB=1, 4,5, NLP=4 , NT P=5 , ACKT=2. 5 E- 5 , NCPR= 1 00 0 0 , TMPR = 10 . , 
KALI=1975, 8, 30, 0,0,31=0. , 

IZER0=2, ZD AT (4) =-4. 201 823E+6 , 1 . 9823 1564E + 8, - 7. 118753E+6 , 
KTrP=l,2,3,-l, KMXQ=3 ,3,4, KT IM= 1 ,4,2,3, PERV=1 . E-5 , 5 . E-5 , 
IBAPS = 2,2, LVLS=2,2, AC ( i , 1 ) =5 . E-4, 2 . 5E-5 , AC ( 1 , 2> =5. E- 4 ,2. 5E-5 , 
KAIG<1,2)=1976„1,5,9,59, GS<2)=4.901, T IMG ( 3 ) =2 . 0 , . 5 , IZER(2)=3, 
KT ARC 1,1) =10, 11, 12, I ARC 1,1) =-4 2018 23. , 198231 564 ., -71 187 53 . , 
KTARC1,2)=8,7,3, T AR ( 1 , 2 ) =40 . 92 , 5000., T0L(1,2>=1., 10., .001, 

KT AR (1, 3 ) = 2, T ARC 1,3) =20 42 8. ,.76,77. ,40. ,50., TOL Cl, 1) =3*10 0 . , 
KALT=1976,1,5,9,59, TS=4.901, KAL T < 1 , 2 ) = 1976 , 7 , 1 9 , 0 , 0 , TSC2)=0., 


The exact data as read in for the Broken Plane Viking case are given above. 

A detailed explanation of this data follows. 

The first line defines the nominal trajectory propagation between guidance 
events. The sun, Earth, and Mars are the gravitational bodies with the 
Earth as the launch planet and Mars as the target planet. The integration 
level of 2.5E-5 is a moderate accuracy level. The trajectory will be recorded 
at intervals of ten days with no printouts occurring on integration increment 
counts . 

The initial date is given on the next line. It is specified only to a 
calendar day; the hours, minutes, and seconds at injection will be computed 
in the zero iterate computation using the internally set launch profile with 
a Cape Kennedy launch. 

The third line defines the zero iterate computation. The option specified 
(IZER0=2) specifies the launch planet to prescribed point option. The helio- 
centric ecliptic coordinates of that point are provided in ZDAT. 

The next lines of input define the guidance events. The data defining each 
event will be discussed in the order of the indices of the events. 
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The first event will be a targeting event (KTYP=1) to occur at 0 days 
(TIMG=0,) after the initial time (KTIM=1) . The correction is to be computed 
and executed (KMXQ=3) using the impulsive model (MDL undefined, hence set 
to 1).' The target values are x^, y^, z (KTAR=10 , 11 , 12) with values 


identical to ZDAT and tolerances of 100 km (T0L=3*100) . The final time is 
read in the KALT, TS arrays (and used incidentally in the zero iterate 
computation). The Newton-Raphson scheme is to be used (METH not set, hence 
equal to 0) with a perturbation size of 10 - ^ km/sec (PERV) at two levels 
(LVLS) defined by AC to be 5x10"^, 2.5x10“^. Bad-step checks will be made 
at the high level only (IBADS=2) . 


The second event will be a retargeting event (KTYP=2) to occur on the 
calendar date (KTIM=2) specified by KALG, SG which is the same time as the 
target time of the first event. The correction is to be computed and 
executed (KMXQ=3) using the impulsive model (MDL=1 since undefined). The 
target parameters are i, r ^ , t^ g (KTAR=8,7,3) with target values of 40,92°, 

5000 km (TAR) and 7/19/1976 (KALT) respectively and tolerances of 1°, 10 km, 
.001 days (TOL) . The scheme to be used is identical to that of the first 
event with the exception of the velocity perturbation which is now set to 
5x10 km/sec (PERV). 


The third event is an orbit insertion event (KTYP=3) occurring .5 days (TIMG) 
after intersection of the Martian SOI (KTIM=2) . The insertion is to be the 
nonplanar option (KTAR=2) with target conic elements of 

a = 20,428. i = 40° 

e = .76 Q = 50° 

CJ= 77° 

The fourth event is a termination event (KTYP= -1) occurring at 0 days 
(TIMG=0 since undefined) after closest approach to Mars (KTIM=3) . Thus the 
nominal trajectory will be integrated and recorded to the Martian closest 
approach . 

b. Sample Output 


Selected pages of the actual output from this run are supplied in the 
Appendix to this volume. 


c. Discussion 

The Broken Plane Viking Mars '75 mission data'*' may be summarized as follows: 
Launch date: 8/30/1975 

Broken plane date: 1/5/1976 9-59 m -4 . 901 sec 

Encounter date: 7/19/1976 


The authors are indebted to R. T. Gamber of the Martin-Marietta Corporation 
who generated this minimum-delta V data based on mass less- planets trajectories. 
The sum of the broken-plane velocity and orbit insertion velocity is minimized. 
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Launch planet: Earth 

Broken plane point (heliocentric ecliptic) ? 
Target planet: Mars 


1.38780(8), -5.96016(7), 
1.49455(2) km 


The massless planet trajectory generated for the first leg of the mission 
led to the following heliocentric conic: 


SMA: 

198.29 M km 

TAL: 

348 

.90° 

ECC: 

.24094 

TAP: 

104 

08° 

INC: 

2.2722° 

TOF: 

128 

42 days 


The injection conditions computed using the massless planet heliocentric 
trajectory and the internally stored launch profile generated the following 
near-earth conic: 

C3 : 14.242 SMA: 

RAD: 6567.6 ECC : 

VEL: 11.646 , 

INJ TIME: 18 hr -26 m -51 SeC 

When the injection state consistent with this conic is integrated by NGMNAL 
to the broken plane time, the error in position is 2.18xlCr km. Three 
iterations are made at both the first and second accuracy levels to obtain 
a trajectory that has an error of 25 km. The A v required and the elements 
of the corrected near earth conic are 


Av = 

13.6 m/sec 

SMA =' 

-28085.0 km 

ECC = 

1.23356 


The corrected nominal trajectory is now integrated to the time of the broken 
plane point. The second guidance event occurs at this time. First, a massless- 
planet trajectory is determined between the current position and the location 
of Mars at the encounter time. The elements of this conic are 


SMA: 

200.21 M km 

TAL: 

101.83° 

ECC: 

.24069 

TAP: 

191.76° 

INC: 

2.4785 

TOF: 

195.58 days 


-27988.6 

1.23438 


The velocity thus generated at the broken plane point by the virtual mass 
trajectory was altered by the massless planet correction (Av = 222 m/sec) 
before integrating to the target planet for the second guidance event. The 
target errors on the first propagation are A i = 40°, Ar = 3906 km, 

At = .003 day. Three iterations at the first accuracy level and two 

iterations at the second are required to reduce the errors to Ai = .000, 

far = 1.1 km, At = .000. The additional correction in velocity was 
CA CS 

A v = 4.2 m/ sec . 
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The corrected nominal trajectory is now integrated to the time of the 
orbit insertion decision (.71 days before CA) . The elements of the approach 
hyperbola, the target orbit, the modified orbit, and the orbit actually 
achieved upon the later execution are compared below: 


Orbit 

a 

e 

U) 

i 

Q 

Approach hyperbola 

-8149 

1.6009 

80.2 

40.26 

51.45 

Target orbit 

20428 

.7600 

77.00 

40.00 

50.00 

Modified orbit 

20471 

.7563 

77.00 

40.00 

50.00 

Achieved orbit 

19967 

.7513 

76.13 

39. 77 

50.53 


The program computed the "best" modified orbit to be the one in which r 

P 

alone is modified. The insertion velocity required was 886 m/sec. The 
time interval between decision and execution was .707 days. 

After computing the time of the execution event and the velocity correction 
to be made at that time NOMNAL returns to the propagation of the nominal 
trajectory. At the required time, the insertion velocity is added impulsively 
and the resulting conic elements relative to the target planet are computed. 

The resulting conic is described in the previous chart. The achieved and 
modified orbits would have been improved if the decision event had been 
entered later. 

5.1.2 Planetary Explorer Venus '78 Mission 
a. Sample Data 

N800=3,NB=1,3,4, NCPR=2Q0, TMPR=500., NLP=4, NTP=3, ACKT=2.5E-5, 
KALI=1978, 8, 17,0,0, SI=0., IZER0=1, 

T IMG = 0 • , 90., .5, .29, 

K T I M = 1, 1, 2, 3, 

KTYP= 1, 1, 3, -1, 

KMXa= 3, 3, 4, 3, 

MDL = 1, 2, 1, 1, 

LVLS= 2, 1, 

PERV=. 90001, .00005, 

0VMAX=. 01, .01, IBA0S=1, 

KT AR = 7 ,8,3, KTAR< 1 , 2 ) =7 , 8, 4, KTAR<1,3>=1, 

TAR=7000. ,-50. , T AR ( 1 , 2 ) = 750 0 . , -60 . , TAR ( 1 , 3) =27 00 0 . , . 75 ,5 . , 

T0L= 100. ,1.,. 01, TOL(1,2)=5Q.,.50,. 005, 

KALT=1973,12,16,0,0, TS=0., 

KALT (1,2) =197 8, 12, 16,5,0, GS(2)=0. , 

PULrtAG=.001, PULMA S = 1 . , DUR = 1., DTI = .l, 

AC=l.E-4,2.5E-5, AC ( 1 , 2 ) =2 . 5£-5 , 

N0IT=12, 

• 

IZ£RO=0, Zi)AT = 1.2244485426 £ + 8, -8. 913915 32Q5E+7 , — 4. 6158128562E+3, 
6.1928578,22. 7126246,-3. 2751748, 

KALI =197 8, 8, 17, 4, 49,31=15.201, 
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The data defining the Planetary Explorer Venus '78 mission is given above. 
This data will now be explained in detail. 

The first line defines the nominal trajectory propagation between guidance 
events. The sun, Earth, and Venus are the gravitational bodies used in 
the trajectory integration, the Earth acting as launch planet and Venus 
acting as target planet. A moderate accuracy level of 2.5E-5 is used. 
Printouts of trajectory information are given every 200 integration incre- 
ments so that the frequency of output is a function of the nearness to the 
virtual mass. 


The initial date is read in as 8/17/1978. Since IZER0=1 the zero iterate 
will be based on the planet- to-planet option. Thus for the zero iterate the 
initial position will be Earth at the initial time and the final position 
will be Venus at the time given by KALT(5,1), TS(1). 


Four guidance events are to be processed during the trajectory. The four 
events are generally defined columnwise on the input for clarity. Each 
event will be discussed separately. 


The first event is a targeting event (KTYI^l) occurring at 0 days (TIMG=0.) 
after the initial time (KTIM=1). The correction is to be computed and 
executed (KMXQ=3) using an impulsive model (MDL=1) for the execution. The 
target parameters are r 

O 


t, t (KTAR=7,8,3) with desired values 7000 km, 

- 


-50 , 1978/12/16 (TAR and KALT) and tolerances 100 km, l", .01 days .(TOL) 
respectively. The Newton- Raphs on scheme is to be used (METH=0) with the 
perturbation size 10 km/sec (PERV) and a maximum allowable step of 10"^ 
km/sec (DVMAX) during the progressive accuracy levels of 10 - ^, 2.5x10”^. 


After targeting and executing the first guidance event, the trajectory is 
to be integrated at the second targeting event (KTYP=1) occurring 90 days 
(TIMG=90.) after the initial time (KTIM=l). This event has target parameters 
°f r cA’ i ’ £ CA (KTAR=7,8,4) with slightly different target values and 

tolerances. The targeting scheme parameters are identical to the first 
event except that now only one accuracy level is used to do the targeting. 
After determining the correction to be made, the execution is to be done 
(KMXQ=3) using the pulsing arc model (MDL=2). The pulsing arc parameters 
are set as thrust magnitude: .001, nominal mass: 1, thrust duration: 1, and 
time interval between pulses: .1 day (PULMAG , PULMAS , DUR, DTI). This 
determines that 1 m/sec of velocity will be imparted on each pulse. 

The third guidance event is a coplanar orbit insertion event (KTYP=3, KTAR=1) 
to be processed a half-day (TIMG= .5) after encountering the sphere of 
influence of Venus (KTIM=2). The desired orbit is to have a semimajor 
axis of 27000 km, an eccentricity of .75, and a periapsis shift of 5° (TAR). 
After generating the time of execution and the correction to be executed, 
the trajectory is to be integrated to that time and then added impulsively 
(KMXQ=4 , MDL=1) . 
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The final guidance event is a termination event (KTYP = -1) to be performed 
at .25 days (TIMG) after closest approach to Venus (KTIM = 3). After 
integrating the trajectory from the insertion execution to this time, the 
program is ended. 

It should be noted that the data as recorded above represents two successive 
runs. The data to the blank line is the first run. The first run while not 
targeting the first event in the allowable number of iterations did signif- 
icantly improve the zero iterate. Therefore a second run was made in which 
the last three rows of data were added to the original data. This has the 
effect of storing the later values over the earlier values. 

Therefore in the second run IZERO was set to zero to permit the direct input 
of the partially targeted initial position and velocity vectors in ZDAT. The 
initial date KALI, SI was also updated to the exact time of injection rather 
than the launch date. Thus the results of the first run are used to good 
advantage in making the second. 

b. Sample Output 

Selected pages of the actual output from this run are supplied in the Appendix 
to this volume. 

c. Discussion 

The Planetary Explorer Venus '78 mission is launched from Earth on 8/17/1978 
and arrives at Venus on 12/16/1978. The massless planet trajectory generated 
for the zero iterate for these dates has the following properties: 


SMA: 

128. 38M km 

TAL 

185.78° 


ECC: 

.18117 

TAP 

328.54° 


INC: 

2.8123° 

TOF 

121 days 


The hear earth conic based on the 

input launch profile and 

the departure 

asymptote of the 

above described 

heliocentric conic may be 

summarized : 

C3 : 

8.801 

RAD: 

6567.3 


SMA: 

-45,599 

VEL: 

11.410 


ECC: 

1.1439 

INJ 

TIME: 4 hr -49 m -15 

sec 


The injection state was computed from this conic. When integrated to the 
SOI of Venus the radius of closest approach is 592,130 km with a time error 
of 1.43 days. Counting the two runs, twelve iterations at the low level 
and four at the high level are required to reduce the errors to 39 km in 

r (=7000), . 7° in i (=-50°) and .005 days in t . The velocity 
CA CA CA 

correction needed to accomplish this and the elements of the refined near 
earth conic are 

Av = 10.4 m/sec 

SMA = -45593. 

ECC = 1.1439 
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The corrected trajectory is now integrated to the midcourse maneuver ninety 
days later. It is desired to increase the r CA by 500 km, vary the inclin- 
ation by 10°, and delay the arrival by 5 hours. This targeting requires 
five iterations to generate the /\v of 34 m/sec which yields errors of 

Ar = .26 km, Ai = .001°, and At = .000 days. The execution model 
CA CA CA 

to be used in the implementation of this correction is the pulsing arc. 

The velocity increment imparted per pulse is 1 m/sec with pulses occurring 
at intervals of 2.4 hours. Thus the pulsing arc requires 34 pulses and 
lasts 3.3 days. The current state as generated by the virtual mass trajectory 
is propagated backwards over half the pulsing arc interval. Then alternately 
single pulses are added and the resulting state propagated until the final 
pulse has been executed. All propagations in this process use a conic 
propagation corrected by the direct term effects of the launch and target 
bodies. A comparison of the heliocentric conics corresponding to the 
uncorrect^d trajectory, the impulsively corrected trajectory, and the pulsing 
arc trajectory is provided below. 


a 


i 


Uncorrected 
Impulsive 
Pulsing Arc 


128 383 969 
128 343 645 
128 341 692 


.18071 2.87435 
.17996 2.87120 
.17997 2.87122 


The n-body trajectory then picks up at the end of the pulsing arc to integrate 
to the next guidance event. ' 

A coplanar orbit insertion occurs .82 days before closest approach. The 
elements of the approach hyperbola, the target orbit, the modified orbit, 
and the orbit actually achieved upon the later execution are compared below: 


Orbit a 

Approach Hyperbola -12625 
Target Orbit 27000 
Modified Orbit 27916 
Achieved Orbit 28327 


e 

w 

i 

Q 

68285 

-102.50 

63.53 

-161.06 

75 

-97.50 

63.53 

-161.06 

69256 

-97.50 

63.53 

-161.06 

69424 

-97.51 

63.59 

-161.16 


The approach hyperbola and original target ellipse do not intersect: hence 
the target orbit is modified to obtain a tangential intersection by varying 
r^. The time until execution is computed along with the velocity impulse 

to be added at that time (Av = 2080 m/sec). The trajectory mode is then 
reentered and the nominal trajectory is integrated to the time of the execution 
of the insertion. Upon addition of the velocity increment, the orbit achieved 
is given in the last row of the above table. 

The integration of the nominal trajectory continues after the addition of 
the impulsive insertion until the termination event .25 day after closest 
approach . 
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.1.3 Lunar Viking '76 Mission 
a. Sample Data 

N300=2, N3=4,L1, NLP=4, NTP=11, I9ARV=i, ACKT=2.5E-5, NCPR=1Q0, 
I2ER0=LQ , RP =656 3.365, ZD AT=- 72 0 0 . , 1 9 2D . , 8 9 . , l TAR G=l> TMPR=2(3 » , 

<Tyo = lf 3, -i, 

KTIH= 1, 2, 3, 

TIHG= 3., .3, 1., 

KMXQ= 3, 4, 

L VLS= 1 , A C = 2.5E-5, 

SPHFAC=, 25, 

KT AR= 7,8,3, 

7 AR=1 82 0 . , 39. , KALT=1976,6,20,12,0 , TS = 0., 

Tni=25.,l.n, .805, 

KT AR ( 1 , 2 ) = 1, T ARC 1,2) = 3000. ,.4, +5., 

The data listed above for the Lunar Viking '76 mission is typical of the 
data required for lunar targeting. 

The first line defines the nominal trajectory. The earth and the moon 
are the only bodies in the integration, the earth being the launch planet 
and the moon the target body. The inertial coordinate system is to be 
the earth-moon barycentric ecliptic system (IBARY=1). The accuracy level 
of 2.5E-5 is a reasonable figure. Printouts of the nominal trajectory data 
will occur every 100 integration increments. 

The second line defines the zero iterate ocmputation. Lunar targeting is 

specified (IZERO=10) with target conditions at closest approach to the moon 

given as a = -7200, r = 1820, and i = 89°. The input parking orbit 

CA CA 

radius is specified as r = 6563.365 km. 

P 

Three guidance events are specified. A targeting event (KTYP=1) is requested 
at the initial tiiiie (TIMG=0 . , KTIM=1) . The velocity refinement is to be 

both computed and executed (KMXQ=3) . The target parameters are to be r , 

CA 

i , and t (KTAR) with desired values of 1820 km, 89°, and 1976/6/20/12 
CA CA 

(TAR,KALT) and tolerances of 25 km, 1 , and .005 days respectively. For 
this event the SOI of the moon is to be reduced to one- fourth its usual 
value (SPHFAC=.25) so that the conditions at the (reduced) SOI will be 
nearly identical to those at closest approach. The Newton- Raphson scheme 
will be used with the internally stored perturbation size PERV=. 00001 and 
maximum step DVMAX=1. Only one level of accuracy will be used. 

A coplanar insertion event (KTYP=1, KTAR=1) will be processed at .3 day 
(TIMG=.3) after intersecting the (normal) SOI. The desired conic elements 
are a semimajor axis of 3000 km, eccentricity of .4, and a periapsis shift 
of 5 . The impulsive model will be used for its execution at the required 
time. ’ 

The program will be terminated upon reaching the termination event (KTYP=-1) 
one day after lunar closest approach. 


a.3 
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b. Sample Output 


Selected pages of the actual output frju this run are supplied in the 
Appendix to' this volume. 


c. Description 


The Lunar Viking '76 mission discussed here is a 100-hour trajectory to the 
moon, arriving at closest approach to the moon on 1976/6/20/12. Because of 
the relatively long flight time involved, this is a reasonable test of the 
lunar trajectory targeter. 


The first guidance event involves determining the injection time, position, 

and velocity required to yield a trajectory with a radius of closest 

approach of 1820 km, an equatorial inclination of 89°, and the time of 

closest approach defined above. NOMNAL first generates a patched conic 

trajectory meeting the targeted conditions and a semi-major axis relative 

to the moon of -7200. Seven iterations are required to produce a patched 

conic trajectory having errors of A a = 3.4 km, A* = 7 km, Ai = .1 , 

LA LA 

and A t . = 0. The injection conditions of this targeted patched conic are 
LA 

then input to the multi-conic targeter. The zero iterate of this process 

has errors of A a = 1150 km, Ar = 630 km, and At = .1 days. One 

CA CA 


hundred and one iterations are required in the multi-conic targeting to 


reduce these errors to Aa = 1 km, 


A r^ A = 10 km, 


A i. 


= .r 


and 


""CA “ * “““ & fc CA 
.001. The targeted multi-conic trajectory when integrated in the virtual 
mass model has errors of A r^ A = .15 km, A i r4 = .4°, and Atp A = >002 


CA 


CA 


and hence falls within the allowable tolerances. The elements of the targeted 
near earth conic are 


a = 214325 km i = 46.46° 

e = .96941 42= 2.23° 

U) = -166.53° 


This trajectory is then integrated to the time of the second guidance event, 
an orbit insertion decision, occurring at .3 days after SOI intersection. 

The elements of the approach hyperbola, the target orbit, the modified 
orbit and the orbit actually achieved upon the later execution are compared 
below. 


Orbit 

a 

e 

OJ 

i 

Q 

Approach hyperbola 

-7224 

1.25206 

135.50 

90.70 

-138.53 

Target Orbit 

3000 

.40000 

140.50 

90.70 

-138.53 

Modified Orbit 

3008 

.39608 

140.50 

90.70 

-138.53 

Achieved Orbit 

3028 

.39406 

135.57 

88.75 

-138.39 

The discrepancy between 

the modified 

orbit and 

the achieved is caused by 

the fact that the time 

interval over 

which the 

insertion 

maneuver 

was 

predicted was so large 

(.3 days) that the conic 

approximation was 

bad. 

Therefore if the orbit 

insertion is 

required more accurately, the 

time of 

the orbit insertion decision should 

be adjusted 

later. 
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5.2 ERRAN Sample Cases 


Two error analysis cases performed by the error analysis program 
ERRAN will be described in this section to illustrate the operation and 
versatility of ERRAN. The two cases to be discussed are: 

Case E-l. Planetary Explorer Venus '78 Mission 
Case E-2. Lunar Viking '76 Mission 

5.2.1 Planetary Explorer Venus '78 Mission 

a. Sample Data 

X I = 6.55455973324E2» 4. 625391 52599F3.-4.615 8128562E3, 

-lo0856 02 36 161F], - 1. 25168 693609, -3. 27287460527, ICOOR =2. N90 0=3. 
Nl_P = 4 »NTP=3»NR(1 ) = 1 ,3 »4»ACC=?.5E-5»TTNJ=0. » 

TN0M3=1 .366 335904 27 F4 ,6.06308 86 80 ? 1 F 3 , 1 1 9 . 4 8 80 5958 » 

TN0MC = -5. 46 45 381 5 909E3 ,4.5978234725762 .-4. 8246241 1161F3. 
2.21797795136,-9.90977032363.-3.45655168685. 120.79425159, 
LMO=8,LOAY=17,LHR=4,LMIN=49,SECL=15.201.LYR=1978,IMO=12. 

IDA Y=16 , IHR=5.TMTN = 53 .SEC 1=1 P.4 65 , T Y R = 1 9? P , I S TMC= 1 » 

IAUGIN(1»=1. 1,1,1 MJG I N ( 11)=2»1.0,1,0.0.1,0.2. NENT=7.NEV1=1» 
Tl=60. »NFV2=1 ,T2= 100. 0 1 , TP T2= 1 20 . 794 . NE V 3 = 4 , 7 3 = 0 . . 1 0 . . 1 00 . . 
X20.4sP(l»l)=2.»P(2»2)=2.»P(3»3)=2.» D (4»4)=19.E — 6»P(5»5)=18.E-6» 
P(6,68=18.£-6«PS(1»1 > = 1 . »PS(2,2) = 9.E-12,U0(1 »1 > = 100. , 
U0J2 b 2)=1.F— 14, U 0(3,31=82.8 IE— 1 6 . V0 ( 1 » 1 ) = 1 . 7F- 7 , V C ( 2 . 2 ) =2 . 2E -1 5 . 
V0C3.3J =2.5F-14,T0NF=0,MNCN(3)=25.E-6,9.E-12,25.E-6t9.E-12» 
25oE-6»9.E-12,STGRFS-16.E-10,SIG°R0=16.E-6,SIGALP=20.E-5, 
SIGRET=20.F-5 ,PSIGS=4.E-1 0,PSTGK=4.F-6 ,P 5 I G A =1 0 . F-5 . 
PSIGB=10oE-5,IPRTNT=3.KPRINT=l.IPRT=0,0.0»l,PULMAG=.001. 
PULMAS=lo,DUR=l. , 0TI=.1»XTAR(1, 41=27000. ,.75,5. ,PR0BI=1. F-5, 
IGUID(lo2S=2»0,l»l»3,IGUID(l.?)=2fO»l»2#3tIGUID(l»3)=2»0»l*2,3» 
IGUID51 ,41=4, 0,0.1 ,4. 


The first two lines of the above sample data define the dynamic model 
assumed in the error analysis. The spacecraft position/velocity injection 
state (XI) is referenced to geocentric ecliptic coordinates as indicated 
by the value of IC00R. The motion of the spacecraft will be subject to the 
influence of the three (NB0D) celestial bodies indicated by the NB array, 
namely, the Sun, Venus, and the Earth, and will be launched from the Earth 
(NLP) toward Venus (NTP) . A moderate trajectory accuracy (ACC) will be 
employed in the generation of the spacecraft trajectory. Since the initial 
time is the injection time, TINJ is set to 0. Nominal B-plane and closest 
approach target conditions are defined by the TN0MB and TN0MC arrays. 

These arrays, which are obtained from the N0MNAL run which generated the 
injection conditions, are required in this error analysis run because several 
biased aimpoint guidance events are to be performed. The launch and final 
dates are defined by the series of variables LM0 through IYR. State tran- 
sition matrices will be computed analytically using the patched conic tech- 
nique as indicated by the value of ISTMC. 
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The IAUGIN array defines the parameter augmentation for this run, 
and indicates that there are two solve- for parameters: target planet 
gravitational constant bias and range-rate biasj three dynamic consider 
parameters: target planet semi-major axis, inclination, and mean anomaly 
biases; and three measurement consider parameters: station 1 radius, 
latitude, and longitude biases. The number of entries in the measurement 
schedule, which is presented below, is specified by NENT. 


The variables NEV1 and T1 indicate that an eigenvector event will 
occur at 60 days, while NEV2, T2, and TPT2 indicate that a prediction 
event will occur at 100.01 days and will predict to 120.794 days. Four 
guidance events are scheduled at 0., 10., 100., and 120.4 days in this 
run, as indicated by NEV3 and T3. The characteristics of these guidance 
events are specified by the IGUID arrays which appear in the last two lines 
of the previous sample data. IGUID(1,1) indicates that the first guidance 
event will be 2VBP, linear, subject to planetary quarantine constraints., 
impulsive, and both computed and executed. The second and third guidance 
events differ from the first in that they will employ a pulsing thrust 
model. The final guidance event is. an impulsive planar orbital insertion, 
not subject to quarantine constraints, computed, and executed at the 
appropriate time. The thrust characteristics of the pulsing engine are 
defined by the variables PULMAG, PUIA1AS, DUR, and DTI. The variable PR0BI 
indicates that the probability of impact with the target planet Venus 
must not be more than l.xlO - "’. The XTAR array indicates that we desire to 
insert the spacecraft in an orbit about Venus having a semi-major axis of 
27000 km., an eccentricity of .75, and a periapsis shift of 5 degrees. 


The spacecraft injection covariance matrix is assumed to be diagonal 
with position variances of 2 km7 and velocity variances 18 m. /sec. . The 
parameter covariance matrices are also assumed to be diagonal. Solve- for 
parameter variances are defined in the PS array; dynamic consider parameter 
variances, in the UO array; and measurement consider parameter variances, 
in the VO array. The arrangement of the elements in the PS, UO, and VO 
arrays must correspond to the structure of the solve-for, dynamic consider, 
and measurement consider parameter vectors. 


\ 


Dynamic noise will be absent from this run, as indicated by the value 
of IDNF. Noise corrupting range and range- rate measurements from all three 
tracking stations will be assumed to have the statistics described by the 
MNCN array. Execution error statistics are described by the variables SIGRES 
through PSIGB. 


Variables not defined in the previous sample data take on internally- 
specified values. For example, the internally-specified values of all required 
numerical differencing factors (DEIHUP, DELAXS, DELICL, and DELMA) were 
deemed satisfactory for this run, and so were not defined in the sample data. 

The measurement schedule input for this error analysis run is presented 
below : 
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.2 

9.2 

1 . 

3 

.6 

9.6 

1 . 

5 

.9 

9.9 

1 . 

7 

19.5 

99.5 

5. 

4 

100.3 

120.3 

2. 

3 

70.1 

90.1 

10. 

9 

110.4 

114.4 

2. 

10 


The first row indicates that range- rate measurements from station 1 
will be taken once a day beginning at .2 days and ending at 9.2 days. The 
other rows are interpreted in a similar fashion. 

b. Sample Output 

Selected pages of the actual output from this run are supplied in the 
Appendix to this volume. 


c. Discussion 

Planetary quarantine constraints require that the nominal aimpoint 
be biased at injection since the probability of impact (POI) exceeds the 
allowable POI. The linear theory indicates that a bias AV of nearly 
15 m/sec. is required to reduce the POI to l.xlO . The bias in B*T is 
computed as 8674 km. , the bias in B*R, as 119,447 km. At 10 days we re- 
compute the target conditions and find that they do not agree with the 
desired bias aimpoint. This indicates that nonlinear guidance should have 
been employed at injection to re- target the trajectory to the bias aimpoint. 
The linear theory, however, does provide us with valid bias AV's at guidance 
events occurring in the heliocentric phase. After the midcourse at 10. days 
the spacecraft is once again heading toward the nominal aimpoint since 
planetary quarantine constraints are no longer violated. 

Shortly before encountering the sphere of influence of Venus, space- 
craft position uncertainties (la) have values of 5.6, 6.1, and 12.6 km. in 
the x,y, and z directions, respectively. Velocity uncertainties have values 

of 3.3x10“ 6, 4. 5x10“ 6, and 6. 7x10“ ^ km. /sec. At this time the uncertainty 
(1(7 ) in the station 1 range- rate bias has-been reduced from 3.xl0“^ to 

.62x10“^ km. /sec. 

At the orbital insertion decision event at 120.4 days we compute an 
insertion Av of 1.96 km. /sec. to be executed at 120.77 days. After execution 
we are close to the desired Venus orbit. 
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5.2.2 Lunar Viking '76 Mission 
a. Sample Data 


X I=-6. 21904 171 F3*-l. 992 9084 8F 3. -6. 5§4 35808E2« 3.162899754* 

-9.6011 53 294 . -4. 17 327 24 04® IC 00 R = 2* NR 00=2 ♦ NLP =4 t NT Prl 1 * IBARY-1* 

ACC = 2.5E-S»NR = 4. 1 1 . LM0 = 6 *L 0 & Y = 1 6 »L HP = 1 1 « L MI N = 7 . SE CL = 50 . 055 » 
LYR=1976*IM0=6. IDAY=20. IHR=i9 e IMIN=30sSECI=43«022.IYR=1976® 
FACP=.01 »FACV=l.E-6* IAU6XNf7)=l»Iol®0»2»2*0®09l9l»0»2»0»l * 
NENT=7*NEV2=1 *T2 = 1.01 . TP T2 = 2 . 99 . NE V 3 r 3 . 1 3= 1 . *3. * 3.83* 
IGUID<l*ll=3*0*0*l*3*IGUTD(l*2)=l*P*0»l*3*I6UIDfl»3)=4*0»0»l»4* 
XTARIlt3>=3000.*.4»5.»P(4»4>=9.E-6»P<5»5)=9.f-6*P«6»6)=9.E-6* 
PSfl*l»-3.6E-3»PSC2*2> = 100<,»PS(3»3)=9.E-F»U0fl*l) = l.E-12» 

U0(2»2 > = 1 .E-12»V0< 1 .1 ) = . 169E-6.V0f 2 *2 i-.223E-14*VOI3*3)=2.46E-14® 
V0«4*4 > = 1 ,E-8*STGPR0= 25 . E- 6 * S XGRE 5 = 9 . E -1 0 . 51 6 & LP= 5 .E -5 * 
SIGBET=5.£~5*IPRINT=3eKPRINT=lfIPRT=0«0*0*l» ISTMC=3* 
TN0MB=-2.36489848784E2»-5.47059458514E3»3 e 433385* 

TNOMC -5 . 2 53 575073E2 * ~1 o 1 1 1 622 873E 3 • 1 • 359994681 E 3 * 1 .05 330 1 945 » 
-1.498851835*-!. 63200 4962*4. 037704*NOftCCrl*ACCNO=;l.F -4* 


The first two lines of the above sample data define the dynamic model 
assumed in the error analysis. The spacecraft position/velocity injection 
state (XI) is referenced to geocentric ecliptic coordinates (IC00R) . Only 
two (NB0D) bodies the Earth and the Moon, will govern the motion of the 
spacecraft. The IBARY code indicates that the barycentric ecliptic coordin- 
ate system will be used as an inertial reference coordinate system. A 
moderate trajectory accuracy (ACC) will be used to generate the spacecraft 
trajectory. 

The ISTMC code indicates state transition matrices will always be 
computed using numerical differencing. Numerical differencing factors 
FA'CP and FACV different from the preset values were selected since the 
preset values are more suitable for interplanetary trajectories than for 
lunar trajectories. The last two variables indicate that a reduced 
trajectory accuracy will be employed in numerical differencing. 

The IAUGIN array defines the parameter augmentation for this run, and 
indicates that there are three solve-for parameters: gravitational constant 
and semi-major axis biases of the Moon, and range bias; two dynamic consider 
parameters: longitude of the ascending node and argument of periaps*is biases 
of the Moon; and four measurement consider parameters: station 3 radius, 
latitude, and longitude biases and star-planet angle 1 bias. 

A single prediction event occurs at 1.01 days as indicated by NEV2 
and T2. The prediction is made to 2.99 days (TPT2) . The variables NEV3 
and T3 define a schedule of three guidance events occurring at 1., 3., and 
3.83 days. According to the IGUID(1,1) array, the first guidance event is 
3VBP, linear, not subject to planetary quarantine constraints, impulsive, 
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and is to be both computed and executed. The second guidance event is 
identical to the first except that it is FTA rather than 3VBP. The third 
guidance event is a planar orbital insertion, applied impulsively, and not 
subject to quarantine constraints. It is to be computed and executed 
later at the appropriate time. The desired lunar orbit has a semi-major 
axis of 3000. km., an eccentricity of .4, and a periapsis shift of 5 degrees. 
This target orbit is defined by the XTAR array. 


The spacecraft injection covariance matrix is assumed to be diagonal 
with position variances equal to the pre-set values of 1 km.^ and velocity 
variances of 9 m. /sec.^, as indicated by the P array. Solve-for parameter 
variances are defined in the PS array; dynamic consider parameter variances, 
in the UO array; and measurement consider parameter variances, in the VO 
array. The arrangement of the elements in these latter three arrays conforms 
to the structure of the three parameter vectors. 


Measurement noise variances will take on internally-specified values 
and so do not appear in the previous sample data. Execution error variances 
are defined by the variables SIGPR0 through SIGBET. 

b. Sample Output 

Selected pages of the actual output from this run are supplied in the 
Appendix to this volume. 


c. Discussion 


Prior to the first guidance event at 1. days from injection, spacecraft 
position knowledge uncertainties (la) have values of 1.85x10”^, 6.36x10“^, 
and 4.31x10“^ km. in the x, y, and z directions, respectively. Velocity 
uncertainties have values of 1.55x10 , 1.24x10”^, and 1.53x10”^ km. /sec. 

The control covariance at this time indicates position control uncertainties 
(la ) of 1075., 579., and 279. km. and velocity control uncertainties of 
1.86x10” , 5.44x10" , and 2.87x10“^ km. /sec. These control uncertainties 
are a measure of the dispersion of the actual trajectory from the nominal 
and are useful in determining the efficacy of the guidance (control) process. 
Propagating these control uncertainties forward to the target indicates 
that if no guidance correction is applied target condition dispersions will 
have la values of 8106 km. in B*T , 1889 km. in B*R , and .057 days in 


SI 


With the execution of the guidance correction, we can reduce target 


condition dispersions to Iff values of 48 km. in B*T , 39 km. in B*R , 
and 3.5x10” days in t . The expected value of the velocity correction 

at 1. days is 22.2 m./sec. and was computed using the Hof fman- Young formula. 


The orbital insertion decision event at 3.83 days computes an orbital 
insertion AV of .525 km. /sec. to be executed at 4.038 days. After execution 
of the insertion AV the spacecraft is in a lunar orbit close to the desired 
orbit . 
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5.3 SIMUL Sample Cases 


Two simulation cases performed by the simulation program SIMUL will 
be described in this section to illustrate the operation and versatility 
of SIMUL. These cases correspond to the two ERRAN sample cases. The two ' 
cases to be discussed are: 

Case S-l. Planetary Explorer Venus '78 Mission 
Case S-2. Lunar Viking '76 Mission 

5.3.1 Planetary Explorer Venus '78 Mission 

a. Sample Data 

X Ir6 0 55 45 597 3 324F?«>4c,62 57915Z599E3.-4.61581285 62E3» 

-1=0856 02 36161E1 . - 1 . 2 51 8869 36 09 « -3 . 27 28 7 4 605 27 . 1 C OORr 2 » NB 0 0= 3 . 
NLP = 4 «NTP=3»N9 tl )rl»3»4»aCC = 2.5E-5?>TTNjrO.» 

T M0MB rl = 366 3359 04 27E4„6«,0630 S T868021E 3. 119=48 80 595 8 » 

TW0WCr-5 = 46453@15909E3» 4.59 78234 7 257E 2 9-4.8246 241 1161 E3» 

2 = 2179779 5136 o -9 o 90 97 70 32 363. -3. 4 5655 1686 8 5. 120.794 2515 9* 
LM0r8 s L0AYrl7,LHPr4.LMINr49»SECL-15.201.LYPr 1978*IM0r 12* 
IDfiYrl6»IHRr5oYMTNr53»SECTrl0.465*IYRrl978»IST4C-l» 

I a UGINU ) = 1 • 1 « lfia UGTN(ll>-2. 1.0. 1.0* 0.1*0. 2. NENT = 7 .NE Vl= 1. 
Tl=60..NEV2=l.T2=10r,Gl . TPT 2r 1 20 . 794 .N E V3 = 4 , T3r 0 . . 10 . . 100 . . 
120o4*Pil.l)r?o»P?2*2)-2. . P ( 3 . 3 ) - ? . . P ( 4 . 4 ) = 1 8 . E-6 . P ( 5 . 5 > = 1 8. E- 6. 
Pd6o6lrl8.E-6.PS! 1. 11=1.. °S<2. 2)r9.E-12.UC(l.l) = l 00.. 

UC8292» = l.E-14,lJ0(3.3)r82.81E-16.V0( 1 . 1) - 1 . 7E- 7 . V 0 ( 2 . 2 > - 2 . 2E -1 5 . 
V053.3)-2.5F-l4»TPMFr0»MNCN(3lr25.F-6.9.F-12»25.r-6»9«E-12. 

25 = E-6.3oF-12.SI5c>F6rl8.E-ia.SIGPP0rl6.E-F.ST9ALP = 20.E-5. 
SIGP.ET=20.E-5.PSTGS = 4 » E - 1 0 » PS TGK = 4 . F -& .PS I G A = 10.F-5. 

PSIGB-10.F-5 . IPRJNTrJ ^KPRTNTr 1 t IPRTrO.C.O. 1 .PULMAG-, 001 » 

PULM&S=l.»0UP-l..nTT=ol.XTaR(l,41=27000...75.5.»PR0RI=l.E-5» 

a0EVXrl„®- o 5.1o.-2=F-3.2oE-3»-1.5F-3»ACCl-2.5F-5.NB0Dlr3» 

NR 1 o 3o4 oOMUPPr .8.08 0=15. .n T 9rl.5F-7,DMflRr-12.F-8.RI&f4)z6.E-6. 
SL4=8.E-4.-in.E-8 .3«E-7» IAMNPrl .aV»R“ C 3 ) r 2 5 . r- 8 . 9 . E- 14 . 25 . E- 8 . 
9<>E — 14925.E— 8*9»E~14t>4*2o5E-l 1 , A 9FS = . 0 1 , -2 . E -5 . -2 . F- 5 * -3 .E -5 » 
fl PRO - « 0 1 » 2 . E - 3 » 2. E _ 3 . 4. E -3 * A AL P- . 0 1 » 1 . E-2 . 1 . F-2. 1 • E-2 • 
aBETr B Qlo-l.E-2»-loF-2®-loE-29NEV5r2»T5=4.»50.»NDACC-l. 

ACCND = 1 .E-4 . TGUTn U .1 5 r2. 1 .2. 1 . 3 * TGUI0 < 1 . 2 }- 2 , 1 , 2 . 2. 3 ♦ 

IGU ID ( 1 o3)r2* 1 »2»2*3. TGUI0( 1 » 4 ) r4 » 1.0*1. 4»XT0L(1 .1 ) = 100. .100*. 
o00 5 . XT0L ( 1 . 2) rlOP. . 1 CO. ..005. XTOL < l .3 1 = 100. .100. . -005. 

XAC U .1 ! = 2.5E-5. XAT( 1 . 2 ) = 2 . 5E-5 . X A 0 ( 1 . 3 } =2 . 5F- 5 . XPER V = 3 * 1 . E- 5 . 
XOVMflXrls.nS.LLVLSrZ*!, 


The first half of this sample data has already been discussed in the ERRAN 
Sample Cases section. The data peculiar to the simulation program is concerned 
primarily with the definition of actual dynamics and actual error statistics, 
and begins with the ADEVX array above. 
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The ADEVX array defines the actual initial spacecraft position/ 
velocity deviation from the nominal trajectory. Actual trajectory accuracy 
(ACC1) and actual celestial bodies (NB0D1 and NB1 array) are not different 
from the assumed values (ACC, NB0D, and NB) , although normally the actual 
dynamic model would be more refined than the assumed model. Some differ- 
ences between these two models have been defined: actual dynamic biases in 
the target planet gravitational constant, semi-major axis, inclination, and 
mean anomaly will be added to the nominal values in the generation of the 
actual trajectory. These dynamic biases are specified by the variables 
DMUBP, DAB, DIB, and DMAB. 


Actual station location errors have been defined in the SLB array 
for station 1. Range-rate measurements from this station will be biased 
by 6.xKT^ m./sec. as indicated by the variable BIA(4). Actual measurement 
noise variances are different from their assumed values as indicated by 
the value of IAMNF. The actual variances are defined in the AVARM array. 
Actual execution errors corresponding to the four guidance events are 
defined in the ARES, APR0, AALP, and ABET arrays. The first value in each 
of these arrays is a dummy value for this particular run since the first 
guidance event occurs at injection. The object of a guidance event at 
injection is solely to change the injection velocity state; injection 
statistics are assumed to remain unchanged. 


The characteristics of the first 3 guidance events differ from the 

first 3 guidance events in the corresponding ERRAN sample case in that 

nonlinear guidance will be used to re-compute the guidance correction. 

This necessitates the specification of targeting tolerances for these three 

events in the three XT0L arrays. In this case all three XT0L arrays are 

identical and indicate that the re-targeted trajectories have tolerances 

of 100 km. in B*T and B*R , and .005 days in t . A single accuracy 

o X. 

level, which is usually chosen to be identical to the trajectory accuracy 
ACC, is employed in nonlinear guidance as indicated by the variables LLVLS 
and XAC. The XPERV array indicates that a velocity perturbation of l.xl0“^ 
km. /sec. will be used to compute the targeting matrix at each nonlinear 
guidance event. Velocity steps of .05 km. /sec. in the targeting process 
are permissible, as indicated by the XDVMAX array. 

b. Sample Output 

Selected pages of the actual output from this run are supplied in the 
Appendix to this, volume. 


c. Discussion 

Nonlinear guidance is employed at each midcourse guidance event in the 
case under consideration. At 0. days the aimpoint is biased in order to 
satisfy planetary quarantine constraints and the bias velocity is re-computed 
using nonlinear guidance. The nonlinear bias velocity has a magnitude over 
twice as large as the linear bias velocity, which indicates that linear 
guidance is not valid during the initial phase of the trajectory. 
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The measurement data presented at 9.9 days shows a convergent navigation 
process at that point in the trajectory. All orbit estimation errors fall 
within their predicted ( 1 cr ) standard deviations „ 

At the guidance event at 10. days the nominal aimpoint satisfies the 
planetary quarantine constraints. The velocity required to remove the 
bias (DVRB) has a magnitude of about 15 m./sec. Linear guidance is valid 
in this region since the sum of the commanded correction (to null out target 
errors) and the velocity correction required to remove the aimpoint bias 
agrees quite well with the total velocity correction computed using non- 
linear guidance. The z-components of the linear and nonlinear AV's will 
always differ for 2VBP guidance since the z-component has been constrained 
to be zero for nonlinear 2VBP guidance. 

The guidance correction at 10. days is executed using the impulse series 
thrust model operating over a 1.3 day arc. The re-computed target conditions 
at 100. days show, as one might expect, that the impulsive guidance policy 
will not satisfy the target conditions if an impulse series is used to execute 
the impulsive AV. 

5.3.2 Lunar Viking '76 Mission 

a. Sample Data 

X I=— 6. 21904171E3.-1.99290848E3.— 6. 55435808E2.3. 162899754. 
-9.601153294.-tt.l7327 24 04»ICO0R-2»NBOO=2*NLP: < *,NTPrll,IBflRYrl, 
aCC=2.5E-5,NB=4,ll.LM0r6,LDAYrl6.LHR=ll.LMIN=7.SECL=50.C55» 
LYR=197S.IM0=S»inAY=20,IHR=19.IMlN=30,SECI=43.022.TYR=1976» 
FACPr.Ql.FftCV-lcE-S® IAUGIN(71::1. 1.1. 0.2.2. 0.0. 1,1. 0.2. 0.1, 
NENT=7»NEV2=1 »T2=1.01 .TPT2=2. Q 9»NFV3=3»T3=1. » 3 . i3.83i 
I GU ID<1. ll=3»lor, 1, 3.T GUIOfl. 21=1. 1.0.1. 3. TGUIDCl. 31=4.1.0.1. 4. 
XTAR < 1 , 3 > = 30Q0<». .4 »5. ,P( 4 ,4 1 = 9.E-6,P (5 .5 1=9„F-6,P (6.61 =9oF-6 » 
PS(l.ll=3«6E-3»PS(2»21=inn.»PSC3.3)=9.E-6*UC(l,l)=l.E— 12. 

U0C2 .2 1 = 1.£- 1? »VP< 1.1 )-.169E-G»Vn(2o2)-o2?7F-14oVO(3*3l-2.46E~14 * 
VOI 4»4)=l.E-a«STGPP0=25.E-6»ST6RES=9.E-10*SIGALP=5.E-5* 
SIGBET-5eE-5»IPRTNT-3.KPRINT-l»XPRTrn?0909 1 » ISTMC-3. 
A0EVX=.3»-.3«.3*>1.5E-3» 1.5E-3.-1 . 5E - 3 , AC C 1 =? . 5E-5 » 
N90Dl-2.NRl-4.il. 1 

DMUP8=.n4,DAB-5..nN08=5.E-6.0WB=~S.F-6»8IA(3l=~5.F-3. 
BI«(9)r6.E-5.SLB(7)=9,E-4.-7.F-8.2.5F-7.IAMNFrl, 
AVARM{3)=25. r -8»9.E-14»25.E-8»9.E-14»25.F-8»9.E-14»4*2.5E-ll , 
ARES=-2.E-5 . -2.F~5,-3.E-5.APR0=2.E-3.2.E~3»4.E-3. 

A ALPrl„F-2. 1 «F-2. 1 . E -2 . A BF T r 3*- 1 .E-2»NEV5r4, T5r,2. .7, 1 .5.2. . 

TN0MB :-2. 364 99848 784E 2. -5. 470 594585 14E 3. 3.033385. 

TN0MC -5.253 575073F2.-1. 11 1622 873F3.1 .35999468 IE 3. 1 .053301945 . 

-1. 498851 8 3 5. -1.63 2004962. 4. 9 377 04, NO ACCrl , ACCN0- 1 .F -4 , 

XT0L (1.11=100. .100. ♦. 005 .XT0L ( 1 • 2 ) = 1 00 . , 1 00 . . 1 00 . , 

X AC ( 1 » 1 ) =2.5E-5 . XAC( 1 .21 =2 . 5F-5 ♦ XPERV = 2* . 000 1 . 

XDVMAX = 2*.010,LLVLS=2*1 , 
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The first half of this sample data has already been discussed in 
the ERRAN Sample Cases section. The data peculiar to the simulation 
program is concerned primarily with tin- definition of actual dynamics 
and actual error statistics, and begins with the ADEVX array above. 

The ADEVX array defines the actual initial spacecraft position/vel- 
ocity deviation from the nominal trajectory. Actual trajectory accuracy 
(ACC1) and actual celestial bodies (NB0D1 and NB1 array) were chosen to he 
identical to the assumed values (ACC, NB0D, and NB) . However, some differ- 
ences between the actual and assumed dynamic models were defined: actual 
dynamic biases in the gravitational constant, semi-major axis, longitude 
of the ascending node, and argument of periapsis of the Moon are to be 
added to the nominal values to generate the actual trajectory. These 
dynamic biases are specified by the variables DMUPB, DAB, DN0B, and DWB. 

Actual station location errors have been defined in the SLB array 
for station 3. Range measurements from station 1 will be biased by -5 
meters, while star-planet angle 1 measurements will be biased by 6.xlO~ J 
radians. These biases are defined in the BIA array. Actual measurement 
noise variances are different from their assumed values as indicated by the 
value of IAMNF. The actual variances are defined in the AVARM array. 

Actual execution errors corresponding to the three guidance events are 
defined in the ARES, APR0, AALP, and ABET arrays. 

The characteristics of the first 2 guidance events differ from the 
first 2 guidance events in the corresponding ERRAN sample case in that 
nonlinear guidance will be used to re-compute the guidance correction. 

This requires the specification of several targeting variables for these 
two events. At the first event tolerances of 100 km. in B‘T and B*R, 
and .005 days in t are imposed; at the second event tolerances of 100 
km. in each of the final position components, are imposed. These tolerances 
are defined in the XT0L arrays. A single accuracy level (LLVLS) is usually 
selected for a nonlinear guidance event. The accuracy level itself (XAC) 
is usually set to the trajectory accuracy (ACC). The XPERV array defines 
the velocity perturbation to be used in the computation of the targeting 
matrix, while the XDVMAX array defines the maximum permissible velocity 
changes that can be used in the targeting process. 

b. Sample Output 

The output from this run is provided in the accompanying document, 
Tabulated Runs. In the final documentation selected pages from that output 
will be reduced to standard size and included in the User's Manual. 

c. Discussion 

All midcourse guidance velocity corrections have been recomputed using 
nonlinear techniques. Comparison of the linear and nonlinear AV ' s at the 
1st guidance event indicates very good agreement. No nonlinear AV is 
computed at the guidance event at 10. days since the target tolerances of 
100. km. in each position component are satisfied by the uncorrected 
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trajectory. This is further substantiated by the target condition 
standard deviations just prior to the guidance event at 3 days which 
indicate standard deviations of 24, 31, and 43 km. in the x, y, and z 
target position components, respectively. 

Comparing actual orbit estimation errors with predicted 1 a position 
and velocity uncertainties indicates that the navigation process is generally 
convergent. In a convergent process the actual errors should be bounded by 
the ± 3 a predicted uncertainties. Prior to encountering the sphere of 
influence, of the Moon at 3.3 days our error in the range bias estimate has 
been reduced from 3 meters to .16 meters. Errors in the estimation of the 
gravitational constant and semi-major axis biases of the Moon, however, 
have not been reduced. This may be due to the strong correlation between 
these two parameters in the trans- lunar phase of the mission, making it 
difficult for the estimation process to separate the effects of biases in 
these two parameters. 


\ 
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5.4 Multiprobe Sample Cases 


5.4.1 NOMNAL Multi-Probe Sample Case 

The exact data as input for a representative 1977 Planetary 
Explorer targeting problem is given below: 

C VIRTUAL-MASS INTEGRATOR DATA 
ACKT=2 . 5E-05 , 

NB0D=3 , 

NB=1 ,3,4, 

NLP=4 , 

NTP=3 , . 

TMPR=50 . , 

NCPR=500 , 

C INITIAL TRAJECTORY TIME 
KALI=19 77,1,4,6,49 
SI=38. , 

C GUIDANCE EVENT SCHEDULE 
KTYP=1 ,4, 1,5, 1,-1, 

KTIM=1, 1,1, 1,1,3, 

TIMG=0 . ,122., 122. 5, 123., 123. 5,1., 

C TARGET SPECIFICATION 
KTAR=13 ,14,1,0,0,0, 

0,0, 0,0, 0,0, 

13.14.1.0. 0.0, 

11 . 11 . 0 . 0 . 0 . 0 , 

13,14,15,0,0,0, 

TAR=0 . ,68.2,0. ,0. ,0. ,0. , 

0. ,0. ,0. ,0. ,0. ,0. , 

-14. ,100. ,0. ,0. ,0. ,0. , 
30.3,-14.6,-34.4,71.3,137.7,66.6, 

20. ,30. ,0. ,0. ,0. ,0. , 

T0L=1. ,1. , .001,0. ,0. ,0. , 

0. ,0. ,0. ,0. ,0. ,0. , 

1. ,1. , .001,0. ,0. ,0. , 

1 . ,1 . , 1. ,1 . ,0. ,0. , 

1. ,1. , .001,0. ,0. ,0. , 

KALT=1977 , 5,17,6,49, 

0,0, 0,0,0, 

1977.5.17.6.49, 

0,0, 0,0,0, 

1977.5.17.5.49, 

TS=38. ,0. ,38. ,0. ,38., 
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C INITIAL INTERATE DATA 
IZER0=1, 

C TARGETING SCHEME DATA 
LVLS-6*'2 f 

AC=1 .0E-04 ,2 . 5E-05 ,0.0, 0.0, 0.0, 

0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 

1 .0E-04 , 2 . 5E-05 ,0 .0 ,0 . 0,0.0, 

0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 

1 . 0E-04 , 2 . 5E-05 , 0 . 0 , 0. 0 ,0 . 0 , 

PERV=6*1 .OE-05 , 

C0NTR=6*2 . , 

DVMAX=6*5 .0E-01, 

N0IT=6*12 , 

MAXB=6*12, 

WGHTM=1 . 0E+05 ,0 . , 1 .0E+05 , .5.1.0E+05, 

C PROBE TARGETING DATA 
IPCS=6*1, 

RPS=6*6200. , 

A brief description of the nominal mission profile is given to 
motivate the selection of guidance events . A detailed explanation 
of the data is then provided. Finally a discussion of the target- 
ing output is given. 

a. Mission Description 

The trajectory is to be of the Type I category with a launch 
data of 1/4/77 and an arrival date of 5/17/77. The entire con- 
glomerate vehicle consisting of the bus, the main probe, and the 
three miniprobes is targeted at injection to impact the main probe 
target site at 0° declination and 68.2° right ascension in the 
planetocentric subsolar frame. The injection-targeted trajectory 
is then flown uneventfully until 11 days prior to entry. At this 
time the main probe is released to impact its target site without 
any further velocity correction. Twelve hours after the main 
probe release the conglomerate spacecraft now consisting of the 
bus and miniprobes is retargeted to a pseudo-impact site at -14 
and 100° declination and right ascension respectively. Since the 
bus as well as the miniprobes will subsequently be retargeted, 
this site will never be reached by any portion of the spacecraft. 
It serves the purpose instead of shifting the ballistic trajectory 
of the spacecraft so that at the time of miniprobe release its 
point of contact with the p.lnnet is centered among those of the 
miniprobes. Ten days before entry the miniprobes are deployed by 



a simultaneous release from the spinning spacecraft. A sufficLent 
number of release controls is not available to permit exact tar- 
geting of all three miniprobes to their respective target sites. 
Hence controls are used which minimize a miss index of the impact 
site distribution; they are tabulated in the Discussion section. 
Finally, nine and one-half days before entry, the bus is retargeted 
to its desired impact site at 20° declination and 30° right as- 
cension. The bus is accelerated at this time to impact the planet 
one hour before the probes thereby easing the data management 
burden . 

b . Sample Data 

Consider the sample input of Table 1. First observe the input 
to the virtual-mass propagator. ACKT sets the VMP accuracy level 
for propagation of the nominal trajectory between guidance events 
at the moderate accuracy level of 2.5E-5. NBOD indicates that 
3 gravitational bodies are to be considered in the virtual-mass 
integration. The array NB identifies these as the sun, Earth, 
and Venus. NLP and NTP further identify the Earth as the launch 
planet and Venus as the target planet. TMPR triggers trajectory 
status printouts every 50 days while NCPR initiates them every 
500 integration steps. 

Examine next the scheduling of guidance event times. The 
array KTYP states that there are to be 6 guidance events in all. 

The first, third and fifth are ordinary targetings; the second 
is a main probe propagation; the fourth is a miniprobe targeting, 
and the sixth is a termination event. The array KTIM specifies 
that the first five event times in the array TIMG are to be re- 
ferenced to the initial trajectory time at injection while the 
entry for the last event is referenced to closest approach of the 
target planet. Observing that the trajectory is 132 days long 
the reader should satisfy himself that the event times listed in 
TIMG will indeed produce the nominal trajectory described in 
the mission profile. 

Consider next the specification of the targets. The array 
KTAR specifies that all of the targeting events have the time, 
right ascension and declination at impact as their target variables. 
The first and third events, however, permit extrapolation of these 
target values from the integrated state at the SOI, while the 
fifth requires virtual-mass propagation all the way to impact. 

For the fourth or miniprobe targeting event KTAR dictates that 
the spin-axis orientation is to be of mode 1 (both the ecliptic 
right ascension and declination of the spin axes are free release 
controls) and that for targeting purposes the miniprobes are to 
be propagated according to the conic model . The TAR array contains 
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the right ascensions and declinations of the respective target 
sites. The reader can easily verify that they agree with those 
given in Table 1 for the intended mission. The array TOL specifies 
that for all the targeting events the desired right ascensions and 
declinations must be achieved to within 1 deg. and the desired 
impact times to within 0.001 day. For the miniprobe targeting 
event, TOL states that equal unity weighting factors should be 
applied to the B.T and B.R errors for all the miniprobe target 
sites and that the weighted sum of the change in length of the 
release control vector and the change in magnitude of the miss 
index be less than 1 for convergence in the least-squares routine. 
The target times are input through the variables KALT and TS in 
accordance with the nominal mission profile discussed above, 
assuming the trajectory initial time given in KALT and SI. 

Actually, the initial date need only be specified to a calendar 
day; the hours, minutes, and seconds at injection are computed in 
the zero iterate computation using the internally set launch , 
profile with a Cape Kennedy launch. 

Next consider the zero iterate data. IZERO specifies for an 
initial iterate, a Lambert massless-planet conic from the launch 
planet at the initial time to the target planet at the target 
time. 

Now study the targeting scheme data. The array LVLS specifies 
that all of the targeting events are to have two levels. Accord- 
ing to AC the first level should be propagated at a VMP accuracy 
level of 1.0E-4 and the second at 2.5E-5. The velocity perturba- 
tion size used in approximating the sensitivity matrices of the 
target variables to velocity controls is given by PERV to be 
1.0E-5 throughout the run. The launch— planetocentric velocity 
controls are to be used exclusively according to CONTR. The max- 
imum permissible velocity correction is universally fized at 0.5 
km/sec by DVMAX. The maximum number of iterations and bad steps 
are both set wherever applicable to 12 by the variables NOIT and 
MAXB, respectively. The weighting factor of timing errors to dis- 
tance errors for bad step calculations is assigned the value 
1.0E5 by WGHTM for all targeting events. For the miniprobe tar- 
geting event the length of pseudo-inverse release control correc- 
tion is bounded above by 0.5 according to WGHTM. 

Consider the special probe targeting data. The array IPCS 
sets the plantocentric probe-sphere frame to subsolar orbit-plane 
coordinates while RPS sets the radius of the probe sphere to 6200 
km for all probe related targeting. 



c. Sample. Output 


Selected pages of the actual output of this run are supplied 
in the appendix to this volume . 

d. Discussion 

The interpretation of the output from the NOMNAL programs is 
aptly illustrated in other sample cases. Only the printout from 
the 'miniprobe targeting algorithm and the main probe propagator 
are unique to this case. This "new" output is thoroughly described 
in the Output Description . 


Table 1 Sample 1977 Planetary Explorer Targeting Results 


Body 

Declination of Impact 
Site in Degrees 

Right Ascension of 
Impact Site in 
Degrees 

Julian Date Epoch 
1900 of Impact Time 

Desired 

Achieved 

Desired 

Achieved 

Desired 

Achieved 

Main 

Probe 

0.0 

0.3 

68.2 

71.1 

28260.784 

28260.78460 

Bus 

(Prior to 
Miniprobe 
Release) 

-14.0 

-12.90 

100.0 

102.4 

28260.784 

28260.78450 

Miniprobe 

1 

30.3 

30.9 

71.3 

74.1 


28260.78731 

Miniprobe 

2 

-14.6 

-11.9 

137.7 

139.8 


28260.78422 

Miniprobe 

3 

-34.4 

-34.5 

66 . 6 

70.8 


28260.79070 

Bus 

(Final) 

20.0 

20.0 ' 

30.0 

30.0 

28260.743 

28260.74292 
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Table 2 Comparison of Results of Conic and 
Virtual-Mass Miniprobe Propagation 
Models 


Miniprobe 

Declination of Impact 
Site in Degrees 

Right Ascension of 
Impact Site in 
Degrees 

Julian Date Epoch 
1900 of Epoch Time 

Conic 

N-Body 

Conic 

N-Body 

Conic 

N-Body 

1 

-12.279 

-11.932 

139.718 

139.798 

28260.78397 

28260.78422 

2 

-34.649 

-34.451 

70.944 

70.785 

28260.79087 

28260 .79070 

3 

31.622 

30.922 

73.562 

74.054 

28260.78746 

28260. 7873L 


The results of the targeting as summarized in Tables 1 and 2 
deserve some comment. First, the disparity between the desired 
and achieved impact sites for the main probe and the bus prior 
to miniprobe release as compared to the bus after final retarget- 
ing is caused by the use of target option 1 for the former cases 
and 15 for the latter cases. Thus, in the first two instances, 
the trajectory was integrated to the SOI and then conically extra- 
polated to impact, while in the latter cases, it was integrated 
over its entirety. Second, the promixity of the various miniprobes 
to their respective sites is more than satisfactory. The initial 
control estimate was relatively accurate and the least-squares 
iterations proceeded entirely by pseudo-inverse steps. The sym- 
metrical distribution of the miniprobe target sites about the 
bus pseudo-impact point no doubt facilitated the miss minimization. 
Table 2 compares the respective impact sites and times of the 
miniprobes propagated under both the conic and virtual-mass 
models. The respective impact sites will be observed to agree 
to within a degree and the times to within 0.025 days illus- 
trating the accuracy of the conic miniprobe propagation model. 
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5.4.2 ERRAN Multiprobe Sample Case 


a. Input data 


The input data for this sample case consists of the namelist 
ERRAN and three measurement schedules. The namelist defines 
the mission and the filter design and consists of the follow- 
ing cards: 

XI=-l,03633510620E+8, -6. 16 23486 7839E+07 , -1 . 27614089726E+06 , 
2.15882811727E+01, -2.62615614273E+01, 1.04039769246E-01, 
LM0=4 , LDAY=14„ LHR=5 , LMIN=22 , SECL=19.879, LYR»1977, 

IM0=5, IDAY=16, IHR=23, IMIN=54, SECI=40.788, IRY=1977, 
IC00R=O, 

NLP=4, NTP=3, NB=1 ,3,4, 

P(1,1)=1.E4, 

P(2,2)=1.E4, 

P(3,3)=1.E4, 

P(4 8 4)=2.5E-3» 

P(5 ,5) =2 . 5E-3, 

P(6,6)-2.5E-3, 

MNCN(4)=1. 315E-14, MNCN(5) =0 . 5E-14 , MNCN(8) =0 . 561E-14 , 
IAUGIN=9*1, 

NEV1=8, 

Tl=21. 383, 23. 383, 25. 383, 27. 383, 28. 383, 29. 383 ,30. 383, 30. 883, 
NEV2=1 , 

T2=19 . 383 , 

TPT2=21. 383 , 

ACC=2 . 5E-5 , 

DTMAX= . 5 , 

KPRINT=1, 

IPRINT-10, 

V0(1,1)=. 1495568E-5, 

V0( 1, 2) =-. 166533E-9 , 

V0(2 ,1) =-. 166533E-9 , 

V0(2 , 2) = . 1854361E-13 , 

VO(3,3)=.3328075E-12, 

V0(4, 4) =. 130 420 7E -5 , 

V0( 4, 5) =- . 1741507E-9 , 

V0(5 , 4)=-. 174150 7E-9 , 

V0(5 , 5) = . 2325433E-13, 

VO(6,6)=.3817565E-12, 

V0( 7 , 7) = . 1498273E-5 , 

VO(7»8)=. 16641E-9, 

VO(8,7)=. 16641E-9, 

VO(8,8)=. 1848279E-13, 

V0(9 ,9) = . 3323088E-12 , 
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NENT=9 , 

NEV3=1, T3=18.5 , IGUID(1 , 1)=2 ,0 ,0 , 1 ,3 , 

SIGPR0=2 . 25E-4, SIGRES=1.0E-8, 

SIGALP=6 . 85E-4 , SIGBET=6 . 85E-4 . 

T6=19 . , 

T7=20. , 

SMN(4)=1. 315E-10, 

SMN(6)=0.5 E-10, 

SMN(8)=0.561E-10, 

ISMN=1 , 

IUTC=0 , 

WFLS=3*1. , 

DCTP=2*0. ,-60. , 

RATP=165. ,105. ,135. , 

S0=O . 5 , 

ISA0=1, 

IPCSK=2 , 

YYL=1. 5E-3 , 

TIMPCT=35.0, 

RPS=6200. , 

XEE(1)=0.01, 

XEE(2)=2.E-12, 

XEE(3)=4.E-4, 

XEE(4)=4.E-4, 

XEE(5)=4.E-4, 

NENT1=9 , 

NENT2=9 . 

The first section of namelist ERRAN variables, beginning with 
XI and ending with SIGBET, defines the bus or primary vehicle 
Venus approach trajectory, the bus event schedule, and the navi- 
gation filter design. The initial bus state at the beginning 
of the approach phase is defined by XI and IC00R. The initial 
bus position/velocity covariance is defined by the P array. 
Doppler measurement noise variance are specified by the MNCN 
variables. The navigation filter solves for no dynamic or 
measurement parameters, but considers the nine station loca- 
tion biases as indicated by the IAUGIN vector. The covariance 
matrix for these consider parameters is given by the VO array. 
Eight eigenvector events and one prediction event are scheduled 
for the bus. The NEV3, T3, and IGUID variables indicate that 
a linear 2VBP midcourse guidance event will occur at 18.5 days. 
Since the variable IGEN does not appear in namelist ERRAN, a 
standard error analysis will be performed for the bus and all 
four probes. A generalized covariance analysis could be per- 
formed for the bus, although the four probes currently can only 
be treated in the standard error analysis mode. 
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The final section of namelist ERRAN variables, beginning the 
T6 and ending with NENT2, defines the main probe and miniprobe 
release events. The variable T6 indicates that the miniprobe 
will be released at 19. days, shortly after the final bus mid- 
course correction, while all three miniprobes will be released 
at T7=20. days. A different set of doppler measurement noise 
variances will be used for miniprobe tracking, as indicated 
by variable 1SMN. These new variances are defined by the SMN 
variables. The absence of the IPMN and PMN variables indicate 
that the measurement noise variances used for the bus will also 
be used for the main probe. The miniprobe release controls 
have not been specified by the user in this example, as is in- 
dicated by IUTC. Thus, the miniprobes must be targeted using 
the ERRAN program. Variables WFLS through IPCSK are required ’ 
to perform this targeting. The miniprobe release execution er- 
ror variances are specified by the XEE. 


The measurement 

schedules for the 

bus , the 

main probe, and the 

three miniprobes 

are each defined by nine 

cards according to 

variables NENT, 
cards are listed 

NENT1, and NENT2 . 

. • 

The bus 

measurement schedule 

0.0167 

32.7000 

1.00 

3 

0.0958 

32.7790 

1.00 

3 

0.1750 

32.7580 

1.00 

3 

0.1760 

32.7590 

1.00 

7 

0.3500 

32.7330 

1.00 

7 

0.5450 

32.7280 

1.00 

7 

0.5710 

32.7540 

1.00 

5 

0.7870 

32.7740 

1.00 

5 

0.9880 

32.7710 

1.00 

5 

The main probe measurement schedule cards 

are listed next. 

19.0167 

32.8000 

2.00 

3 

19.0958 

32.8790 

2.00 

3 

19.1750 

32.8580 

2.00 

3 

19.1760 

32.8590 

2.00 

7 

19.3500 

32.8330 

2.00 

7 

19.5450 

32.8280 

2.00 

7 

19.5710 

32.8540 

2.00 

5 

19.7870 

32.8740 

2.00 

5 

19.9880 

32.8710 

2.00 

5 
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Finally, the measurement schedule cards used for all three 
miniprobes are listed. 


20.0167 

30.4000 

1.00 

3 

20.0958 

30.4790 

1.00 

3 

20.1750 

30.5580 

1.00 

3 

20.1760 

30.5590 

1.00 

7 

20.3500 

30.7330 

1.00 

7 

20.5450 

30.9280 

1.00 

7 

20.5710 

30.9540 

1.00 

5 

20.7870 

31.1740 

1.00 

5 

20.9880 

31.3710 

1.00 

5 


b. Output discussion 

Selected pages from the output of this sample case appear in 
the appendix, where it is referred to as case MP-2 . Only por- 
tions of miniprobe 1 output are presented since the output 
formats for the main probe and the remaining miniprobes are 
essentially the same. Output associated with the targeting 
of the miniprobes occurs on page 4. Targeting was successful, 
as indicated by KKWIT = 0. The final target controls are 
given by the UCNTRL vector. Immediately following this infor- 
mation, the release execution error covariance matrix for mini- 
probe 1 is given. The control correlation matrix partitions 
and standard deviations at entry time 32.781 days for miniprobe 
1 are obtained by adding this execution error covariance matrix 
to the bus covariance matrix immediately before releasing the 
probe and progating the result to entry. The matrix represent 
the l-o dispersions of the miniprobe deviations about the nomi- 
nal target site. The probe state relative to Venus is written 
out in plane to centric ecliptic and subsolar orbital-plane co- 
ordinates before transforming this state and the previous con- 
trol covariance matrix to entry parameter coordinates. The 
following page shows the output for measurement 10 of mini- 
probe 1. The format is identical to bus measurement output 
and requires no further explanation. The output for this sample 
case was generated on the IBM 360 computer at GSFC. 
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5.5 Generalized Covariance Analysis Sample Cases 

Two generalized covariance analysis sample cases generated by 
ERRAN will be described in this section. The nominal trajectory 
for both cases is an approach trajectory to Venus beginning at 
about 30 days before encounter. The spacecraft is tracked from 
three earth-based stations using doppler measurements only. In 
each case the objective is to examine the sensitivity of the navi- 
gation filter to off-design conditions. 

5.5.1 Spectral Mismatch Sample Case 

a. Input data 

The input data for this sample case consists of two namelists 
and a measurement schedule. The first namelist, entitled 
ERRAN, defines the mission and the filter design and consists 
of the following cards: 

XI=-1.00905168154E8,-6 . 40 1660 755 81E 7,-1. 72 4087669 17E6 , 
22.4153491439,-25.9718887728,0.248319225059, 
LM0=4,LDAY=15 ,LHR-14,LMIN=42 .SECL-48. 82 ,LYR=1977 , 

IM0=5 , IDAY=16 , IHR=23, IMIN=54 ,SECI=40 . 788 ,IYR=1977 , 

IC00R=O , 

NLP=4,NTP=3,NB=1,3, 4, 

P( 1, 1) =1.E4, 

P(2,2)=1.E4, 

P(3, 3) =1.E4, 

P(4,4)=2.5E-3, 

P(5,5)=2.5E-3, 

P(6,6)=2.5E-3, 

MNCN94) =1. 315E-14,MNCN(6) =0. 5E-14,MNCN(8)=0 . 561E-14 , 
IAUGIN=9*1, 

NEV1=8, 

Tl-21, 383, 23. 383, 25. 383, 27. 383, 28. 383, 29. 383, 30. 383 ,30. 883, 
NEV2=1, 

T2=19 . 383, 

TPT2=21 S 383, 

ACC=2.5E-5, 

DTMAX=1. , 

KPRINT=1 , 

IPRINT=10, 

V0( 1,1)=. 149556 8E-5, 

V0(1 ,2) =-. 166533E-9 , 

V0(2 . 1) =-. 166533E-9 , 

VO(2,2)=. 185436 IE-13, 

V0(3, 3) = . 3328075E-12 , 
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VO(4,4)=. 130420 7E-5, 
V0(4 , 5) =- . 1741507E-9 , 
V0(5,4)=-. 1741507E-9, 
V0(5 , 5) =. 2325433E-13 , 
VO(6 , 6) =. 3817565E-12 , 
V0( 7 , 7) =. 1498273E-5, 
VO(7,8)=.16641E-9, 
V0(8, 7)=. 16641E-9, 
VO(8, 8) =. 18482 79E-13, 
V0(9,9)=. 3323088E-12, 
NENT=9 , 

IGEN=1, 


The spacecraft state at the beginning of the approach phase is 
defined by XI relative to the heliocentric ecliptic coordinate 
system as indicated by IC00R. The computation of the nominal 
trajectory will include the gravitational effects of the sun, 
Venus, and the earth. The initial assumed position/velocity 
covariance matrix is defined by the P-array. The assumed meas- 
urement noise variances for doppler measurements for each of 
the three tracking stations are given by the MNCN variables. 

The navigation filter design involves the nine consider param- 
eters indicated by the IAUGIN vector. These consider parameters 
are the station location biases for all three stations. The 
assumed covariance matrix for these consider parameters is 
specified by the VO array. The variable IGEN set to 1 indi- 
cates that a generalized covariance analysis will be performed. 
The measurement schedule is defined by nine measurement cards, 
as indicated by NENT. These measurement cards are listed. 


.0167 

30.4000 

1.00 

3 

.0958 

30.4790 

1.00 

3 

.1750 

30.5580 

1.00 

3 

.1760 

30.5590 

1.00 

7 

.3500 

30.7330 

1.00 

7 

.5450 

30.9280 

1.00 

7 

.5710 

30.9540 

1.00 

5 

.7870 

31.1700 

1.00 

5 

.9880 

31.3710 

1.00 

5 


The second namelist, GENRAL, which follows the above measure- 
ment schedule cards, defines the "actual" statistics of rele- 
vant parameters. The namelist GENRAL cards for this sample case 
consist of the following. 
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GV(1,1) = . 13460 HE -4, 
GV(1, 2) =- . 1498797E-8 , 
GV ( 2 , 1) =- .1498797E-8 , 



N 


GV(2,2)=.1668925E-12, 

GV(3,3)=.2995268E-11, 

GV(4, 4) = . 1173786E-4 , 

GV(4,5)=-. 1567356E-8, 

GV(5 , 4) =-. 1567356E-8, 

GV(5,5)=. 209289E-12 , 

GV(6,6)=.3435808E-11, 

GV( 7 , 7)=. 1348446E-4 , 

GV(7,8)=. 149769E-8, 

GV(8,7)=. 149769E-8, 

GV( 8, 8) =. 1663451E-12 , 

GV(9,9)=. 2990779E-11. 

For the sample case under consideration, the only difference 
between actual and assumed (by filter) error statistics occurs 
for the station location biases considered by the filter. This 
is indicated by the appearance of only the GV array in name- 
list GENRAL. Actual standard deviations for the station loca- 
tion biases are defined to be three times as large as the cor- 
responding standard deviations assumed by the filter to describe 
these error sources. Cases involving differences in the ac- 
tual and assumed statistical distributions of an error source 
acknowledged (i.e., solved-for or considered) by the filter 
can be referred to as "spectral mismatch" cases for convenience. 

b. Output discussion 

Selected pages from the output of this sample case appear in 
the appendix, where it is referred to as case G-l. The output 
corresponding to measurements 220 and 270 is presented. The 
quantities of interest are the assumed and actual position and 
velocity standard deviations after the measurement is processed. 
At measurement 220, which occurs several days before the space- 
craft pierces the sphere of influence (SOI) of Venus, actual 
standard deviations are about two to three times as large as 
the assumed standard deviations. However, at measurement 270, 
which occurs about a day after the spacecraft has entered the 
SOI, the actual standard deviations range from one to two times 
the magnitude of the assumed standard deviations. These re- 
sults indicate that the navigation filter design is less sen- 
sitive to spectral mismatch for (considered) station location 
biases when the spacecraft is tracked inside the SOI. The out- 
put for this sample case was generated on the CDC 6400/6500 
computer at Martin Marietta's Denver Division. 
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5.5.2 Ignore Parameter Case 
a. Input data 

The mission and filter design have not changed from the pre- 
vious case discussed. Consequently, the namelist ERRAN , and 
measurement schedule cards for both cases are essentially the 
same and need not be reproduced in their entirety. Namelist 
ERRAN contains one additional card, namely. 

IAUGIN(19)=3, 

which indicates that the doppler bias for station 1 will be 
treated as an ignore parameter in the generalized covariance 
analysis . 

The actual variance for this single ignore parameter is de- 
fined by 


GW(1,1)=1.E-10, 


which is the only card appearing in namelist GENRAL for the 
case under discussion. 

b. Output discussion 

Selected pages from the output of this sample case appear in 
the appendix, where it is referred to as case G-2. Output 
for measurement 90 and an eigenvector event occurring at 30.383 
days is shown. As in the previous sample case, we are pri- 
marily interested in comparing actual and assumed position 
and velocity uncertainties. At measurement 90, which occurs 
some distance from the Venusian sphere of influence (SOI) , ac- 
tual standard deviations are nearly an order of magnitude larger 
than the assumed standard deviations. Inside the SOI, however, 
the ratio of actual to assumed standard deviations has been 
much decreased, as is evident in the output for the eigenvector 
event occurring at 30.383 days. This behavior appears plausi- 
ble. Prior to penetrating the SOI the spacecraft velocity is 
not changing rapidly and as a result (doppler) observability 
is reduced. Consequently, ignoring the doppler bias in the 
filter design during this phase can be detrimental. But after 
the SOI has been penetrated, the spacecraft velocity begins to 
change rapidly both in magnitude and direction so (dbppler) 
observatibility increases. In this situation, neglecting the 
doppler bias in the filter design is of less consequence. The 
output for this sample case was generated on the CDC 6400/6500 
computer at Martin Marietta's Denver Division. 
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APPENDIX: SELECTED SAMPLE CASE OUTPUT 



Case N-l. Broken Plane Viking Mars '75 Mission 
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GUIDANCE EVENT AT 121.013 DAYS 

CALENDAR DATE 1978 12 16 5 8 33.046 JULIAN DATE 28838.71427 
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OYNAMIC CONSIDER PARAMETERS 

NODE 1.5259197986EHJ1 -B.4172886457E+00 4. 5567341519E-01 3.4230278412E-03 -1.8979989437E-03 1. 2204699464E-04 

OMEGA 1.5391054792E+01 -8.4734037883E*00 6. 1180794213E-01 3. 4583697505E-03 -1. 9094248684E-03 1. 5620655896E-04 



ORRELATION MATRIX PARTITIONS AND STANDARD DEVIATIONS AT EVENT TIME 1.000 DAYS 
PROPAGATED FORWARD FROM TIME .900 DAYS 
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SENSITIVITY MATRIX TARGETING MATRIX AUX ERROR VEL COR DES AUX VAL DES TAR VAL TAR TOL 

04E+07 -2.31E+05 0. -2.76E-08 -1.56E-09 0. -1.26E*09 3.27E-04 22337.371 22337.371 100.000 
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CORRELATION MATRIX PARTITIONS AND STANDARD DEVIATIONS AT TIME 9.900 DAYS, JUST AFTER THE MEASUREMENT 
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VIRTUAL MASS PROGRAM MILL INTEGRATE UNTIL REACHING ANORMAL STOPPING CONDITION 



DYNAMIC CONSTANT BIASES TO BE USED IN THE DETERMINATION OF THE ACTUAL TRAJECTORY 
GRAVITATIONAL CONSTANT OF THE SUN. ...... 0. 

GRAVITATIONAL CONSTANT OF TARGET PLANET. . 4.0000000000000 E- 02 
SEMI-MAJOR AXIS OF TARGET PLANET . „ o . . 5 . 0000 0 000 000 00 E+ 0 0 
ECCENTRICITY OF TARGET PLANET . 0. 
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OYNANIC CONSIDER PARAMETERS 

NODE 1.88487 42 106E+01 -7. 4240 865992E* 00 7. 1410730340E-01 3 .4086853366E- 0 3 -1.9048 147881E- 03 1. 2068663794E-04 

OMEGA 1. 973491435 3E +01 -7. 2643147963E+00 9. 2878042182E-01 3. 4384529843E-0 3 -1.9162264598E- 03 l. 5561400787E-04 



CORRELATION MATRIX PARTITIONS AND STANDARD DEVIATIONS AT EVENT TIME 1.000 DATS 
PROPAGATED FORWARD FROM TIME .900 DAYS 
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CONTROL CORRELATION MATRIX PARTITIONS ANO STANDARD DEVIATIONS JUST BEFORE GUIDANCE CORRECTION AT TIME 1.000 DAYS 
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MEASUREMENT CONSIDER PARAMETERS 



NON-LINE AR-GJI DANCE 
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COMP. Y - COMP. 2 - COMP. RESULTANT 
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VIRTUAL MASS CATA 

VIRTUAL MASS POSITION -0.153430305750 08 -0 . I 0759 434573D 09 - 0 . 60 7564 32 e I 60 06 0.108684501910 
VIRTUAL MASS VCLJCITY 0 . 3 44 2 1 2 1 0 06 3D 02 -0 . 5 1 692 835499D 01 -0.205671640840 01 0 . 3486791 1 853D 
SPACECRAFT f Ob. REL. TU V.M. -0.103556638140 05 O . 1 075966553 ID 04 - 0 • 2577 5 30683 2D 04 0.107 25723 35 40 
SPACECRAFT VfcL. REL • TO V.M. 0.6922 9050 94 00 01 0 .3 1 722526029D 01 0 . 47625608693D 01 0 . 898 1 74 74664D 
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STAT6= 58648144. <57 72546. -1054235. 13.49123 9.291 11 0.53727 

I0r>T= 4 CCJN= 0 .46956 7 €7930- 03 0. 20762o64 1 1 D-04 -0 . 23760 9627 l D-05 

DV= 0.19227043710-03 0 . 5 48 J 0 684 6 1 0- J3 0.14397552230-04 



SENSITIVITY MATRIX TARGETING MATRIX AUX ERrtUK VLL COR OES AUX VAL DES TAR VAL TAR TOL 
2<+l> 05 -0»1 20 03 0. 1 30 07 -O. 190-05 0.430-06 -0.440 00 -0.250 03 0.190-03 -15053.652 20.000 1.000 
4 U> 06 -O.oSD 06 -0.100 05 -0. 930-07 -0.30D-03 -0.290-02 0.210 03 0.55D-03 -2436.790 30.000 1.000 
190 01 0.5 IL 02 -0.15D 02 -U. 750-07 -0.61D-07 -0.19D-01 0.190-03 0. 140-04 26260.743 28260.743 0.001 
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case HP-2. EMU 8 Multiprobe Sample Case 



MiMl PKUHL f* VENT AT TRAJLCTORY TIME 20.000 DAYS 


— — O C 

— — — ~ o 

O I I I I 

O O O O Q 

O Nlfl ftl h 
• O IS O h 

O N — N <f 

Al O AJ O O' 

o O' <* n 

— — n in in 

n h « in s 

■> ** rn n n 

<M — <0 <M 
« N O >0 


O Q 

4 ® 

- O' 

in 4 

© O' 
C\J O' 

- O' 
© O' 
O O' 

- O' 

- O' 


O O O O O o 


UJ * 
J * 

o * 

□ * 


© n 

M - n - '] c o 

o o o o o 

O G 

o o o o o r •© 

© rc o -* w © 4 

r t m ^ o ,n 

CON o © © AJ 

© n © n n m o 

f J 1 ? O .} 4 S if 

H in (t N O N 

M 4 M O N ‘C N 

© — © c o* o © 

lo ifl ") o in w O' 

N N O "J - © 

• . • • • o o 


— a> 

V © 
> 4 


' 0 1 — N - © 

> o - c c © 

III I 

3 0 0 0 0 0 

> O © O 4 O 

) — © 4 © O 

I N « O' O ifi 

J — iv © o © 

© O' N O <VJ 

O'NfOo — 

io - r o in 

CM © © O O 

AJ — 4 O — 

- N O - - 


o o o o c o 


O 0* — O N O' 

— c — o v o 

III II 
Q O O O O O 

n o n © 4 © 

aj i© © © n © 

4 aj n © o o 

O AJ o O' © O' 

n©oO'Nn 

o © — O' © © 

4 © — O' I© N 

O' 4 — O' ON 

•n Ai <M 0 ' 4 <M 

AJ - — O © 

0 o o o o o 

1 


O' © o 

— o o o 

O I I 

o o c o 

o O' © © 
• 4 © 0* 
Off s 4 
ai © r O' 

— — 7 
*- -4 — O' 
M C* © O' 
N O O' 
i\J — O' 


3 © 4 

0 o o 

1 I 

O Q O 
O — © 
N O' © 
■”m © o- 
<c m r> 
o 4 o 
o © n 
o^o 

© ai n 

-TO 

f\i n — 


o o o o o o 


i I i 

O 3 O C 
0< © AJ © 
■ O' AJ AI N 

o o m 

n — o n 

— — ON 


i© 


1 «n 


© a 
o a i m — 

— m o © 

— \j — r* 


0 0.0 
ON© 
— At © 

n © aj 

V © U' 
4 N — 
« "1 4 
U O N O 
%. a» m ii 
“ ill 4 yi 
N AJ N 4 

in n m 

n n 


< in 


c 


m aj © 
n 4 n 
"3 — AJ 


■X O O c o O O 

O I I I I 

I- X > NX > N 


< D 
I • Z Z 
: u z o u 

J D © > 


->-11 


' fr i 

c o 


T S <r - o O 


i- u < i i < < 

uj a x x 

u. O U l/> to 

O IU - - z z 

< * ■ 1 * x o © 

< < - — 


Q- a uj z x i 


■ j) < < < < 


o r 'u a. lj 


-J z v. z z 


- Z Z _J J 

N Cl Aj O G 

Y N k- U U 

4 Z z o u 

U < 4 c u 


r> N C 4 4 C 
— O O O O O O 


o o o r r n o 

c © © n .,© e © 


— o © <? © -*> n © 

t— 'vj js N N 7- © •'> 

— *r -r - ?• © n 


c 4 *>i o a> — tn 

c © © — o: aj i 


o o o o o 


Oj © © © 0/ © 

o O' O' o c c 


7- O O' © ?• O' 




A -124 


CJRPbLAT ItJN MATRIX PARTITIONS AND STANDARD DEVIATIONS AT CVENT TIME 20.000 DAY 
OR0MAGATLD FOR WARD FRO 4 T Mb i DAYS 
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TARGETEO NOMINAL TRAJECTORY ENCOUNTERED SPHERE OF INFLUENCE AT TRAJeCTORY TIME 32.77735 OAYS 
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X Y Z RESULTANT 

POSITION RELATIVE TO TARGET PLANET -0.2745748636D 04 -0*16375747170 04 -0 • 5 3 1 2 1 7 5964D 04 0.62000000000 04 

VELOCITY RELATIVE TO TARGET PLANET 0 .664 26091 66D 01 0.51478853640 01 0 • 7 3 1 43 1 78 1 40 01 0.11141158880 02 
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RADIUS 2 0.05577962 0.05706960 0.06061731 0.06122813 0.05820295 0.05925948 
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(LTk) COVARIANCE MATRIX 



RANGE-RATE WAS MEASURED PROM STATION 1 AT TRAJECTORY TIME 21.01670 DAYS 
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DIAGONAL OF DYNAMIC NU I Sri MATRIX 
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JUST BEFORE THE MEASUREMENT 
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RANGE-RATE WAS NEASURED FROM STATION 3 AT TRAJECTORY TIME 24.17600 OATS 
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OBSERVATION MATRIX PARTITIONS — TRANSPOSES SHOWN 

RANGE -RATEC3 ) 
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GAIN MATRIX PARTITIONS 
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MEASUREMENT CONSIDER PARAMETERS 

RADIUS 1 -.06780882 -.12545936 .10620438 -.01225188 -.02120595 .13615625 
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RANGE-RATE WAS MEASURED FROM STATION 2 AT TRAJECTORY TIME 29.98900 DAYS 
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diagonal of dynamic noise matrix 
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gain matrix partitions 



65904522 52E-01 
2454433 046E-01 
2354494762E+3Q 


OMO ° 
ONfJ N 
«=» K O' 
O^MQ 

a ffl O' a 
C3HO » 

o a >c n 


e 5(J'>ONi 
rs eo U> <U ® 
OWN W J 

o & cr< oo 

O 95 O' o O' 
Ofjtf NO 
aro JO>»c 

o O' O' *H s. 


a O' O N fO 
e> O' CO CM UN 
o K « O' W 

o Q O' 
o pg O' w CM 
O o O O' D 
° O' » H W 
0^1150 id 


O N s 9 ) O W 

0^li5WHS 
O rl N. N O' IJ5 
a © fO O W O 
X o ro 4 o vD o 
t=> \D UN f* <M vD 
O N k W kO « 

O CVJ ro If! 0O M 


a in © ^ tO j- 
o o O' h fo S 
a WniftK ^ 

c? U) eo ro 

X O HU5N«1) o' 
oMftOKK 
O O M O' O' N 
oO'CKO'N. 


♦ ♦ ♦ • I I 

Ul U) U) UJ UJ Ul 
© © iC O' O' r< 
^ BWOlOM 
UN « vDS r( « 
>W HNNO^O 
tucocMvor-O'j- 
o O fO W <0 K? UN 

00 '«O(MON 


n n a iO 

□ o o o o o 

♦ * + • I I 

Ul Ul UJ UJ UJ Ul 

o n kj os 

to CT \0 m oo oo 

lAWO'r^NW 
> ^ « CT S ^ M 
uj © a iftn s o' 
Q © Lft >o n N s 
-0-CTU5s.ro <m 
Q Ift o (\j O' O' O' 

M © ^ \C M N O' 


A -143 



DYNAMIC CONSIOER PARAMETERS 



A -144 


3469642 6E + 01 1.000 00000 

17746652E*00 .72878617 1.00000000 

45542539E+02 -.75186776 -.98897175 1.00000000 

41605733E-06 -.79719975 -.99131587 .19563020 1.00000000 

3984589 2E- 06 .37158926 -.34117213 . 33139369 .25749002 1.00000000 
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MEASUREMENT CONSIOER PARAMETERS 

RADIUS 1 -•20575754 -.08294901 ,08487222 ,18912295 -.21399302 -.10044967 
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RANGE-RATE WAS MEASURED FROM STATION 2 AT TRAJECTORY TIME 9.98800 DAYS 
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DIAGONAL OF DYNAMIC NOISE MATRIX 
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2.12579382E*02 -.98B41725 1.00000000 

6.0 17354 2 5E* 0 2 .94134184 -.98005086 1.00000000 

1. 83715121E-05 .98841140 -.99847980 .97226523 1.00000000 

7.6 7717 7 23E- 05 .99658092 -.99715249 .96452503 .99607252 1.00000000 

5.490341 31E-0 4 .97204716 -.99640630 .99125767 .99438851 .98742483 1.00000000 
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MEASUREMENT CONSIOER PARAMETERS 

RS0IUS 1 .0171655*. -.01119957 -.00427859 .01947249 .01313470 .00989542 
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ERROR ANALYSIS MODE-EIGENVECTOR EVENT AT TRAJECTORY TIME 30.333 DAYS 
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VELOCITY 1 EIGENVALUES SQUARE ROOTS OF EIGENVALUES 

1 3. 78816760 26 622E- 14 1 1.9463215569E-07 

2 3 *48 5466 79411 97 E-12 2 1 . 8669404902E-06 

3 1.151417 8 118116E-10 3 1. 0730413840E-05 



VELOCITY EIGENVECTORS 

1 9.7698769G15582E-Q1 2. 1 094560 Q121 05E-O1 -3, 1575365057773E-02 

2 -2.1168983005613E-01 9. 40 81758220 625E-01 -2. 64669 02512 827E-0 1 

3 -2.61241098457&6E-02 2. e526256318155E-01 9, 638222364415QE-0 1 
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IGNORE PARAMETERS 
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